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Abstract: 


During  firing,  propellant  residues  are  scattered  onto  the  soil  surface  where  their  energetic  com¬ 
pounds  can  be  dissolved  by  precipitation.  The  residues,  like  the  unfired  propellants,  are  com¬ 
posed  of  nitrocellulose  imbibed  with  either  2,4-DNT  (single-base),  nitroglycerin  (NG)  (double¬ 
base)  or  NG  and  nitroguanidine  (NQ)  (triple-base).  Although  nitrocellulose  is  insoluble,  2,4- 
DNT,  NG,  and  NQ  are  soluble;  and  2,4-DNT  and  NG  are  also  toxic.  Consequently,  data  on  how 
quickly  2,4-DNT,  NG,  and  NQ  are  dissolved  from  propellant  residues  are  needed  to  determine 
the  flux  of  these  compounds  to  soil.  Once  in  soil  solution,  the  partition  coefficient,  Kd,  and  deg¬ 
radation  rate,  k  values  are  needed  to  predict  the  transport  of  energetics  through  the  vadose  zone 
and  to  groundwater. 

We  measured  the  2,4-DNT,  NG,  and  NQ  dissolution  rates  for  different  propellants  using  labora¬ 
tory  batch  and  drip  tests  where  no  soil  was  present  and  soil  column  studies,  which  used  similar 
propellant  and  residues  as  source  terms,  to  determine  partition  coefficients  and  degradation 
rates.  Because  the  surfaces  of  propellants  and  residues  may  play  an  important  role  in  dissolution 
of  the  energetic  constituents,  we  studied  these  using  both  light  and  electron  microscopy. 

We  found  that  2,4-DNT  is  well  bound  to  NC  and  dissolves  out  slowly,  but  that  both  NG  and 
NQ  have  fast  initial  dissolution  followed  by  slower  mass  loss.  The  amount  of  NG  dissolved  is  a 
function  of  the  NG/NC  ratio  in  the  propellant  and  both  our  mass  loss  data  and  our  microscopy 
results  suggest  that  NG  exists  as  fine  liquid  droplets  within  an  NC  matrix  rather  than  as  dis¬ 
persed  molecules.  NG  droplets  near  the  grain  surface  are  quickly  dissolved  and  once  this  layer 
of  liquid  NG  is  depleted,  NG  diffuses  through  the  NC  matrix  slowly  (~10“14  cm2  s_1).  The  NQ 
also  dissolves  rapidly  initially  but  quickly  mass  loss  for  the  NQ  becomes  smaller  than  that  for 
NG,  despite  higher  NQ  concentrations  in  the  studied  triple-base  propellants.  NQ  is  added  as  a 
crystal  during  manufacturing,  and  was  observed  to  remain  solid  in  the  propellant,  so  dissolution 
of  the  NQ  crystal  would  have  to  precede  its  removal  by  water.  Both  2,4-DNT  and  NG  are  added 
as  liquids  and  cannot  be  distinguished  from  the  NC  matrix.  Therefore,  their  distribution  and 
movement  within  the  nitrocellulose  matrix  is  poorly  understood,  hampering  our  ability  to  derive 
a  physically  based  dissolution  model  that  can  predict  energetic  losses  from  a  variety  of  propel¬ 
lant  types. 

Different  interactions  between  2,4-DNT,  NG,  NQ,  and  the  soils  were  seen  in  both  the  soil  batch 
and  column  studies.  The  2,4-DNT  interacted  strongly  with  soils  and  had  the  highest  adsorption 
and  transformation  rates  measured.  As  a  result,  no  2,-4-DNT  was  detected  in  column  outflow. 
NG  experienced  both  adsorption  and  transformation  in  the  soils,  resulting  in  retardation  of  the 
breakthrough  curve  and  decreased  concentrations  in  the  outflow.  The  short  half-life  of  NG  in 
most  soils  suggests  that  it  should  rarely  reach  groundwater.  NQ,  on  the  other  hand,  does  not 
readily  adsorb  to  soil,  and  does  not  degrade  or  transform.  We  would  expect  that  NQ  dissolved 
from  propellants  would  reach  groundwater. 
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1  Background 

Propellants  are  chemical  mixtures  formulated  to  burn  at  a  controlled  rate, 
to  contain  the  oxygen  needed  to  burn  and  to  evolve  gas.  It  is  the  expanding 
gas  that  propels  the  bullet  or  the  round.  Modern  propellants  are  generally 
composed  of  a  polymer,  a  plasticizer,  and  a  stabilizer.  These  three  compo¬ 
nents  provide  the  structure,  contribute  the  oxygen  and  fuel,  and  slow  the 
deterioration  of  the  propellant,  respectively.  Most  common  are  the  nitro¬ 
cellulose-based  propellants.  About  12  million  kg  (26  million  lb)  of  nitrocel¬ 
lulose  is  manufactured  each  year  in  the  U.S.  (Richie  2012;  a  partial  listing 
is  given  in  Table  At).  This  wood-like  compound  is  used  as  the  structure  of 
the  propellant.  Propellants  are  nitrocellulose  (NC)  impregnated  with  ei¬ 
ther  2,4-dinitrotoluene  (2,4-DNT  single-base),  nitroglycerin  (NG  double¬ 
base),  or  NG  and  nitroguanidine  (triple-base).  Although  nitrocellulose 
does  not  dissolve,  2,4-DNT,  NG,  and  NQ  are  soluble  in  water  (Rosenblatt 
et  al.  1991;  Windholz  1976).  Furthermore,  screening  levels  for  2,4-DNT 
and  NG  are  low  in  residential  water  (2,4  DNT  =  0.2  pg  LT1;  NG  =  1.5  pg  L1; 
NQ  =  1.6x103  pg  Lr1;  U.S.  EPA  2012)  and  2,4-DNT  is  classified  as  carcino¬ 
gen  (Human  Health  Screening  Values  2012)  making  their  fate  and 
transport  of  interest.  Limited  toxicological  data  for  NQ  suggests  it  is  far 
less  toxic  (Nipper  et  al.  2009). 

We  know  that  the  processes  controlling  the  dissolution  of  propellant  com¬ 
pounds  are  quite  different  from  those  controlling  explosive  dissolution 
(Lever  et  al.  2005;  Taylor  et  al.  2009a,  b).  Unlike  TNT  and  RDX,  nitrocel¬ 
lulose  (the  main  constituent  of  most  propellants)  does  not  dissolve  in  wa¬ 
ter.  While  solid  pieces  of  high  explosives  decrease  in  size  as  they  dissolve, 
many  nitrocellulose-based  propellants  do  not  change  shape  as  the  energet¬ 
ics  leach  from  their  NC  matrix.  Furthermore,  propellants  do  not  break 
apart  and,  because  propellant  grains  are  made  to  specification,  they  have  a 
narrower  size  distribution  than  explosives. 

How  much  2,4-DNT,  NG,  or  NQ  reaches  groundwater  depends  on  four 
main  factors.  First,  how  much  energetic  is  deposited  on  the  soil  surface  af¬ 
ter  firing  a  given  weapon  (deposition).  Second,  the  number  of  rounds  ex¬ 
pended  for  each  weapon  during  a  training  session  and  the  number  of 
training  sessions  per  year  (accumulation).  Third,  the  rate  at  which  2,4- 
DNT,  NG,  and  NQ  dissolve  from  the  residues  when  wetted  by  precipita- 
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tion.  Fourth,  the  fraction  of  the  dissolved  energetics  that  is  chemically  or 
biologically  degraded  in  the  soils  during  their  transit  to  groundwater. 

Previous  work  focused  on  quantifying  factors  l,  2,  and  4.  These  studies 
found  that  propellant  residues  are  scattered  onto  the  soil  surface  at  firing 
points  and  that  the  residues  consist  of  grains  or  fibers  of  nitrocellulose  that 
contain  energetic  compounds  (Walsh,  M.E.  et  al.  2007).  Jenkins  et  al. 
(2008)  quantified  the  concentrations  of  2,4-DNT,  NG,  and  NQ  in  soils  at 
firing  points  and  Walsh,  M.R.  et  al.  (2007)  published  robust  estimates  of 
the  amount  of  the  energetic  compounds  deposited  on  clean  snow  from 
single  firings  of  different  weapon  systems.  In  all  of  these  studies,  the  sam¬ 
ples  were  chemically  extracted  to  measure  the  concentration  of  energetics. 
Although  this  ensures  that  all  of  the  energetic  contained  within  the  nitro¬ 
cellulose  matrix  is  analyzed,  it  does  not  provide  information  on  how  quick¬ 
ly  or  how  much  of  the  energetic  is  dissolved  from  the  nitrocellulose  matrix 
by  rainfall. 

Once  in  solution,  2,4-DNT,  NG,  and  NQ  can  be  carried  by  infiltrating  pre¬ 
cipitation  to  the  groundwater.  During  transit  through  the  vadose  zone,  the 
energetics  can  interact  with  soil  in  a  variety  of  ways:  they  can  be  reversibly 
bound  to  the  soil,  irreversibly  bound,  or  transformed  altogether.  Kd,  the 
adsorption  or  partition  coefficient,  characterizes  the  affinity  of  the  energet¬ 
ics  to  sorb  onto  soil  organic  matter  or  mineral  surfaces.  When  the  relation¬ 
ship  between  the  concentration  on  the  soil  and  in  water  phases  is  linear 
with  respect  to  total  concentration,  the  partition  coefficient  can  be  de¬ 
scribed  as 

s  =  Kdc  (1) 

where 

s  =  soil  concentration  in  mg  kg_1 
Kd  =  adsorption  or  partition  coefficient  (L  kg  1  or  cm3  g  1) 
c  =  water  concentration  (mg  LT1) 

As  the  values  of  Kd  increase,  the  chemical  resides  mostly  on  the  solid  phase 
and  little  is  transported  downward  in  the  moving  pore  water. 

The  energetics  can  also  be  broken  down  chemically  or  biologically  while  in 
solution.  The  rate  at  which  the  energetic  is  mineralized  or  transformed  is 
related  to  its  concentration  in  solution  by  k,  the  reaction  rate  constant.  In 
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this  case  the  concentration  of  the  constituent  in  the  soil  water  and  the  rate 
of  loss  from  the  aqueous  phase  can  be  described  as 

dc/dt  =  -kc  (2) 


where 

c  =  concentration  of  the  specific  chemical  in  solution 

t  =  time 

k  =  reaction  rate  constant. 

Thus,  if  the  partition  coefficient  ( Kd )  and  the  rate  constant  k  are  known, 
the  fate  of  the  chemical  in  the  soil  system  can  be  modeled.  As  both  Kci  and 
k  vary  for  different  soils,  laboratory  batch  and  soil  column  experiments  are 
run  to  estimate  these  parameters. 

Different  types  of  laboratory  experiments  provide  estimates  for  these  pa¬ 
rameters.  To  estimate  Kd,  batch  experiments  shake  the  energetic  in  a  soil 
suspension  until  the  solution  concentration  remains  constant  (Pennington 
et  al.  1999,  2001;  Brannon  et  al.  2002;  Speitel  et  al.  2002;  Yamamoto  et  al. 
2004;  Mirecki  et  al.  2006)  (see  Tables  A2  and  A3).  Similar  experiments 
are  run  to  estimate  k,  but  the  concentrations  are  measured  as  a  function  of 
time  (Miyares  and  Jenkins  2000;  Jenkins  et  al.  2003;  Mulherin  et  al. 
2005).  In  soil  column  tests,  energetics  in  the  outflow  are  measured.  Some 
researchers  have  dripped  clean  water  on  columns  containing  contaminat¬ 
ed  soils  (Hewitt  and  Bigl  2005),  while  others  use  clean  soils  onto  which 
solutions  of  energetics  are  dripped  or  propellant  grains  or  residues  are 
leached  with  clean  water  (Dontsova  et  al.  2009).  The  Kd  values  measured 
for  NG  on  a  variety  of  soils  range  from  near  o  to  about  4  cm3  g“i.  These  low 
values  suggest  that  NG  should  not  be  retarded  moving  through  most  soils. 
However,  estimates  of  the  half-life  of  NG  in  various  soils  (t1/2  =  0.693  A) 
vary  from  much  less  than  an  hour  to  about  0.49  days  indicating  fast  deg¬ 
radation  (Jenkins  et  al.  2003).  Figure  1  summarizes  the  processes  we  think 
need  to  occur  if  propellant  energetics  are  to  reach  groundwater. 
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Figure  1.  Conceptual  model  showing  processes  needed  to  transport  energetic  compounds  to 
groundwater. 

The  NG,  2,4-DNT  and  NQ  are  always  found  in  NC,  indicating  they  are  not 
deposited  as  pure  product  and  must  dissolve,  diffuse,  or  desorb  from  their 
NC  matrices.  To  model  this  process,  we  need  information  on  how  these 
compounds  are  bound  in  and  released  from  the  nitrocellulose.  For  exam¬ 
ple,  does  the  NC  form  a  solid  through  which  the  energetics  have  to  diffuse 
or  do  the  NC  fibers  sorb  the  energetics  and  release  them  when  in  contact 
with  water?  The  latter  process  implies  two  or  more  distinct  phases.  We 
used  microscopic  techniques  to  help  distinguish  between  these  processes. 
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2  Materials  and  Methods 

We  used  field  collected  propellant  residues  and  unfired  propellant  grains 
for  laboratory  drip  and  batch  tests  (involving  no  soil)  and  for  soil  column 
studies.  Together,  these  data  allow  us  to  better  understand  the  dissolution 
rate  of  propellant-derived  2,4-DNT,  NG,  and  NQ  given  the  rainfall  rate 
and  provided  insight  into  how  well  vadose  transport  models  predict  disso¬ 
lution  rates  for  these  energetic  compounds. 

2.1  Propellant  samples 

Table  1  lists  all  the  propellants  tested  during  this  study.  A  list  of  nitrocellu¬ 
lose-based  propellants  used  by  the  U.S.  military  is  given  in  Table  Ai.  We 
examined  three  single-base  propellants,  eight  double-base  propellants  and 
three  triple-base  propellants.  Many  of  these,  Ml,  M2,  M9,  M10,  M45, 
WC844,  WC846,  WC860,  WPR289,  and  AT-4  propellants  were  collected 
from  live  fire  exercises  using  aluminum  trays  set  near  the  gun  muzzle  or  by 
melting  and  filtering  snow  that  contained  residues.  A  live-fire  test,  where 
triple-base  propellants  were  used,  produced  little  residue  in  the  particle 
traps  (Walsh,  M.R.  et  al.  2011). 

For  the  Ml  single-base  propellants,  we  collected  both  residues  and  unfired 
grains  from  single-  and  seven-perforation  Ml  propellants.  The  fired  single- 
perforation  propellant  deposited  sand-sized,  clear,  potassium  sulfate 
grains  that  contained  no  2,4-DNT  (Walsh,  M.E.  et  al.  2007)  and  conse¬ 
quently  the  residues  were  not  included  in  our  dissolution  tests.  The  seven- 
perforation  propellant,  on  the  other  hand,  left  slivers  of  NC  with  2,4-DNT 
concentrations  similar  to  those  of  the  unfired  propellant  (Walsh,  M.R.  et 
al.  2007).  We  have  tested  the  unfired  single  perforation  Mi  propellant  and 
both  the  unfired  seven-perforation  propellant  and  its  residues.  We  also 
have  unfired  and  fired  M10  propellant,  also  single-base,  but  as  it  only  con¬ 
tains  NC,  which  is  insoluble,  this  propellant  was  not  included  in  our  disso¬ 
lution  tests. 
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Table  1.  Propellants  tested  and  types  of  studies  performed  on  each  propellant.  All  samples  were 
photographed  optically. 


Propellant 

Weapon  System 

Energetic 

% 

Mass  (mg) 
per  grain 

Drip 

B 

s 

T 

M10-  grain 

60-,  81-mm  Mortar 

NC 

100 

X 

X 

Ml-1  hole 

155-mm  Howitzer 

NC  &  DNTs 

10 

10  ±  0.4 

U 

X 

X 

Ml- 7  hole 

105-mm  Howitzer 

NC  &  DNTs 

10 

-100 

u 

F 

U 

M9-1  hole 

81-mm  III.  Mortar 

NC  &  NG 

40 

1.23  ±  0.1 

u 

u 

X 

X 

AT-4-  sheet 

Rocket 

NC  &  NG 

36.5 

NA 

F 

M2-1  hole 

40-mm  grenade 

NC&  NG 

19.5 

0.32±0.06 

F 

F 

X 

X 

WPR289-grain 

9-mm  pistol 

NC&  NG 

12-18 

0.05±0.02 

u 

F 

U 

X 

X 

M45-1  hole 

120-mm  mortar 

NC&  NG 

10 

4.0±0.4 

u 

F 

u 

X 

WC844-grain 

M-16  Rifle 

NC&  NG 

9-11 

0.16±0.06 

u 

F 

u 

X 

WC846-grain 

7.62-mm  mach.  gun 

NC&  NG 

8-11 

0.18±0.06 

u 

F 

u 

X 

WC860-grain 

0.50-cal  mach.  gun 

NC&  NG 

8-11 

0.45±0.10 

u 

F 

u 

X 

X 

M31-1  hole 

Heavy  artillery 

NC&  NG  &  NQ 

20&55 

47  ±  1.2 

u 

u 

X 

X 

M31A1-7  hole 

Heavy  artillery 

NC&  NG  &  NQ 

20&55 

1405  ± 13 

u 

M31A2-7  hole 

Heavy  artillery 

NC&  NG  &  NQ 

20&55 

874  ±  7 

u 

u 

X 

U=Unfired,  F=  Fired,  B=  Batch,  S=SEM,  T=TEM 


For  the  double-based  propellants,  we  collected  unfired  grains  and  fired 
residues  for  most  of  the  formulations  (Table  l).  We  were  not  able  to  collect 
grains  of  unfired  M2  propellant  or  a  sheet  of  unfired  propellant  used  to 
fire  AT-4  rocket  as  the  unfired  propellant  could  not  be  safely  separated 
from  either  munition.  We  collected  unfired  M9  propellant  grains  (used  to 
fire  the  81-mm  illumination  rounds)  but  no  residues. 

For  the  triple-base  propellants  we  have  only  unfired  grains.  Although  we 
sampled  from  a  live  fire  test  that  used  M30  to  fire  a  self  propelled  Howitz¬ 
er,  the  test  yielded  no  residues  we  could  measure  (Walsh,  M.R.  et  al.  2011). 
We,  therefore,  obtained  three  triple-based  unfired  propellants  from  Rad¬ 
ford  ammunition  plant  and  tested  these  in  drip  and  batch  tests. 

For  all  of  the  field  tests,  we  collected  unfired  propellants  from  bullets  or 
bags,  whenever  possible,  and  included  them  in  our  dissolution  tests.  Alt¬ 
hough  unfired  propellants  are  not  intentionally  dumped  on  the  ground, 
excess,  unfired  propellant  bags  are  often  burned  in  the  field  a  process  that 
scatters  propellant  grains  (Walsh,  M.R.  et  al.  2009).  Furthermore,  com¬ 
parisons  of  the  dissolution  behavior  of  unfired  and  fired  grains  provide  in¬ 
sight  on  how  firing  changes  their  composition  and  dissolution  rate. 
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2.2  Microtome  sections  and  microscopy 

We  used  a  Leica  Ultracut  UCT  to  produce  thin  (2.5-15  pm)  slices  of  the 
propellants  for  optical  microscopy  and  TEM  analyses.  The  goal  was  to  map 
the  different  components  in  single-  (NC  with  or  without  2,4-DNT),  double- 
(NC  and  NG)  and  triple-  (NC,  NG,  and  NQ)  base  propellants. 

We  took  both  optical  and  scanning  electron  microscope  (SEM)  images  of 
the  propellants  and  their  residues.  The  color  of  the  propellants  can  be  seen 
in  the  optical  microscopy  images,  whereas  the  SEM  images  show  fine  de¬ 
tails  that  can  only  be  seen  at  high  magnification.  The  optical  images  were 
taken  at  CRREL  using  a  Leica  DMLM  microscope  equipped  with  a  digital 
camera.  We  used  an  FEI XL-30  field  SEM  at  Dartmouth  College  to  obtain 
high  magnification  images  of  the  propellant  surfaces.  The  SEM  had  both 
secondary  and  backscatter  electron  detectors  and  an  X-ray  microanalysis 
light  element  Si(Li)  detector  for  analyses  of  all  elements  heavier  than  car¬ 
bon. 

To  quantify  the  carbon,  nitrogen,  and  oxygen  and  map  their  distributions 
in  the  propellant  grains,  we  used  the  scanning  transmission  electron  mi¬ 
croscope  (SuperSTEM)  capable  of  ultra  high-resolution  imaging  (0.1  nm) 
located  at  Lawrence  Livermore  National  Laboratory.  The  instrument  im¬ 
ages  a  specimen  using  a  high  angle  annular  dark  field  (HAADF)  and  can 
analyze  sub-Angstrom  spots  with  an  electron  energy  loss  spectroscope 
(EELS).  The  specimens  are  generally  ultra-microtomed  into  thin  sections. 
The  SuperSTEM  is  an  80-300  keV  FEI  TITAN  equipped  with  a  mono¬ 
chromator,  dual  spherical  aberration  (Cs)  correctors,  a  EDAX  Genesis 
4000  Si(Li)  solid  state  energy  dispersive  x-ray  (EDX)  spectrometer,  and  a 
high-resolution  electron  energy-loss  spectrometer  (Tridiem  Gatan  Imag¬ 
ing  Filter).  TEM  brightfield/darkfield  images  and  STEM  high-angle  annu¬ 
lar  dark  field  (HAADF)  images  were  acquired  at  magnifications  of  10,000 
to  100,000 x.  Compositions  were  measured  using  EDX  and  electron  ener¬ 
gy-loss  spectroscopy  (EELS)  at  accelerating  voltages  of  80,  200,  and  300 
keV.  Energy-filtered  imaging  or  “EFTEM”  was  used  to  map  the  distribu¬ 
tion  of  C,  N,  and  O. 
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2.3  Dissolution  tests 

2.3.1  Drip  tests 

Before  beginning  the  drip  tests,  we  weighed  -50  mg  of  each  type  of  unfired 
propellant  and  residue  on  a  Mettler  A230  balance  (0.1  ±  0.1  mg)  and  ex¬ 
tracted  and  analyzed  these  to  determine  how  much  of  the  energetic  com¬ 
pound  was  present  as  a  function  of  their  mass.  Once  we  determined  the 
concentration  of  energetics  in  our  starting  materials,  we  placed  a  collec¬ 
tion  of  either  unfired  propellants  or  fired  residues  into  separate,  l-cm- 
diameter  Buchner  funnels  fitted  with  glass  frits  (Fig.  2).  A  syringe  pump 
was  used  to  drip  Milli-Q  water  (pH  =  6)  at  0.5  or  1.0  mL  hr  1  onto  the  pro¬ 
pellants.  The  water  flowed  through  the  frits  into  20-mL  scintillation  vials. 
We  replaced  the  vials  every  weekday,  measured  the  water  volume,  and  de¬ 
termined  the  concentration  of  the  propellant  constituents  using  high- 
performance  liquid  chromatography  (HPLC).  At  the  end  of  each  test,  we 
again  weighed  the  unfired  propellants  or  the  residues  and  extracted  them 
in  acetonitrile  to  determine  the  energetic  mass  remaining  in  the  NC  matrix 
and  estimate  a  mass  balance  for  the  test. 


Figure  2.  Laboratory  set-up  for  drip  tests  measuring  the  dissolution  of  energetics  from 
propellants. 


2.3.2  Batch  tests 

A  single-propellant  grain  was  massed  using  a  Mettler  Toledo  MX5  micro¬ 
balance  (1  ±  lpg)  and  imaged  using  an  optical  microscope  before  begin¬ 
ning  the  test.  Each  sample  was  then  placed  in  its  own  20-mL  glass  scintil- 


ERDC/CRREL  TR-12-9 


9 


lation  vial  and  deionized  water  (3  or  6  mL  depending  on  the  grain  size) 
was  added  by  volumetric  pipette.  The  vials  were  capped  and  placed  on  a 
Lab-Line  Instruments  Junior  Orbit  Shaker,  set  to  100  rpm. 

At  intervals  throughout  the  experiments,  the  water  in  the  vials  was  re¬ 
moved  and  replaced  with  new  deionized  water.  The  removed  water  (sam¬ 
ple)  was  then  massed  by  weighing  the  vial  it  was  decanted  into  before  and 
after  addition  of  the  water.  The  energetic  compounds  present  in  the  water 
were  detected  and  measured  using  HPLC.  The  tests  were  stopped  when 
the  concentration  of  the  energetic  in  the  water  samples  neared  the  detec¬ 
tion  limit.  The  propellant  samples  were  then  allowed  to  air  dry,  and  were 
re-imaged  using  the  same  optical  microscope. 

2.4  Propellant-soil  interactions 

2.4.1  Batch  soil  studies 

To  determine  how  2,4-DNT,  NG,  and  NQ  interact  with  soils,  we  measured 
the  changing  concentrations  of  aqueous  solutions  of  2,4-DNT,  NG,  and  NQ 
in  contact  with  clean  soils.  For  these  sorption  studies,  0.5  mL  of  a  10,000 
mg  L  1  stock  solutions  of  2,4-DNT,  NG,  and  NQ  were  added  to  250  mL  of  a 
0.05M  CaCl2  solution  to  create  a  2  mg  L_1  input  solutions.  Triplicate  tests 
were  conducted  using  4:1  solution/soil  mixture  of  the  three  propellant 
constituents  and  three  soil  types.  Each  propellant  constituent  was  studied 
separately.  The  three  soils  selected  for  the  study  were  Catlin  silt  loam 
(fine-silty,  mixed,  mesic,  superactive  Oxyaquic  Argiudolls),  Plymouth 
sandy  loam  (mesic,  coated  Typic  Quartzipsamments),  and  Sassafras  loam 
(fine-loamy,  siliceous,  mesic  Typic  Hapudults)  (Appendix  B).  Two  of  these 
soils,  Plymouth  sandy  loam  and  Sassafras  loam,  were  collected  on  military 
installations,  Massachusetts  Military  Reservation  and  Aberdeen  Proving 
Ground,  respectively.  We  also  used  these  soils  in  the  column  studies. 

Samples  were  placed  on  a  reciprocating  shaker  and  equilibrated  for  1,  2,  6, 
12,  24,  48,  72,  96,  and  120  h  and  then  centrifuged  for  30  min.  The  super¬ 
natant  solution  was  analyzed  for  target  compounds  by  HPLC  using  U.S. 
EPA  Method  8330b  for  2,4-DNT  and  NG,  and  the  Walsh,  M.E.  (1989) 
method  for  NQ.  Linear  distribution  coefficients  (IQ)  and  degrada¬ 
tion/transformation  constants  (k)  were  calculated  (eq  1  and  2). 
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2.4.2  Soil  column  tests 

To  determine  how  2,4-DNT,  NG,  and  NQ  dissolving  from  fired  propellants 
interact  with  soils,  we  conducted  saturated  flow  column  experiments  using 
propellant  residues  collected  at  the  same  time  as  those  used  in  drip  stud¬ 
ies.  Our  flux-controlled  flow-through  columns  consisted  of  Supelco  (Belle- 
fonte,  PA)  glass  tubes  (7-cm  high  with  1.18-cm  internal  diameter)  with 
PTFE  caps  (Fig.  3).  Silanized  glass  wool  was  used  on  top  and  bottom  to 
prevent  migration  of  the  particles.  A  Cole-Parmer  (Vernon  Hills,  IL)  Mas¬ 
ter  Flex  peristaltic  pump  was  attached  to  the  top  of  each  column  to  supply 
solution  at  a  given  flux.  Outflow  samples  were  collected  continuously  into 
4-mL  HPLC  vials  using  a  Teledyne  ISCO  (Lincoln,  NE)  Foxy  200  Fraction 
collector  with  a  200-vial  capacity.  This  experimental  setup  allowed  us  to 
run  up  to  10  columns  at  the  same  time. 


Figure  3.  Experimental  setup  used  for  the  soil  tests. 


Columns  were  packed  with  soil  and  saturated  from  the  bottom  up  with  a 
0.005-M  CaBr2  solution,  to  avoid  air  entrapment.  Pore  volume  (PV)  was 
determined  during  saturation  as  the  volume  of  solution  necessary  to  fill 
the  packed  column.  Propellants  were  then  placed  on  the  soil  surface,  the 
tops  of  the  columns  were  connected  to  the  pump  and  the  flow  was  started. 
Two  flow  rates  were  used:  0.01  mL  min  1  (0.55  cm  h_1  flux)  and  0.02  mL 
min"1  (1.1  cm  h_1  flux).  These  were  similar  to  fluxes  used  in  the  drip  stud- 
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ies.  We  tested:  Ml,  fired  (0.005  g  Per  column)  and  unfired  (one  particle  or 
0.115  ±  0.003  g  Per  column);  WC  860,  fired  (0.0166  g  per  column)  and 
unfired  (0.1  g  per  column);  and  M31  unfired  (two  particles  or  0.091  ± 
0.002  g  per  column).  After  about  35  to  60  PV,  the  propellants  were  re¬ 
moved  and  the  flow  continued  for  another  20-30  PV  to  evaluate  desorp¬ 
tion  of  the  energetics  from  the  soil.  We  used  a  conservative  tracer  (Br  )  an¬ 
alyzed  using  ion  chromatography  (Dionex  ICS  5000  with  electrical 
conductivity  detector)  to  distinguish  between  physical  and  chemical  non¬ 
equilibrium  processes  and  to  determine  dispersion  and  diffusion  in  the 
columns.  Outflow  from  the  columns  was  analyzed  using  HPLC  as  de¬ 
scribed  below. 

To  determine  how  much  2,4-DNT,  NG,  and  NQ  remained  in  the  soil  at  the 
end  of  the  tests,  we  extracted  the  soil  samples  for  18  h  in  acetonitrile  at  1:2 
soil/acetonitrile  ratio.  Extracts  were  diluted  with  water  to  achieve  1:3  ace¬ 
tonitrile  to  water  ratio,  filtered  through  0.45-pm  PTFE  filter,  and  analyzed 
using  HPLC,  as  described  below.  For  sets  IV  and  V,  we  performed  soil  ex¬ 
tractions  in  three  layers:  for  top,  middle,  and  bottom  of  the  soil  column. 

2.5  Analytical  methods 

We  quantified  2,4-DNT  and  NG  in  the  propellants  and  water  samples  us¬ 
ing  SW-846  Method  8330B  (EPA  2006)  and  NQ  using  a  separate  method 
(Walsh,  M.E.  1989).  For  method  8330B,  3  mL  of  the  water  sample  are 
added  to  the  1  mL  acetonitrile,  extracted,  and  filtered  through  a  0.45-iLim 
Millipore  cartridge.  High  Performance  Liquid  Chromatography  separates 
the  NG,  2,4-DNT,  and  their  co-contaminants  using  a  Water  NovaPak  C8 
column  eluted  at  1.4  mL  min"1  (28°C)  with  85:15  water:  isopropanol  mix 
and  detected  by  UV  at  254  or  210  nm.  For  NQ  the  water  sample  is  injected 
straight  onto  a  mixed-mode  PRi8/cation  column.  This  column  is  eluted 
with  1.5  mL  min"1  water  and  NQ  detected  by  UV  at  263  nm.  The  certified 
reporting  limit  for  water  was  5.0  pg  L  1  (Walsh,  M.E.  1989). 

Commercially  available  standards  (Restek),  specifically  developed  for 
Method  8330,  are  used  for  calibration.  We  prepared  1  and  10  ppm  8095A 
standards.  Sample  concentrations  of  around  1  ppm  are  expected,  but  if 
their  concentrations  are  outside  of  the  linear  range  for  the  calibration,  ~20 
ppm,  we  diluted  and  reanalyzed  the  samples.  The  l-ppm  standard  was 
used  to  recalibrate  the  instrument  after  every  tenth  samples.  Blanks  were 
run  before  each  standard  run  to  minimize  the  possibility  of  carryover, 
which  would  produce  a  poor  calibration.  The  10-ppm  standards  were  in- 
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terspersed  with  the  samples  as  unknowns,  and  a  blank  was  run  after  each 
to  minimize  carry  over.  The  concentrations  of  energetics  in  all  the  water 
samples  were  lower  than  solubilities  for  the  energetic  compounds. 
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3  Results  and  Discussion 

3.1  Propellant  characteristics 

3.1.1  Optical  images 

We  found  that  the  shape  of  the  original  propellant  grain  and  the  presence 
or  absence  of  holes  (perforations  made  to  increase  the  burn  rate)  dictated 
the  appearance  of  the  residue  (Fig.  4-8).  For  example,  propellants  without 
perforations  tended  to  leave  residues  that  were  smaller  versions  of  the 
original  propellant.  Propellant  grains  with  a  single  perforation  left  rings  or 
crescent  shaped  residues  (Fig.  6),  while  multiple  perforations  left  slivers. 


Figure  4.  Side  view  (a)  and  end  view  (b)  of  Ml  single-perforated  propellant.  Clear 
grains  deposited  from  firing  this  propellant  (c)  contain  potassium  and  sulfur,  but  no 
2,4-DNT. 
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a 


1  mm 


Figure  5.  Seven-holed  Ml  propellant  used  to  fire  the  105-mm  howitzer  rounds  (a) 
and  the  resultant  fiber  residues  (b).  Tetra  butyl  ammonium  hydroxide  placed  on  a 
fiber  (c)  shows  that  it  contains  2,4-DNT. 


ERDC/CRREL  TR-12-9 


15 


Figure  6.  Unfired  (a)  and  fired  (b)  M45  propellant  used  to  fire  a  120-mm  mortar  from 
a  Stryker  vehicle. 


Figure  7.  Unfired  (a)  and  fired  (b)  M9  propellant  used  to  fire  81-mm  illumination 
mortars. 


Figure  8.  Unfired  grains  (a)  and  fired  residues  (b)  used  to  fire  a  0.50-cal.  machine 
gun  (photographed  at  the  same  scale). 

For  propellants  composed  mainly  of  nitrocellulose  (>85%),  wetting  or 
submersion  did  not  change  the  shape  of  the  grain.  This  finding  is  reasona¬ 
ble  given  the  fact  that  nitrocellulose  is  hydrophobic  and  those  samples 
measured  have  extremely  small  pores  (10-30  nm  radii)  making  it  difficult 
for  water  to  penetrate  the  NC  (Ksiazczak  et  al.  2003).  For  the  M9  propel- 
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lant  that  contained  <60%  nitrocellulose,  however,  we  saw  a  collapse  of 
some  of  the  grains.  After  being  dripped  on  for  40  days,  the  M9  propellants 
developed  a  “skirt”  where  the  propellants  touched  the  glass  frit  (Fig.  9). 
This  “skirt”  was  all  that  remained  of  two  of  the  eight  M9  grains  tested.  We 
think  that,  as  the  NG  was  removed  from  the  propellant,  the  remaining  NC 
was  not  strong  enough  to  retain  its  shape  and  it  collapsed  into  a  disk  with 
a  diameter  slightly  larger  than  that  of  the  original  propellant  grain. 


Figure  9.  M9  unfired  propellants  used  in  drip  tests.  Note  “skirt”  that  formed  where 
the  propellant  contacted  the  glass  frit  (arrows). 


3.1.2  Thin  sections 

We  microtomed  thin  sections  for  optical  microscopy  and  TEM  analyses. 
We  found  that  single-  and  double-base  propellants  are  transparent  to 
translucent  and  show  signs  of  having  been  extruded  (Fig.  10). 


Figure  10.  Thin  section  of  M10  propellant  viewed  in  reflected  light 
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The  double-base  propellants  examined  (0.50-cal.,  9-mm,  M2,  and  M9) 
contained  micron-sized  inclusions  dispersed  throughout  (Fig.  11a  and  c). 
However,  the  outer  surfaces  of  these  propellants  are  devoid  of  inclusions  if 
the  propellant  has  been  exposed  to  water  during  drip  or  batch  tests  (Fig. 
lib  and  d).  We  think  these  inclusions  are  liquid  droplets  of  NG  contained 
within  the  nitrocellulose  matrix.  One  of  the  thin  sections  shows  what  ap¬ 
pears  to  be  a  liquid  phase  exiting  the  matrix  via  a  crack  (Fig.  lib). 


Figure  11.  Images  of  thin  sections  of  an  unfired  9-mm  propellant  placed  in  water. 
Optical  reflected  light  images  (a,  b),  optical  transmitted  light  image  (c)  and 
transmission  electron  image  (d).  The  dark  bars  in  b  and  d  are  the  TEM  grids  and  the 
white  areas  near  the  surface  of  the  propellant  in  d  are  rips  of  the  thin  section.  Note 
fluid  like  yellow  phase  moving  toward  propellant  surface  (arrow).  Note  also  that  areas 
close  to  the  marked  surfaces  of  the  propellants  do  not  have  droplet-like  features  like 
those  seen  in  the  interior  of  the  grain. 


The  triple-base  propellants  we  examined,  M31,  M31A1,  and  M31A2,  were 
cylindrically  shaped  grains  with  one  (M31)  or  seven  perforations  (M31A1 
and  M31A2)  (Fig.  12a).  Unlike  the  single-  and  double-base  propellants, 
which  contain  no  visible  crystals,  the  triple-base  propellants  have  crystals 
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of  NQ  (-5  |um  in  diameter)  running  parallel  to  the  long  axis  of  the  cylinder 
(Fig.  12b). 


Figure  12.  Triple-base  propellants  studied,  (a)  1.  M31,  2.  M31A1  and  3.  M31A2.  (b)  SEM  image  of 
M31A2  surface.  Dark  patches  are  the  graphite  coating,  long  crystals  are  NQ  and  the  light  grey  smooth 
looking  material  is  the  NC  and  NG  mix. 


Thin  sections  cut  perpendicular  to  the  long  axis  of  the  propellants  show 
that  triple-base  propellants  are  heterogeneous  in  appearance  and  brown  to 
orange  in  color  (Fig.  13a).  The  color  is  produced  by  2-nitrodiphenylamine, 
a  minor  ingredient  (1.5%  by  mass).  The  visual  heterogeneity  results  be¬ 
cause  the  2-nitrodiphenylamine  stains  NC  but  not  NQ.  The  NQ  crystals 
look  like  round,  grey  circles  at  high  magnification  (Fig.  13b). 


Figure  13.  Thin  section  of  M31  triple-base  propellant  viewed  in  reflected  light. 


We  sectioned  two  unfired  M31  propellants,  one  used  in  the  drip  tests  and 
the  second  stored  dry.  Images  of  the  sections  in  reflected,  darkfield,  and 
transmitted  light  (Fig.  14)  show  that  the  propellant  used  in  the  drip  test 
appear  to  have  lost  mass;  they  have  a  raggedy  looking  interior  surface.  We 
think  holes  around  the  exterior  surfaces  of  M31  grains  formed  when  NQ 
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crystals  dissolved  during  the  drip  tests  (Fig.  15a).  NQ  readily  dissolves  as 
evidenced  by  transparent  NQ  crystals  that  form  on  a  glass  slide  where  wa¬ 
ter,  that  had  wet  the  triple-base  propellant,  had  evaporated  (Fig.  15  b  and 
c). 


r 


500|im 


Figure  14.  Reflected,  dark  field  and  transmitted  light  of  (a)  thin  section  of  an  M31  not 
exposed  to  water  and  (b)  thin  section  of  an  M31  used  in  drip  test. 


Figure  15.  Appearance  of  M31  thin  section  exposed  to  water:  (left)  dark  field  view  of  yellow  NC+NG 
and  missing  NQ  crystals;  (center)  image  showing  wetting  front  and  NQ  crystals  that  formed  from  water; 
and  (right)  close  up  of  NQ  crystals. 


3.1.3  Scanning  electron  microscopy  images 

SEM  images  of  propellant  residues  show  that  the  residues  have  been  heat¬ 
ed  to  different  degrees  (Fig.  16).  These  images  show  two  interesting  things. 
First,  during  heating,  the  formation  and  release  of  gas  creates  pits  on  the 
residues  surface  (Fig.  16  and  17).  Second,  submicron  metal  grains,  present 
in  the  propellant,  become  concentrated  on  the  surface  of  the  residue  (Fig. 
16  and  17). 
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Figure  16.  Images  of  a  series  of  9-mm  propellant  grains  showing  how  the  grains  change  with 
increased  heating;  (a)  an  unfired  grain;  (b)  pits  visible  on  surface;  (c)  surface  has  begun  to 
bubble;  (d)  gas  bubbles  have  opened  larger  holes  on  surface;  (e)  metals  and  other  elements 
contained  within  the  propellant  are  concentrated  at  the  surface  as  the  propellant  burns;  (f) 
breakdown  of  nitrocellulose  structure. 


Figure  17.  High  magnification  images  of  surfaces  showing  (a)  the  initiation  of  pits  and  (b)  the 
bright  metal  particles  that  are  concentrated  on  the  residue  surface. 


SEM  images  of  small  arms’  unfired  propellants  immersed  in  water  for  5 
months  (batch  test  1)  showed  little  difference  from  those  stored  dry  (Fig. 
18)  despite  the  loss  of  between  7  and  21%  of  their  NG.  The  only  difference 
was  slightly  smoother  surfaces  attributable  to  dissolution  of  constituent 
salts  (e.g.,  sodium  sulfate,  calcium  carbonate,  potassium  nitrate)  present 
on  their  surfaces  (Fig.  18). 
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Figure  18.  SEM  image  of  WPR289  propellant  used  to  fire  9-mm  pistol:  (a)  surface  of  grain 
immersed  in  water  for  5  months,  (b)  surface  of  a  propellant  grain  not  placed  in  water. 


Figure  19.  Images  of  the  Ml  propellant:  (a)  and  (b)  were  not  wetted;  (c)  and  (d)  were 
immersed  in  water  for  500  days. 


However,  SEM  images  of  propellants  after  they  were  immersed  in  water 
for  500  days  (batch  test  2)  show  clear  changes  to  the  surfaces  of  the  grains 
For  example,  the  wet  Ml  single  perforation  grain  had  hollows  between  the 
ropey  strands  of  NC,  features  not  seen  on  dry  grains  (Fig.  19).  The  surface 
of  the  submerged  0.50-cal.  grain  had  rough  sections  that  look  sponge-like 
and  appear  to  extend  into  the  grain  (Fig.  20),  features  not  seen  on  dry 
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grains.  The  M9  grain  had  deep  holes  in  its  surface  where  material  was 
missing  (Fig.  21).  The  M31  surface  had  lost  much  of  its  carbon  coating  and 
the  outlines  of  dissolved  NQ  crystals  could  be  seen  over  the  entire  surface 
of  the  grain  (Fig.  22). 


Figure  20.  Images  of  unfired  0.50-cal.  propellants;  not  wetted  (a  and  b)  and  immersed  in 
water  for  500  days  (c  and  d). 
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Figure  21.  Images  of  unfired  M9  propellants;  not  wetted  (a  and  b)  and  immersed  in 
water  for  500  days  (c  and  d). 


Figure  22.  Images  of  unfired  M31  propellants;  not  wetted  (a  and  b)  and  immersed 
in  water  for  500  days  (c  and  d). 
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3.1.4  Scanning  transmission  electron  imaging  and  spectroscopy 

We  selected  WPR289  (9-mm  pistol)  and  M2  propellants  for  study  because 
of  their  different  NG  content:  WPR289  had  10%  NG  and  90%  NC  and  M2 
had  20%  NG  and  80%  NC. 

Energy-filtered  images  of  WPR289  show  a  distinct,  lighter  colored  zone 
near  the  graphite  coated  surface  in  an  unfired  grain  used  in  the  drip  tests 
(Fig.  23).  The  width  of  this  zone  is  ~io  pm.  The  propellant  rim  has  a  lower 
N/C  ratio  than  the  interior,  consistent  with  the  loss  of  nitroglycerin.  The 
sharp  boundary  is  also  consistent  with  the  behavior  of  NG  in  NC  and  in 
polymers,  generally  (Levy  1955).  A  study  of  diffusion  of  NG  into  NC  pro¬ 
pellants,  a  process  used  to  manufacture  extruded  impregnated  (El)  propel¬ 
lants,  showed  a  sharp  diffusion  front  (Levy  1955). 


Figure  23.  Energy  loss  image  of  sectioned  unfired  WPR289  propellant  immersed  in  water 
(left)  and  not  wetted  (right).  The  thin  graphite  layer  marks  the  surface  of  the  propellant.  Zones 
1  and  2  seen  in  the  left  image  are  interpreted  as  a  sharp  diffusion  boundary  between  NC 
lacking  NG  (1)  and  NC  containing  NG  (2).  We  think  the  stippled  pattern  seen  in  both  sections, 
but  prominent  in  the  right  image,  results  from  the  manufacture  of  grains  and  is  not  NG. 


Electron  energy-loss  spectroscopy  (EELS)  and  energy-filtered  images  of 
the  M2B  HEDP  propellant  (Fig.  24)  show  compositional  heterogeneities  in 
C,  N,  and  O  abundances  and  distributions.  The  surface  regions  appear 
more  homogenous  than  the  interior.  The  images  suggest  that  the  interior 
of  the  propellant  has  two  different  and  well-mixed  phases  with  essentially 
identical  compositions  but  slightly  different  densities.  We  also  noted  the 
presence  of  tiny  bright  spots  in  energy-filtered  image  that  exhibit  strong 
UV  absorption  at  207  nm  (Fig.  25).  Although  NG  absorbs  well  at  207  nm, 
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it  also  absorbs  at  shorter  wavelengths  and  we  suspect  these  spots  may  be 
constituent  salts  in  the  propellant. 


Region  1  - 
Interior 


Region  2  - 
interior 


Region  3  - 
Surface 


Figure  24.  Energy  loss  maps  made  of  the  surface  and  interior  regions  of  a 
sectioned  M2  propellant. 


Carbon 


Low-loss  EELS  spectrum  6eV(207  nm) 

ultraviolet  image 


Figure  25.  Figure  showing  an  energy  loss  map  for  carbon  (left),  an  EELS  spectrum  (center), 
and  a  map  of  the  sample  at  6eV  showing  very  bright  regions  (right). 
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3.2  Laboratory  drip  and  batch  tests 

The  laboratory  drip  tests  mimic  field  conditions  on  training  ranges,  where 
propellant  residues  are  scattered  on  the  soil  surface  and  whose  constitu¬ 
ents  are  dissolved  by  precipitation.  Here,  the  propellants  were  generally 
wet  but  could  experience  drying  during  the  weekends  when  the  syringe 
pumps  stopped.  The  tests  were  similar  to  those  successfully  used  to  meas¬ 
ure  dissolution  of  explosives  (Lever  et  al.  2005).  We  also  ran  batch  disso¬ 
lution  tests  on  individual,  unfired  propellant  grains  (Table  2).  These  tests 
allowed  us  to  investigate  energetic  loss  from  individual  propellants,  as 
compared  to  the  multiple  grains  needed  for  the  drip  tests.  The  fixed  water 
volume  and  weeklong  sampling  intervals  allowed  the  dissolved  energetic 
concentrations  to  build  up  sufficiently  to  measure.  The  tests  also  let  us 
study  many  different  types  of  propellants  simultaneously. 


Table  2.  Solubility  (mg  L_1  at  20  to  25°C)  and  molecular  formulae  for  energetic  compounds 
studied. 


Solubility 

Mol.  formula 

0/C 

N/C 

Cellulose 

NA 

C6H10O5 

0.83 

0 

Nitrocellulose 

NA 

C6N8H2O9 

1.5 

1.33 

2,4-DNT 

2701 

C7H6N2O4 

0.57 

0.29 

Nitroglycerine 

15002 

C3H5N3O9 

3 

1 

Nitroguanidine 

26003 

CH4N4O2 

2 

4 

iPhelan  and  Barnett  (2001);  2Yinon  (1999);  3Haag  et  al.  (1990). 


3.2.1  Single-base  propellants 

Single-base  propellants  containing  2,4-DNT  are  used  to  fire  artillery, 
mainly  105-  and  155-mm  rounds  from  howitzers.  Figure  26  shows  the  per¬ 
centage  of  DNTs  dissolved  during  drip  and  batch  tests  on  a  variety  of  Ml 
propellants.  Note  that  for  the  unfired  grains,  the  dissolution  rate  was  quite 
linear  and  that  even  after  500  days  the  maximum  DNT  loss  is  only  10%. 
The  larger  propellant  grains  lost  more  DNT  than  the  smaller  grains  but  a 
smaller  percentage  of  what  they  contained  (Table  3).  The  residues  from  Ml 
seven-perforation  propellant  (12  fibers)  lost  the  highest  percentage  of  their 
2,4-DNT  owing  the  their  large  surface  to  volume  ratios,  although  the  total 
mass  was  small.  The  mass  loss  curve  for  the  residues  was  not  linear— it 
rose  rapidly  initially,  but  after  day  20  it  becomes  more  linear,  still  having  a 
positive  slope  (Fig.  26). 
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Table  3.  Mass  balance  calculations  for  single-based  propellants  used  in  the  drip  and  batch  experiments. 


Propellant 

single-base-2, 4-DNT 

10±2 %  DNT 

Mass 
initial  (/) 
(mg) 

Mass 
final  (F) 
(mg) 

Initial 
-  final 

(l-F) 

(mg) 

Mass 

diss. 

(D) 

(mg) 

Mass 

ext. 

(E) 

(mg) 

%  Mass 
balance 

([F+D]/I)x 

100 

%  2,4-DNT 

Recovered 

«D+E]/Exp) 

xlOO 

%2, 4-DNT 
Measured* 

2,4-DNT 
/Prop  xlOO 

M  1-seven  perforations 

Ml-7  unfired  extracted 

115.3 

11.21 

9.72 

Ml-7  fired  extracted 

3.74 

0.30 

8.01 

Ml-7  (drip) 

113.3 

112.9 

0.4 

0.20 

100 

M-7  fired  (drip) 

3.9 

0.07 

0.17 

75.38 

Weathered-A  (batch) 

245.6 

243.5 

2.1 

0.43 

99 

Weathered-B  (batch) 

238.6 

237.3 

1.3 

0.34 

100 

Weathered-C  (batch) 

238.4 

237 

1.4 

0.12 

99 

Ml-one  perforation 

Ml  unfired  extracted 

30.7 

2.78 

9.06 

M1-1A  (batch-unfired) 

10.7 

10.5 

0.2 

0.10 

99 

M1-1B  (batch-unfired) 

10 

9.8 

0.2 

0.11 

99 

Ml- 1C  (batch-unfired) 

10 

9.8 

0.2 

0.10 

99 

From  Technical  manuals,  *  analyzed  at  CRREL.  F  =  final  mass,  D  =  dissolved  mass,  /  =  initial  mass,  Exp  =  expected,  Prop  = 
mass  of  propellant. 


Time  (days) 

Figure  26.  Percentage  of  2,4-DNT  and  2,6-DNT  dissolved  into  water  from 
unfired  seven  perforation  Ml  propellant  grains  (black)  and  from  fired  residues 
(red). 
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3.2.2  Double-base  propellants 

Double-base  propellants  are  the  most  common  and  are  used  to  fire  small 
arms,  mortars,  and  rockets.  In  both  drip  and  batch  tests,  unfired  double¬ 
base  propellants  show  initial  rapid  dissolution  of  NG  followed  by  much 
slower  dissolution.  Most  double-base  propellants  lost  NG  in  proportion  to 
how  much  NG  they  contained.  For  example,  the  M9  propellant  with  40% 
NG  lost  a  greater  percentage  of  its  NG  than  a  propellant  that  contained 
less  (Fig.  27).  We  found  that  the  amount  of  NG  dissolved  was  a  function  of 
the  NG/NC  ratio  in  the  propellant  (Fig.  28).  Clustering  of  the  data  using 
this  normalization  suggests  that  NC  binds  10  to  20%  of  the  NG  and  that 
any  extra  NG  is  not  well  retained. 


0  50  100  150  200  250 


Time  (Days) 

Figure  27.  Percentage  of  NG  dissolved  (mg  dissolved/mg  present  in  propellant) 
for  a  variety  of  double-based  propellants.  We  show  the  single-base  Ml  for 
comparison. 

Exceptions  from  this  trend  are  the  ball  propellants,  used  to  fire  small 
arms,  and  the  M45  propellant,  used  to  fire  mortars  from  the  Stryker.  The 
ball  propellants  all  contain  ~io%  NG,  yet  variable  amounts  of  NG  are  dis¬ 
solved  independent  of  their  NC  content.  The  M45  propellant  also  contains 
10%  NG,  yet  it  loses  less  than  1%  during  both  the  drip  and  batch  tests  (Ta¬ 
ble  4,  Fig.  28).  The  M45  is  a  squat  grain  with  a  central  perforation.  We  ex¬ 
pected  the  central  perforation  to  increase  its  surface  area  and,  therefore, 
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increase  its  NG  loss.  In  the  case  of  M45,  we  speculate  that  the  nitrocellu¬ 
lose  was  not  fully  nitrated  when  it  was  manufactured  so  that  the  NG  effec¬ 
tively  binds  to  the  NC. 
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Figure  28.  Plot  showing  the  %  NG  dissolved  normalized  by  the 
NG/NC  ratio  versus  time.  All  of  the  data  from  NG  containing 
propellants,  except  for  the  small  arms,  are  plotted  on  this  figure. 


Table  4.  Mass  balance  calculations  for  double-based  propellants  used  in  the  drip  and  batch  experiments. 


Propellant-NG 

Double-base  NG 

Mass  initial 

(/) 

(mg) 

Mass  final 

(F) 

(mg) 

Initial- 

final 

(mg) 

Mass  diss. 

(D) 

(mg) 

Mass  ext. 

(E) 

(mg) 

%  Mass 
balance 

([F+DV  l)x 
100 

%  NG 
recovered 

([D+EV  Exp)x 
100 

%  NG 

measuredt 

(ng//; 

xlOO 

%  NG  ex¬ 
pected** 

M2-one  perforation 

19.5±1 

M2  extracted  (R 

3.33 

0.64 

19.25 

M2C  (drip-F) 

29.49 

24.55 

4.94 

1.13 

3.58 

87 

83 

M2D  (drip-F) 

4.56 

3.97 

0.59 

0.19 

0.61 

91 

91 

M2A  (batch-F) 

0.3 

0.2 

0.1 

0.01 

0.03 

71 

76 

M2B  (batch-F) 

0.4 

0.3 

0.1 

0.02 

0.05 

80 

86 

M2C  (batch-F) 

0.3 

0.3 

0 

0.02 

0.04 

106 

106 

M9  solid  grain 

40  ±1.5 

M9  extracted  (U) 

117.3 

46.00 

39.22 

M9  extracted  (F) 

7.3 

1.47 

20.14 

M9  (drip-6-U)* 

7.2 

3.53 

3.67 

1.39 

0.71 

68 

73 

M9  (drip-2-U) 

2.6 

2.08 

0.52 

0.49 

0.46 

99 

91 

M9  (batch-U)  A 

1.3 

1.07 

0.23 

0.32 

107 

sent  to  LLNL 

M9  (batch-U)  B 

1.3 

1.13 

0.17 

0.33 

0.20 

112 

101 

M9  (batch-U)  C 

1.1 

1.05 

0.05 

0.28 

0.18 

121 

105 
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Propellant-NG 

Double-base  NG 

Mass  initial 

(/) 

(mg) 

Mass  final 

(F) 

(mg) 

Initial- 

final 

(mg) 

Mass  diss. 

(D) 

(mg) 

Mass  ext. 

(E) 

(mg) 

%  Mass 

balance 

«F+D]/I)x 

100 

%  NG 
recovered 

([D+E]/Exp)x 

100 

%  NG 

measuredt 

0 NG/I ) 
xlOO 

%  NG  ex¬ 
pected** 

M45-one  perforation 

10  ±2 

M45  extracted  (U) 

4.47 

9.95 

M45  extracted  (U) 

4.79 

10.25 

M45  extracted  (U) 

3.85 

9.57 

M45  extracted  (F) 

4.1 

0.36 

8.89 

M45  extracted  (F) 

13.71 

8.9 

M45  (drip-F) 

14 

13.71 

0.29 

0.01 

0.85 

98 

61 

M45  (batch-5  U) 

19.3 

19.3 

0 

0.02 

1.88 

100 

98 

WPR289  (9-mm) 

12-18 

WPR289  (ext.  U) 

12.20 

WPR289  (ext.  F) 

48.75 

6.24 

12.80 

WPR289  (drip-F) 

38.6 

30.6 

8 

0.22 

5.94 

80 

125 

WPR289  (drip-U) 

14.5 

11.1 

3.4 

0.29 

1.07 

79 

77 

WC844  (5.56mm) 

8-11 

WC844  (ext.  U) 

9.90 

WC844  (ext.  F) 

54 

4.4 

8.10 

WC844  (drip-F) 

68.2 

62 

6.2 

0.33 

5.62 

91 

97 

WC844  (drip-U) 

32.8 

31.2 

1.6 

0.65 

2.09 

97 

103 

WC846  (7.62mm) 

8-11 

WC846  (ext.  U) 

10.20 

WC846  (ext.  F) 

47 

3.70 

7.90 

WC846  (drip-F) 

17.9 

0.10 

lost 

WC846  (drip-U) 

17.7 

0.34 

1.27 

89 

WC860  (0.5  cal) 

8-11 

WC860  (ext.  U) 

9.70 

WC860  (ext.  F) 

61.5 

4.80 

7.80 

WC860  (drip-F) 

57 

49.8 

7.2 

0.27 

4.72 

88 

112 

WC860  (drip-U) 

20.2 

19.1 

1.1 

0.47 

1.21 

97 

86 

WC860  (batch-U) 

2.1 

2.1 

0 

0.08 

0.12 

104 

97 

AT-4-sheet 

36.5-39.5 

AT-4  ext.  F 

0.45 

0.16 

35.56 

AT-4  ext.  F 

71 

26.00 

36.62 

AT-4  drip  F&W 

3.84 

3.52 

0.32 

0.40 

0.63 

102 

75 

AT-4  drip  F 

7.97 

5.31 

2.66 

1.78 

0.86 

89 

92 

AT-4  drip  F 

6.4 

4.21 

2.19 

1.09 

0.80 

83 

82 

*Four  of  6  initial  grains  recovered.  E  =  extracted,  U  =  unfired,  F  =  fired,  W  =  weathered,  F  =  final  mass,  D  =  dissolved  mass, 
/  =  initial  mass,  Exp  =  expected  mass,  tmeasured  at  CRREL,  **  from  propellant  manual. 
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We  collected  both  fired  residue  and  unfired  grains  from  all  the  small  arms 
propellants.  Concentrations  of  NG  in  the  unfired  grains  were  within  the 
variability  given  in  the  technical  manual,  whereas  the  fired  residues  con¬ 
tained  about  80%  of  the  NG  on  a  mass  basis  (Table  3).  Figure  29  shows 
that  more  NG  was  dissolved  from  the  unfired  propellants  (15-20%)  than 
from  their  residues  (3-7%).  For  the  unfired  propellants,  the  shape  of  the 
cumulative  mass  loss  versus  time  curves  were  consistent  with  rapid  loss  of 
the  NG  from  the  surface  of  the  grain,  followed  by  slower  diffusion  of  the 
NG  from  the  interior  of  the  grain.  The  high  aqueous  solubility  of  NG  sug¬ 
gests  that  NG  could  be  rapid  dissolved  by  contact  with  water.  If  NG  existed 
as  fine  liquid  droplets  within  an  NC  matrix  rather  than  as  dispersed  mole¬ 
cules,  droplets  at  the  grain  surface  would  be  in  direct  contact  with  water 
whenever  the  grain  was  wet.  Once  this  outer  layer  of  liquid  NG  was  deplet¬ 
ed,  NG  would  need  to  diffuse  through  the  NC  matrix  to  reach  the  water, 
with  considerably  lower  diffusivity  (-10-14  cm2s_1).  Late  time  dissolution 
would  thus  be  limited  by  molecular  diffusion. 


0  20  40  60  80  100 


Time  (days) 

Figure  29.  Percent  NG  dissolved  versus  time  for  unfired  and 
fired  (F)  small  arms  propellants. 

This  concept  also  qualitatively  accounts  for  the  much  lower  dissolution 
rates  observed  for  fired  grains.  Firing  likely  burns  or  volatilizes  surface  NG 
droplets,  and  we  estimate  that  volatilizing  NG  from  a  surface  layer  about 
5%  of  the  thickness  of  the  fired  particle  would  lower  the  NG  concentration 
by  20%,  the  difference  we  measured  between  fired  and  unfired  propel¬ 
lants.  SEM  images  of  the  small  arm  propellants  show  pits  that  are  about  5 
pm  in  diameter  and  similarly  deep  (Fig.  16),  consistent  with  the  depth 
needed  to  produce  the  drop  in  concentration.  All  dissolution  of  NG  from 
fired  grains  would  thus  be  limited  by  molecular  diffusion.  Indeed  the  data 


ERDC/CRREL  TR-12-9 


32 


for  fired  grains  did  not  show  the  abrupt  roll-off  characteristic  of  the  un¬ 
fired  grains,  but  were  similar  to  the  latter’s  late-time  dissolution  rates.  The 
linear  shape  of  the  cumulative  mass  loss  curves,  the  slower  dissolution  rate 
of  NG  from  the  fired  residues,  and  their  20%  lower  NG  concentration 
compared  to  unfired  grains  are  explained  if  NG  near  the  surface  is  burned 
during  firing. 

The  shapes  of  the  mass  loss  curves  for  most  unfired  double-base  propel¬ 
lants  were  different  from  those  seen  for  single-base,  DNT-containing  pro¬ 
pellants:  they  were  much  steeper  initially  and  then  flattened  out.  We  think 
that  this  results  because  NG  is  added  as  a  liquid  to  the  NC  during  manu¬ 
facture.  Depending  on  the  nitration  level  of  the  NC,  variable  amount  of  NG 
are  attached  to  the  NC.  Any  excess  NG  migrates  toward  the  surface  of  the 
grain  where  water  contacting  the  liquid  NG  removes  it  from  the  propel¬ 
lant. 

Dissolution  from  fired  versus  unfired  are  also  very  different  for  single  and 
double-base  propellants.  For  the  single-base  propellants,  the  fired  residues 
lose  a  larger  percentage  of  2,4-DNT  than  do  the  unfired  grains.  Further¬ 
more,  the  %mass  loss  for  the  residues  was  much  higher  at  the  beginning  of 
the  tests  than  at  the  end.  We  think  this  occurred  because  the  residues  were 
small  and  had  much  larger  surface  to  volume  ratios,  allowing  the  2,4-DNT 
to  diffuse  to  the  surface.  For  double-base  propellants,  except  for  the  M45 
whose  mass  loss  was  similar  to  single-base  propellants,  the  unfired  dou¬ 
ble-base  propellants  lost  a  higher  percentage  of  their  energetics  (15-20%) 
than  their  fired  residues  (3-7%).  This  suggests  that  the  NG  from  the  outer 
surface  of  the  double-base  propellants  is  easily  lost  relative  to  2,4-DNT, 
and  that  once  lost,  the  rate  at  which  NG  diffuses  out  to  the  surface  is  slow¬ 
er  than  the  2,4-DNT  diffusion  rate. 

3.2.3  Triple-base  propellants 

Like  the  single-base  propellants,  triple-base  propellants  are  also  used  to 
fire  artillery.  We  ran  drip  tests  for  140  days  on  eight  individual  M31  grains, 
four  at  0.5  mL  hr  1  and  four  at  1.0  mL  hr  1  flow  rates.  During  this  time  the 
M31  propellants  lost  an  average  of  6.32±o.22  mg  of  NG  (or  -73%)  and 
3.23±o.55  mg  of  NQ  (or  ~13%).  The  mass  loss  scaled  linearly  with  the  drip 
rate,  so  the  data  from  the  0.5  and  1.0  mL  hr  1  flow  rates  were  plotted  to¬ 
gether  (Fig.  30).  Batch  tests  run  on  three  M31  grains  showed  similar  mass 
losses,  6.87  ±  0.23  mg  NG  and  3.03  ±  0.38  mg  NQ  (Fig.  31— black  curves), 
although  the  batch  tests  ran  for  almost  four  times  as  long.  Note  that  most 
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of  the  NG  was  dissolved  rapidly,  as  was  seen  in  the  double-base  propel¬ 
lants.  The  NQ  also  dissolves  rapidly  initially  but  the  late  time  dissolution 
for  the  NQ  was  smaller  than  that  for  NG. 


Time  (days) 


Figure  30.  Drip  test  results  showing  the  percent  NG  and  NQ 
dissolved  (Ave.  ±  1  sigma,  n  =  8)  versus  time  for  unfired  M31 
single-perforation  propellants. 


Time  (Days) 


Figure  31.  Batch  tests  results  showing  the  percent  NG  and  NQ 
dissolved  (Ave.  ±  1  sigma,  n  =  3)  versus  time  for  unfired  M31 
(black),  M31A1  (red)  and  M31A2  (green)  triple-base  propellants. 

The  lower  of  the  pair  of  curves  is  always  NQ. 

Although  there  is  more  NQ  than  NG  in  the  M31  propellant  (55  vs.  20%), 
and  NQ  is  more  soluble  than  NG  (2600  vs.  1500  mg  L  1  at  20°  C,  Haag  et 
al.  1990;  Yinon  1999),  both  by  mass  and  by  percentage,  more  NG  was  dis¬ 
solved  than  NQ.  We  think  this  is  because,  during  manufacture  of  the  pro¬ 
pellant,  NG  is  added  as  a  liquid,  whereas  NQ  is  mixed  in  as  a  solid.  Tests  to 
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measure  how  well  the  components  in  triple-based  propellants  were  mixed 
(Yazici  and  Kalyon  1998)  showed  that  NG  is  not  as  well  mixed  as  NC  and 
NQ  and  that  there  is  more  NG  near  the  surface  and  less  in  the  interior  of 
the  grain.  These  authors  suggest  that,  after  a  certain  threshold  (27%NG  for 
a  12.2%  nitrated  NC),  the  NG  does  not  effectively  bond  to  the  NC  and  it 
migrates  to  the  propellant  surface  as  a  low  viscosity  fluid.  This  migration 
would  make  liquid  NG  available  near  the  surface  of  the  propellant  where  it 
would  be  removed  when  in  contact  with  water.  The  NQ  on  the  other  hand 
would  have  to  dissolve  to  leave  the  propellant. 

The  M31A1  and  M31A2  propellants  are  similar  in  composition  to  the  M31 
(54%  NQ  and  20%  NG)  but  are  larger,  cylindrically  shaped  grains  with 
seven  perforations  (Fig.  12a).  Drip  and  batch  tests  run  on  M31A1  again 
showed  higher  losses  of  NG  (80%)  than  NQ  (20%)  (Fig.  31).  The  multi- 
perforated  grains  showed  a  different  dissolution  pattern  in  that  the  origi¬ 
nal  loss  of  NG  and  NQ  was  less  abrupt  and  the  late  time  dissolution  rates 
for  both  NG  and  NQ  were  higher  than  for  the  single-perforation  grains. 
Both  these  multi-perforated  grains  were  large  (1.4-  and  0.87-g  masses; 
Table  4)  and  had  lower  surface  to  volume  ratios. 

3.2.4  Mass  balance 

Tables  3-5  list  the  initial  and  final  masses  of  the  propellants  tested  in  the 
drip  and  batch  tests  along  with  the  masses  of  energetic  compounds  dis¬ 
solved  and  extracted  at  the  ends  of  the  tests.  The  masses  of  propellants  or 
residues  extracted  to  determine  the  percentage  of  energetics  initially  pre¬ 
sent  in  the  propellant  are  underlined.  We  calculated  a  percent  mass  bal¬ 
ance  for  all  the  propellants  by  summing  the  final  mass  and  the  dissolved 
mass  and  dividing  by  the  initial  mass.  For  the  double-base  and  triple-base 
propellants,  we  also  calculated  an  energetics  percent  mass  recovery  by 
adding  the  dissolved  and  the  extracted  energetic  compounds  and  dividing 
by  the  initial  mass  of  energetics.  The  latter  value  was  determined  by  mul¬ 
tiplying  the  initial  mass  of  the  propellant  or  residue  by  the  fraction  of  en¬ 
ergetic  extracted  from  a  similar  sample. 
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Table  5.  Mass  balance  calculations  for  triple-based  propellants  used  in  the  drip  and  batch  experiments. 


Propellant— Triple 
Base  NQ+NG 

Mass  initial 

(0 

(mg) 

Mass 

final 

(1 F) 

(mg) 

Initial  - 

final 

(mg) 

Mass  diss. 

(D) 

(mg) 

Mass  ext. 

(E) 

(mg) 

%  Mass 

balance 

([F+D]//)x 

100 

%  NG  recovered 

([D+E  ]/Exp)x 
100 

M31-one  hole 

NG+NQ 

NG+NQ 

NG 

M31-1  (drip-U) 

45.61 

35.2 

10.41 

9.90 

9.40 

99 

90 

M31-2  sent  to  LLNL 

M31-3  (drip-U) 

45.9 

35 

10.9 

9.54 

9.34 

97 

85 

M31-4  (drip-U) 

45.56 

34.6 

10.96 

9.52 

9.35 

97 

80 

M31-5  (drip-U) 

42.7 

30.9 

11.8 

10.68 

8.98 

97 

84 

M31-6  (drip-U) 

44.24 

34.4 

9.84 

8.74 

9.42 

98 

87 

M31-7  (drip-U) 

45.21 

34.8 

10.41 

9.18 

9.68 

97 

86 

M31-8  (drip-U) 

43.28 

32.3 

10.98 

9.84 

9.53 

97 

88 

M31-A  (batch-U) 

46 

35.1 

10.9 

10.04 

98 

76 

M31-B  (batch-U) 

46.1 

35.2 

10.9 

10.39 

99 

80 

M31-C  (batch-U) 

48.2 

36.3 

11.9 

11.19 

99 

76 

M31Al-seven  holes 

NG 

M31A1-A  (drip-U) 

1450 

1351.4 

98.6 

198.94 

107 

65 

M31A1-B  (drip-U) 

1400 

1261.3 

138.7 

242.64 

107 

80 

M31A1-C  (drip-U) 

1430 

1308.6 

121.4 

242.17 

108 

78 

M31A1-D  (drip-U) 

1320 

1203.9 

116.1 

215.22 

108 

74 

M31A1-A  (batch-U) 

1392.2 

1119.1 

273.1 

424.16 

111 

85 

M31A1-B  (batch-U) 

1417.8 

1164 

253.8 

394.17 

110 

84 

M31A1-C  (batch-U) 

1404.5 

1104.3 

300.2 

449.07 

111 

86 

M31A2-seven  holes 

M31A2-A  (batch-U) 

875.3 

510.1 

365.2 

413.25 

105 

87 

M31A2-B  (batch-U) 

866.2 

526 

340.2 

380.40 

105 

89 

M31A2-C  (batch-U) 

879.3 

530.1 

349.2 

365.39 

102 

78 

The  percent  mass  balances  for  the  single-base  propellants  average  99±i% 
(Table  3),  for  the  double-base  propellants  93±i3%  (Table  4),  and  for  the 
triple-base  propellants  103±5%  (Table  5).  The  percent  mass  balances  are 
generally  good.  Low  mass  balances  resulted  if  either  water  or  propellant 
samples  were  lost.  For  example,  the  glass  frit  in  our  drip  tests  clogged  from 
time  to  time,  causing  the  funnels  to  overflow,  thereby  losing  both  water 
and  propellant  grains.  Small  propellants  were  subject  to  loss  in  this  man¬ 
ner,  which  resulted  in  low  percent  mass  balances  (e.g..  the  M9  sample 
where  six  grains  were  tested  but  only  four  recovered).  High  mass  balances 
also  were  measured  and  these  occurred  mainly  for  large  propellants.  For 
example,  the  M31A1  propellants  are  ~io%  high.  Although  we  waited  to 
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mass  the  grains  for  about  a  week  after  the  end  of  the  tests,  we  think  the 
perforations  in  these  grains  still  contained  some  residual  water.  Owing  to 
the  residual  water,  their  final  masses  were  high  and  their  initial  minus  fi¬ 
nal  mass  was  smaller  than  it  should  have  been. 

For  the  double-base  propellants,  the  %NG  balance  average  is  gi±i5%,  a 
value  that  seems  reasonable  given  that  we  do  not  know  the  exact  concen¬ 
tration  of  NG  in  our  starting  materials.  For  example,  for  the  sample  that 
had  the  highest  value,  125%  (WPR289,  fired-9mm  residue),  the  %NG 
would  vary  between  133  and  89%  depending  on  if  12  or  18%  where  used  as 
the  starting  composition— the  range  given  for  this  propellant. 

We  calculated  the  %NG  and  NQ  balance  for  the  single  perforation  M31 
grains  but  had  difficulty  extracting  these  samples.  The  energetic  balance 
for  the  seven  M31  grains  averaged  78±4%,  the  NG  recovered  was  86±3%, 
and  the  %NQ  47±3%.  As  the  NQ  is  not  as  soluble  as  NG  in  ACN,  we  think 
that  the  transparent  NQ  crystals  remained  unnoticed  in  the  bottom  of  the 
vial.  A  better  extraction  method  is  needed  for  triple-base  propellants. 

3.3  Diffusion  modeling  of  NG  from  small  arm  propellants 

Lever  et  al.  (2005)  and  Taylor  et  al.  (2009a,  b)  showed  some  success  in 
modeling  dissolution  of  particles  of  high  explosives  exposed  to  rainfall.  In 
those  cases,  explosive  at  a  particle’s  surface  are  in  direct  contact  with  liq¬ 
uid  water,  and  as  they  dissolve  the  particle  becomes  smaller.  However,  dis¬ 
solution  of  2,4-DNT,  NG,  and  NQ  from  propellant  grains  is  quite  different. 
The  NC  is  not  soluble  and  we  needed  to  learn  if  the  energetic  components 
were  transported  via  molecular  diffusion  through  an  effectively  solid  NC 
matrix  or  whether  water  can  penetrate  the  matrix  allowing  the  NG  to  dif¬ 
fuse  through  tortuous  water  channels  to  the  grain’s  surface.  Our  laboratory 
data  on  cumulative  dissolved  mass  loss  versus  time  cannot  differentiate 
between  these  two  mechanisms  and  we  found  no  data  on  these  processes 
in  the  literature.  We  focused  our  modeling  effort  on  understanding  NG 
loss,  as  NG  is  present  in  many  propellants,  including  triple-base  propel¬ 
lants.  For  simplicity  we  chose  the  small  arms  propellants  and  we  modeled 
diffusion  through  a  solid  matrix.  We  modeled  both  constant  diffusion  and 
variable  diffusion. 

First,  we  modeled  propellant  dissolution  as  diffusion  through  a  solid  ma¬ 
trix,  using  simple  grain  geometries  and  a  constant  diffusion  coefficient. 
Figure  32  shows  model  results  for  the  solid  diffusion  model  applied  to 
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small  arms  propellants,  both  from  batch  tests  on  single  grains  and  drip 
test  on  many  grains. 


0  50  100  150  200 

Time  (d) 

b.  Drip  tests  (many  grains). 


Figure  32.  Laboratory  data  from  batch  tests  (single  grain)  and 
drip  tests  (many  grains)  compared  with  model  predictions 
(solid  lines)  of  solid-phase  diffusion  of  NG  from  spherical  NC 
grains.  We  fit  the  model  results  to  the  long-term  trends  in  the 
data  to  obtain  diffusion  coefficients.  Initial  dissolution  rates 
are  higher  than  predicted  for  spherical  particles,  and  we  are 
investigating  effects  of  grain  shape  and  surface  roughness. 


The  model  assumes  spherical  particles  of  equal  mass  to  that  of  the  propel¬ 
lant  grains.  We  fit  the  model  results  to  the  long-term  trends  (tail  slopes)  in 
the  data  to  obtain  diffusion  coefficients.  The  diffusion  coefficients  are  dif- 
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ferent  for  each  propellant  type  (9  mm  =  3.20xio_14  cm2  s_1;  5.56  m  = 
5.90xio_14  cm2  s'1;  7.62  =  i.ooxio_13  cm2  s”1;  and  0.50  cal  =  i.20xio_13 
cm2  s_1)  but  the  same  for  the  two  tests  on  the  same  propellants.  The  results 
suggest  that  dissolution  is  similar  to  diffusion  from  a  sphere  as  time  pro¬ 
gresses.  However,  dissolution  rate  is  higher  at  early  times  than  diffusion 
from  a  sphere. 

The  results  found  by  Levy  (1955)  suggested  that  concentration-dependent 
diffusion  of  NG  from  the  NC  matrix  could  control  the  NG  dissolution  rate. 
To  test  this,  we  examined  the  effect  of  concentration-dependent  diffusivity 
compared  with  constant  diffusivity,  the  role  of  grain  geometry  (spherical 
vs.  cylindrical),  and  the  influence  of  wet/dry  cycling  compared  with  con¬ 
stant  wetting. 

The  governing  diffusion  equation  is 

— -V.(DVc)=  0  (3) 

where  c(x,t )  is  the  time-dependent  NG  concentration  within  the  grain  and 
D  is  diffusivity  (possibly  concentration  dependent). 

We  assume  that  the  NC  matrix  is  isotropic  with  respect  to  NG  diffusion. 

We  also  assume  that  the  initial  concentration  within  the  grain  is  uniform, 
c(x,o)=Co,  and  that  diffusion  of  NG  into  the  surrounding  water  is  rapid  rel¬ 
ative  to  diffusion  through  the  solid  NC  matrix.  This  latter  condition  is 
equivalent  to  specifying  c  =  o  on  the  grain  boundary  when  seeking  analyti¬ 
cal  solutions  (e.g.,  spherical  grain  immersed  in  water).  However,  for  stable 
numerical  solutions  and  to  examine  the  role  of  wet/dry  cycling,  we  found 
it  best  to  prescribe  the  flux,  q(x,t),  along  the  boundary: 

dry  grain  q  =  -DVc*  n  =  0  (4) 


wet  grain  q--DVc»n-h(c -cw)  (5) 

The  condition  set  up  in  eq  4  prevents  NG  from  leaving  the  grain  when  it  is 
dry,  although  NG  can  still  diffuse  within  the  grain.  We  may  simplify  the 
condition  set  up  in  eq  5  by  assuming  NG  concentration  in  the  water  is 
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small  relative  to  the  surface  concentration,  cw«c,  and  that  the  mass- 
transfer  coefficient,  h,  is  large  such  that  diffusion  within  the  grain  controls 
the  mass  flux  into  the  water. 

Crank  (1975)  and  Skelland  (1974)  present  analytical  solutions  for  un¬ 
steady,  one-dimensional  diffusion  through  spheres,  slabs,  and  cylinders 
with  constant  diffusivity,  D0,  uniform  initial  concentration,  c0,  and  con¬ 
stant  surface  concentration.  These  all  take  similar  form  as  infinite  series 
for  concentration  profiles  and  cumulative  mass  fluxes.  We  focus  on  the  lat¬ 
ter  because  they  equate  to  the  cumulative  dissolved  NG  mass  measured 
during  dissolution  experiments.  For  a  sphere  of  radius  a,  the  fractional 
mass  loss  to  time  t  is  given  by 

— -  =  1 --^y.-^-exp  (-D0n27r2t/a2)  (6) 

M »  Kn=\n 

where  M»  is  the  ultimate  mass  loss,  where  the  total  NG  mass  within  the 
grain  is  M0.  A  two-term  expansions  of  eg  6  converges  well  for  short  dimen¬ 
sionless  time,  t  =  Dot/ a2- 


-3  T 


(7) 


Truncation  error  in  eq  7  is  less  than  0.1%  for  r<  0.15  or  equivalently 
M1/M0  <  0.86,  which  exceeds  the  fractional  mass  loss  encountered  for  all 
grains  tested.  Equation  7  approximates  the  fractional  mass  loss  for  a 
spherical  grain  with  constant  diffusivity  immersed  continuously  in  water. 

To  investigate  the  roles  of  gain  shape,  concentration-dependent  diffusivity, 
and  wet/dry  cycling,  we  discretized  eq  3-5  using  the  finite  element  method 
(FEM)  for  the  spatial  variables  and  the  finite  difference  method  for  time. 
The  simple  geometries  reduce  the  problem  dimensions.  Pitt  (2011)  de¬ 
scribes  the  numerical  approach  and  its  validation;  we  summarize  the  key 
points  here. 

We  apply  to  eq  3  the  standard  Ritz-Galerkin  approximation  (e.g.,  Brenner 
and  Scott  1994)  and  arrive  at  a  system  of  ordinary  differential  equations  in 
time: 


[M] 


dC_ 

a 


+  [K(C)]C  +  h[Msl]C  =  0 


(8) 
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where 

[M]  =  mass  matrix 

[Man]  =  mass  matrix  for  the  boundary  terms 
[f<C(C)]  =  concentration-dependent  system  matrix 
C(f)  =  vector  of  time-dependent  nodal  unknowns. 

Note  we  have  assumed  that  cw  =  o  and  can  switch  between  no-flux  and 
flux  boundary  conditions  (eq  4, 5)  by  simply  setting  h  to  zero  or  some  large 
value. 

For  time  integration,  we  use  a  standard  second-order  Crank-Nicholson 
scheme: 


([M]  +  k[K(C"+l)]  +  kh"+\Mj)Cn+ 1  (g 

=  {[M]-k[K{Cn)]-kh"[Mj)Cn 

where  n  denotes  solution  vector  at  previous  time  step,  h  can  change  due  to 
wet/dry  cycling,  and  k  =  At/ 2  is  half  the  fixed  time  step. 

Equation  9  has  a  nonlinearity  in  \_KCn+1 ]  to  allow  for  concentration- 
dependent  diffusivity.  We  wanted  to  explore  this  possibility  because  early 
work  on  the  formation  of  double-base  propellants  revealed  that  NG  diffu¬ 
sion  into  NC  decreased  with  decreasing  concentration,  producing  a  sharp- 
edge  penetration  front  (Levy  1955).  Unfortunately,  we  have  not  found  in¬ 
formation  on  the  form  of  this  dependence,  so  as  a  first  approximation  use 
a  simple  linear  dependence: 


D(c)  =  D0- 
co 


(10) 


To  solve  for  the  required  [KCn+1],  we  use  a  semi-implicit  scheme  to  deter¬ 
mine  D(0=1): 


Z)(C"+1)wZ)(C")  +— (Cn -C”'1)  (11) 

co 

We  implemented  the  numerical  algorithm  in  C++  utilizing  the  deal.ii  dif¬ 
ferential  equation  solver  library  (Bangerth  et  al.  2007).  The  symmetry 
simplifications  allow  use  of  a  desktop  computer  to  run  the  simulations. 
Using  the  method  of  manufactured  solutions  (Roach  2002),  the  code 
shows  at  least  second-order  convergence  with  grid  refinement  (Pitt  2011). 
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FEM  predictions  also  agreed  within  0.1%  with  analytical  predictions  for  a 
spherical  grain  with  constant  diffusivity  (eq  7). 

3.4  Model  and  experiment  comparison 

We  focused  our  attention  on  modeling  the  experiments  that  used  unfired 
propellants  owing  to  the  complex  geometries  of  the  fired  grains.  The  batch 
tests  represent  the  simplest  case,  with  single  grains  in  a  stirred  bath  satis¬ 
fying  our  modeling  assumptions  of  cw  ~  o  and  high  surface  mass-transfer 
coefficient  h.  Figure  33  compares  data  from  the  batch  tests  with  analytical 
predictions  based  on  spherical  grains  with  constant  diffusivity  (eq  7). 

All  grains  have  mass  density  p~  1.6  g  cm"3  from  which  we  calculated 
equivalent  spherical  radii.  This  leaves  diffusivity  D0  as  the  only  tuning  pa¬ 
rameter,  and  we  fit  the  curves  to  the  data  for  the  first  3-4  measurements 
to  estimate  it  for  each  propellant  type.  Table  6  summarizes  the  parameters 
used.  Note  that  best-fit  diffusivities  are  orders-of-magnitude  lower  than 
that  of  NG  diffusing  into  water  at  similar  temperature  (£)0  ~  6x10  6  cm2 
s'1,  Wilke  1949).  This  satisfies  a  modeling  requirement  that  NG  mass 
transfer  from  a  grain’s  surface  is  high  relative  to  its  diffusion  rate  through 
the  NC  matrix. 


Figure  33.  Analytical  model  for  spherical  grains  and  constant 
diffusivity  compared  with  data  from  batch  tests  using  single, 
unfired  grains.  Analytical  curves  were  fit  to  initial  portion  of  the 
experimental  data  and  yield  different  diffusivities  for  each 
propellant  type. 
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Table  6.  Parameters  used  in  analytical  model,  eq  7  derived  for  spherical  grains  with  constant 
diffusivity,  using  radius  for  grain  of  equivalent  mass  ( r  =  1.6  g  cm-3).  To  simulate  drip  tests  on 
multiple  grains,  model  uses  average  grain  mass  and  diffusivity  obtained  from  best  fit  to 
single-grain  data. 


Propellant 

type 

Single  grain  batch  tests 

Multiple  grain  drip  tests 

Grain 

mass 

(mg) 

Radius 
a  (pm) 

Diffusivity* 
Do  (cm2/s) 

Grain 

mass 

(mg) 

Radius 
a  (pm) 

Diffusivity** 
Do  (cm2/s) 

Diffusivity*** 
Do  (cm2/s) 

9-mm 

0.080 

229 

3.2  x  10-14 

0.050 

195 

3.2  x  10-4 

3.2  x  10-14 

5.56-mm 

0.210 

315 

3.0  x  10-13 

0.157 

286 

3.0  x  10-13 

5.9  x  lO-i4 

7.62-mm 

0.210 

315 

4.5  x  10-13 

0.179 

299 

4.5  x  10-13 

1.0  x  10-13 

0.50-cal 

0.720 

475 

2.8  x  IO-12 

0.455 

408 

2.8  x  lO-i2 

1.2  x  10-13 

*  Best  fit  to  first  3-4  data  points. 

**Used  same  diffusivities  as  those  estimated  for  single  grains. 

***Best  fit  to  tail  slope. 

The  model  produces  reasonable  agreement  with  the  data,  particularly  giv¬ 
en  its  simplicity.  However,  Figure  33  reveals  two  important  concerns:  ex¬ 
cept  for  the  9-mm  propellant,  the  model  over-predicts  mass  loss  at  later 
times,  and  the  best-fit  diffusivities  span  two  orders-of-magnitude  although 
propellant  formulations  are  all  similar  (double-base,  NG  within  NC  ma¬ 
trix). 

Figure  34  compares  data  from  drip  tests  on  multiple  grains  with  the  same 
analytical  model.  We  retained  the  best-fit  diffusivities  determined  from 
the  batch  tests  but  used  the  average  mass  of  each  group  to  calculate  equiv¬ 
alent  spherical  radius.  Figure  34  reinforces  our  concerns  that  a  simple  dif¬ 
fusion  model  might  not  capture  key  physics  of  NG  dissolution  from  the 
propellants.  Except  for  the  0.50-cal.  results,  model  agreement  worsens, 
failing  to  track  even  short-time  NG  dissolution  from  the  grains.  Interest¬ 
ingly,  the  dissolution  rates  during  the  drip  tests  were  all  higher  than  those 
from  the  corresponding  batch  tests.  Given  that  mass-transfer  from  the  sur¬ 
face  into  the  water  is  high  in  both  cases,  the  net  diffusion  rates  were  some¬ 
how  higher  for  drip  tests  per  day  of  wetting. 
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Figure  34.  Analytical  model  for  spherical  grains  and  constant 
diffusivity  compared  with  data  from  drip  tests  using  multiple, 
unfired  grains.  The  analytical  curves  used  the  same  diffusivities 
determined  from  best  fit  to  single-grain  batch  tests. 

We  formulated  the  FEM  simulations  to  assess  whether  some  straightfor¬ 
ward  additional  physics  would  improve  model  predictions.  For  example, 
grain  geometries  are  more  cylindrical  than  spherical.  Initial  NG  loss  from 
grain  edges  should  boost  dissolution  rates  relative  to  long-term  rates.  Ad¬ 
ditionally,  if  diffusivity  decreases  with  concentration,  NG  loss  could  drop 
off  rapidly  as  it  is  depleted  from  the  outer  layers.  Lastly,  because  the  drip 
tests  were  stopped  each  weekend,  internal  diffusion  would  steepen  NG 
concentration  profiles,  leading  to  higher  dissolution  rates  when  the  tests 
resumed.  This  could  account  for  higher  drip-test  dissolution  rates  com¬ 
pared  with  batch  tests. 

We  limited  the  FEM  study  to  a  single  propellant  type— unfired  5.56  mm. 
The  model  parameters  were  the  same  as  the  analytical  study  (Table  6)  with 
the  addition  of  the  initial  concentration,  c0  =  0.158  g  cm"3,  and  the  meas¬ 
ured  dimensions  of  the  cylindrical  grain,  radius  and  height,  both  347  pm. 

Figure  35  compares  FEM  predictions  for  continuous  wetting  with  data 
from  the  batch  test.  Modeling  the  grain  as  a  cylinder  indeed  increases  dis¬ 
solution  rate,  owing  to  an  increase  in  surface  area.  However,  dimension¬ 
less  time  is  too  small  (r  <  0.01)  to  reveal  a  reduced  rate  after  depletion  of 
NG  along  the  grain  edges.  Rescaling  based  on  surface  area  would  essen¬ 
tially  collapse  the  cylinder  and  sphere  results. 
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Figure  35.  FEM  predictions  compared  with  data  for  5.56-mm 
unfired  grains  in  batch  test  (constant  wetting). 


Figure  36  does  reveal  the  expect  influence  of  concentration-dependent  dif- 
fusivity  as  modeled  via  eq  8:  dissolution  rates  decrease  much  more  rapidly 
with  time  than  for  constant  diffusivity  as  concentration  (and  hence  diffu- 
sivity)  decreases  in  the  outer  layers.  However,  the  predicted  curves  (sphere 
and  cylinder)  do  show  as  abrupt  a  roll-off  as  the  data  do,  suggesting  that 
the  concentration  dependence  might  be  higher  order  than  linear. 

As  noted,  the  FEM  simulations  can  mimic  the  cyclic  wetting  pattern  of  the 
drip  tests,  namely  5  wet  days  followed  by  2  dry  days.  Figure  36  shows  the 
results  along  with  the  data  from  the  drip  tests  on  5.56-mm  propellant 
grains.  The  data  show  a  stronger  roll-off  than  the  FEM  predictions,  even 
for  variable  diffusivity.  As  with  the  simulations  of  the  batch  test  (Fig.  35), 
the  results  suggest  a  stronger  dependence  of  diffusivity  on  concentration 
than  the  linear  dependence  modeled  here. 
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Time  (wet  days) 


Figure  36.  FEM  predictions  for  constant  and  variable  diffusivity 
with  5-day  wet/2-day  dry  cyclic  wetting  (red  curves)  compared 
with  data  for  5.56-mm  unfired  grains  in  drip  test  and  FEM 
predictions  for  variable  diffusivity  and  continuous  wetting  7- 
day  wet/O-day  dry  and  cyclic  1-day  wet/6-day  dry. 

Figure  36  also  shows  FEM  results  for  continuous  wetting  and  l-day  wet/6- 
day  dry  cyclic  wetting,  both  for  variable  diffusivity.  All  results  are  plotted 
against  the  number  of  wet  days  to  demonstrate  that  wet/dry  cycling  in¬ 
creases  the  dissolution  rate  per  wet  day  owing  to  steepening  of  concentra¬ 
tion  gradients  at  the  grain  boundaries.  This  effect  would  occur  in  the  field, 
as  rainfall  days  would  generally  be  less  frequent  than  dry  days. 

We  modeled  easiest  cases  first  to  determine  whether  molecular  diffusion 
could  account  for  observed  dissolution  of  NG  from  small  arms  propellants. 
The  analytical  model  using  constant  diffusivity  and  continuous  wetting  has 
only  diffusivity  as  a  tuning  parameter.  Results  show  roll-off  with  time  is 
not  as  abrupt  as  data  and  require  best-fit  diffusivities  that  vary  by  two  or¬ 
ders  of  magnitude. 

FEM  extension  to  more  realistic  cylindrical  grain  shape  did  not  improve 
agreement,  indicating  that  the  abrupt  roll-off  is  not  due  to  depletion  of  NG 
near  grain  edges  with  constant  diffusivity.  FEM  simulations  using  linear 
concentration-dependent  diffusivity  show  a  more  abrupt  roll-off  but  still 
not  as  abrupt  as  the  data. 
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The  FEM  simulations  of  wet/dry  cycling  demonstrate  the  expected  effect 
that  dissolution  rate  per  wet  day  increases.  Although  the  results  did  not 
substantially  improve  agreement  with  data  from  drip  tests  with  similar 
wet/dry  cycling,  this  effect  would  occur  for  in-field  dissolution  of  propel¬ 
lant  grains  caused  by  intermittent  rainfall. 

The  FEM  model  has  two  major  tuning  parameters:  the  scale  and  concen¬ 
tration  dependence  of  NG  diffusivity  through  NC.  We  resisted  the  urge  to 
fit  the  data  by  varying  these  independently  for  each  propellant  type.  Ra¬ 
ther,  we  implemented  a  series  of  simple  diffusion  models  to  understand 
the  effects  of  geometry,  concentration  dependence,  and  wet/dry  cycling  on 
dissolution  rate.  We  prefer  to  wait  for  additional  investigations  into  the 
physical  and  chemical  structure  of  these  propellant  grains  and  use  that  in¬ 
formation  to  guide  future  modeling. 

The  analytical  and  FEM  results  presented  above  show  that  if  we  fit  diffu- 
sivities  to  the  initial  portions  of  our  mass  loss  curves  we  cannot  account  for 
the  abrupt  decrease  in  mass  loss  rates,  even  using  concentration  depend¬ 
ent  diffusivities  and  different  grain  geometries.  However,  if  dissolution  of 
near  surface  NG  droplets  is  driving  initial  dissolution  rates  (as  suggested 
by  our  microscopic  results),  then  the  long-term  trends  in  the  data  (tail 
slopes)  may  represent  diffusion-controlled  dissolution.  The  diffusivities 
obtained  by  fitting  the  tail  slopes  of  the  drip  and  batch  tests  differ  by  less 
than  a  factor  of  four  among  the  propellants  tested  compared  with  two  or¬ 
der  of  magnitude  differences  seen  when  the  initial  slopes  are  fitted  (Table 
6).  Furthermore,  the  diffusivities  for  the  batch  tests  and  the  drip  tests  are 
the  same  for  each  propellant. 

3.5  Propellant-soil  interactions 

3.5.1  Batch  soil  studies 

We  performed  kinetic  batch  adsorption  studies  on  2,4-DNT,  NG,  and  NQ 
in  three  different  soils  to  determine  adsorption  and  transformation  rates 
for  these  compounds.  Detailed  information  on  the  soils  is  given  in  Appen¬ 
dix  B.  Treatments  without  soils  (blanks)  exhibited  no  change  in  their  con¬ 
centration  during  the  120-h  experiment  (Fig.  37).  In  the  presence  of  soil, 
both  adsorption  and  transformation  of  these  compounds  occurred  (Fig. 

37).  The  initial  decrease  in  concentration,  1st  time  point  at  1  h,  was  at¬ 
tributed  to  adsorption,  while  further,  generally  more  gradual,  decrease  in 
concentration  was  attributed  to  transformation  (Fig.  37).  When  the  con- 
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centration  of  energetic  compounds  is  expressed  as  a  natural  logarithm  and 
plotted  as  a  function  of  time,  the  slope  of  the  line  can  be  used  to  determine 
transformation  rate  ( k )  (Fig.  38). 


W> 

E 


Q 

<N 


bJO 

E 


<D 

u 

c 

o 

u 


txo 

E 


3 

2.5 
2 

1.5 
1 

0.5 

0 

2.5 

2 

1.5 

1 

0.5 

0 

2.5 

2 

1.5 

1 

0.5 

0 


20  40  60  80 

Time,  h 


100 


120 


Figure  37.  Concentrations  of  2,4-DNT  (top),  NG  (center),  and  NQ  (bottom) 
found  in  suspensions  of  Sassafras,  Plymouth,  and  Catlin  soils  as  a  function 
of  time.  Our  control  sample  that  contained  no  soil  is  also  plotted. 
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Figure  38.  The  natural  logarithm  of  2,4-DNT  (top),  NG  (center),  and  NQ 
(bottom)  in  suspensions  of  Sassafras,  Plymouth,  and  Catlin  soils 
concentrations  plotted  as  a  function  of  time.  Also  plotted  is  our  control  that 
contained  no  soil.  The  slope  of  regression  line  is  the  transformation  rate  ( k ) 

(Table  A2). 

Kinetic  batch  soil  experiments  showed  that  2,4-DNT,  NG,  and  NQ  had  dif¬ 
ferent  behaviors  when  in  contact  with  soil.  For  2,4-DNT  and  NG,  there  was 
an  initial  drop  and  a  steady  decrease  in  solution  concentration,  indicating 
adsorption  and  degradation  in  the  three  soils  tested  (Fig.  37  top  and  cen¬ 
ter).  This  resulted  in  measurable  adsorption  and  transformation  constants 
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(Table  7)  that  were  influenced  by  soil  type.  Nitroguanidine,  on  the  other 
hand,  showed  minimal  interaction  with  soils  and  appeared  not  to  have  de¬ 
graded  during  the  experiment  (Fig.  37  bottom). 


Table  7.  Fate  and  transport  parameters  for  propellant  constituents  in  soils  determined  in 
batch  adsorption  studies. 


Compound 

Soil 

C% 

Kd  cm3  g-1 

Koc 

log  Koc 

k  hour1 

tV2d 

R2 

NG 

Sassafras 

1.30 

0.26 

20 

1.30 

0.0193 

1.49 

0.86 

NG 

Plymouth 

1.72 

1.41 

82 

1.91 

0.0117 

2.48 

0.99 

NG 

Catlin 

3.75 

1.27 

34 

1.53 

0.0317 

0.91 

0.78 

NQ 

Sassafras 

1.30 

0.60 

46 

1.66 

-0.0006 

-48.14 

0.18 

NQ 

Plymouth 

1.72 

0.44 

26 

1.41 

0.0001 

458.43 

0.01 

NQ 

Catlin 

3.75 

0.24 

7 

0.82 

-0.0018 

-16.37 

0.53 

2,4-DNT 

Sassafras 

1.30 

2.34 

180 

2.26 

0.0076 

3.80 

0.95 

2,4-DNT 

Plymouth 

1.72 

5.06 

294 

2.47 

0.0152 

1.90 

0.91 

2,4-DNT 

Catlin 

3.75 

15.30 

408 

2.61 

0.0174 

1.66 

0.93 

^Transformation  constant  (k)  and  t1/2  values  for  NQ  are  not  significantly  different  from  zero. 


The  2,4-DNT  adsorption  coefficients  ranged  from  2.34  cm3  g“i  for  Sassa¬ 
fras  soil  to  15.3  cm3  g~i  for  high  organic  matter  Catlin  soil.  The  Plymouth 
soil  adsorption  coefficient  was  between  these  two  values  (5.02  cm3  g“i). 
Lower  2,4-DNT  coefficients  (0.28  and  3.3  cm3  g  1)  were  measured  for 
Plymouth  soil  using  column  studies  (Yamamoto  et  al.  2004;  Dontsova  et 
al.  2009).  The  difference  between  values  from  our  batch  studies  and  the 
column  studies  is  that  batch  tests  measure  irreversible  and  reversible  ad¬ 
sorption  together,  while  column  studies  measure  irreversible  adsorption 
with  transformation.  Because  2,4-DNT  transformation  products  experi¬ 
ence  irreversible  adsorption  (Pennington  et  al.  2003;  Thorn  et  al.  2008),  a 
systematic  bias  exists  between  these  two  methods. 

Adsorption  coefficients  determined  for  2,4-DNT  showed  a  strong  correla¬ 
tion  with  organic  carbon  (OC)  content  ( R 2  =  1.00)  (Fig.  39a),  a  trend  also 
true  for  Kds  published  in  the  literature  (Fig.  39b).  The  slope  of  our  regres¬ 
sion  function  (5.21)  was  similar  to  the  slope  obtained  for  published  data 
(4.84),  confirming  2,4-DNT’s  affinity  for  organic  matter.  Our  2,4-DNT 
transformation  rates  ranged  between  0.0076  and  0.0174  h_1  and  increased 
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with  soil  organic  matter  OM  content  (Table  7).  The  measured  transfor¬ 
mation  rate  for  2,4-DNT  in  Plymouth  soil  was  0.0152  h_1  within  the  range 
previously  determined  for  this  soil,  0.002-0.0238  h_1  (Dontsova  et  al. 
2009). 


Figure  39.  Adsoprtion  coefficints  ( K<* )  plotted  as  a  function  of  carbon  content  for:  (a)  the  soils  we  studied 
and  (b)  literature  values  (Table  A2). 

Nitroglycerin  adsorption  coefficients  ranged  from  0.26  to  1.41  cm3  g~i  for 
the  soils  studies  (Table  7).  These  are  lower  than  the  ones  determined  for 
2,4-DNT  but  significantly  different  from  zero.  Our  values  are  similar  to  lit¬ 
erature-reported  values  (Dontsova  et  al.  2007,  2008;  Table  Ai).  NG  ad¬ 
sorption  coefficients  did  not  exhibit  a  strong  relationship  with  organic 
matter  content  (Fig.  39a),  indicating  that  other  mechanisms  are  responsi¬ 
ble  for  adsorption.  NG  is  a  polar  molecule  (Winkler  1985)  and  may  form 
dipole-dipole  and  hydrogen  bonds  with  polar  moieties  in  the  soils.  Report¬ 
ed  Kow  values  for  NG  (1.62-1.77)  (Mirecki  et  al.  2006)  vary  depending  on 
the  method  used  to  estimate  them,  but  generally  indicate  a  preference  for 
non-polar  interactions.  Transformation/degradation  rates  determined  for 
NG  were  high,  with  half-life  of  only  0.91  days  for  Catlin  to  2.48  days  for 
Plymouth  soils  (Table  7).  Similar  to  adsorption  rates,  there  was  no  trend  in 
NG  transformation  rates  as  a  function  of  soil  OC  content. 

Nitroguanidine  did  not  show  the  steady  decrease  in  solution  concentra¬ 
tions  that  was  observed  for  2,4-DNT  and  NG.  For  all  three  soils,  the  slopes 
of  the  regressions  lines  in  plot  of  In  NQ  concentration  vs.  time  (transfor¬ 
mation  rates)  were  not  significantly  different  from  zero  (Fig.  38  bottom, 

Table  7)  illustrating  its  conservative  behavior.  This  finding  is  consistent 
with  observations  for  similar  soils  (Dontsova  et  al.  2007;  Table  A2).  NQ 
adsorption  coefficients  (Table  7,  Fig.  39a)  were  all  less  than  one  and  tend- 
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ed  to  increase  with  decreases  in  soil  organic  carbon  content  (R2= 0.92),  in¬ 
dicating  repulsive  rather  than  attractive  interactions  with  soil  carbon.  This 
is  likely  related  to  the  polar  nature  of  the  NQ  molecule  with  log  Kow  values 
between  -0.89  and  0.156  (Mirecki  et  al.  2006),  which  results  in  low  affini¬ 
ty  for  non-polar  organic  matter  in  the  soils. 

Based  on  high  adsorption  and  transformation  coefficients  determined  for 
Catlin  soil  in  batch  studies,  as  well  as  preliminary  column  experiments  us¬ 
ing  2,4-DNT,  2,6-DNT,  NG,  and  NQ  in  solution,  which  indicated  high  re¬ 
tention  of  all  measured  compounds  in  this  soil,  we  decided  to  not  use  Cat¬ 
lin  soil  in  the  column  experiments.  Instead,  we  focused  on  Plymouth  and 
Sassafras  soils,  which  provided  a  more  conservative  estimate  of  soil  reten¬ 
tion  for  2,4-DNT,  NG,  and  NQ  and  allowed  us  to  measure  their  transport 
through  the  soil  columns. 

3.5.2  Soil  column  studies 

We  completed  a  total  of  31  individual  column  dissolution  and  transport 
experiments  (Table  8).  We  used  two  soils,  Sassafras  loam  from  Aberdeen 
Proving  Ground  and  Plymouth  sandy  loam  from  Massachusetts  Military 
Reservation,  two  flow  rates,  0.01  and  0.02  mL  min  \  and  three  different 
propellant  compositions:  Ml,  WC  860  and  M31.  For  two  of  the  propel¬ 
lants,  Ml  and  WC  860,  we  also  studied  fired  residues  (Fig.  40).  The  main 
energetic  component  of  Ml  propellant  was  2,4-DNT,  while  WC  860  con¬ 
tained  NG,  and  M31  had  both  NG  and  NQ  (Table  9).  Each  experiment  was 
conducted  twice,  with  and  without  flow  interruption.  The  0.01  mL  min”1 
flow  rate  corresponded  to  a  0.55  cm  h”1  flux,  and  the  0.02  mL  min”1  flow 
rate  to  1.1  cm  h”1  flux.  The  water  fluxes  used  in  this  study  were  about  12 
and  25%  of  the  average  l-year  rainfall  intensity  (4  cm  h”1)  for  the  eastern 
U.S.  At  the  end  of  the  soil  column  tests  we  extracted  the  propellants  and 
the  soil. 
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Figure  40.  Unfired  propellants  and  residues  used  in  the  column  study:  (a)  unfired  WC  860;  (b) 
WC  860  residues;  (c)  unfired  Ml;  (d)  Ml  residues;  and  (e)  unfired  M31  propellant. 


Table  8.  Soil  column  experiments. 


Propellant 

Energetic 

compounds 

Flow  rate,  mL 
min-1 

Soil 

Set 

Number  of 
columns 

WC  860 

fi  red + unfired 

NG 

0.01 

Plymouth 

1 

4 

0.01 

Sassafras 

III 

4 

0.02 

Plymouth 

II 

4 

0.02 

Sassafras 

IV 

4 

M31 

unfired  only 

NG,  NQ 

0.01 

Plymouth 

III 

2 

0.01 

Sassafras 

III 

2 

0.02 

Plymouth 

IV 

2 

0.02 

Sassafras 

IV 

1 

Ml 

fi  red + unfired 

2,4-DNT 

0.01 

Plymouth 

V 

4 

2,6-DNT 

0.01 

Sassafras 

V 

4 

0.02 

Plymouth 

* 

4 

0.02 

Sassafras 

* 

4 

*No  detectable  2,4-DNT  and  2,6-DNT  in  outflow  for  0.01  mm  h1  flow  rate  so  higher  flow  rate  tests  were  cancelled. 
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Table  9.  Composition  of  propellants  used  in  column  experiments  (Defense 
Ammunition  Center  2006). 


Function 

Compound 

WC  860  (%) 

Ml  (%) 

M31  (%) 

Energetic 
plasticizer  and 
binder 

Nitrocellulose 

78.7  ±  6.5 

85  ±2 

20  ±  1.3 

Nitroglycerin 

9.5  ±  2.1 

19  ±1 

Dinitrotoluene 

10  ±2 

Oxidizer 

Nitroguanidine 

54.7  ±  1 

Plasticizer 

Dibutylphthalate 

8  ±2.8 

5  ±  1 

1.5  ±  0.1 

Stabilizer 

Diphenylamine 

1.1  ±  0.5 

1±0.1 

2-nitrodiphenylamine 

1.5 

Flash  reducer 

Potassium  nitrate 

0.8  ±1 

Sodium  sulfate 

0.5 

Calcium  carbonate 

1 

Lead  carbonate 

1  ±  0.3 

Potassium  sulfate 

Cryolite 

0.3±0.1 

Antistatic 

coating 

Graphite 

0.4 

3.5.2. 1  Column  outflow 

No  2,4-DNT  was  detected  in  the  Ml  o.oi  mL  min”1  column  outflow.  This  is 
reasonable  given  the  slow  2,4-DNT  dissolution  rates  from  Ml  propellants 
(Fig.  26,  27;  Dontsova  et  al.  2009),  and  the  strong  interaction  between  2,4- 
DNT  and  the  soil  matrix,  particularly  the  irreversible  adsorption  of  amino 
products  produced  by  reductive  transformation  (Thorn  et  al.  2008).  Be¬ 
cause  we  found  no  2,4-DNT  in  the  Ml  column  outflow  at  0.01  mL  min”1 
flow  rate,  we  did  not  run  the  Ml  experiments  at  0.02  mL  min  1  flow;  larger 
water  fluxes  always  resulted  in  smaller  outflow  concentrations  of  energet¬ 
ics  in  the  other  experiments. 

A  comparison  of  the  column  experiments  using  fired  and  unfired  WC  860 
showed  that  the  residues  had  lower  outflow  NG  concentrations  than  the 
unfired  propellants  (Fig.  41  and  42,  Table  10).  Part  of  the  reason  was  that 
we  used  a  smaller  amount  of  fired  residues  (0.016  g)  than  unfired  propel¬ 
lant  (0.1  g),  because  we  had  only  a  small  amount  of  the  fired  residues.  The 
fired  residues,  however,  also  had  less  dissolution  per  unit  mass  than  the 
unfired  particles  (Table  10).  The  ratio  between  dissolution  per  unit  mass  of 
NG  in  unfired  and  fired  particles  was  3.9  for  the  0.01  mm  h”1  flow  rate  and 
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1.9  for  the  0.02  mm  h_1  flow  rate.  We  think  that  partial  burning  of  the  NG 
in  the  residue  resulted  in  lower  total  concentration  of  NG  (Table  3)  and 
decreased  the  NG  concentration  in  the  rim  close  to  the  surface  of  the  resi¬ 
due  (e.g.,  Fig.  11).  These  results  are  similar  to  what  we  found  in  the  drip 
tests  where  the  unfired  propellants  lost  15-20%  and  their  residues  lost  3  to 
7%,  also  resulting  in  ~4  per  unit  mass  ratio  for  unfired  and  fired  propel¬ 
lants  for  the  equivalent  flow  rate. 


Figure  41.  Breakthrough  curves  for  continuous  water  flow  onto  unfired  and  fired 
0.50-cal.  (WC  860)  propellants.  The  solid  vertical  line  indicates  when  the  propellant 
was  removed  from  the  soil  surface.  Greater  outflow  concentrations  were  observed  in 
unfired  propellants  indicating  greater  dissolution  per  unit  mass  of  propellant. 
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Plymouth,  WC860  unfired,  IF,  0.01  mL  min-1 
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Plymouth,  WC860  fired,  IF,  0.01  mL  min  1 
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Figure  42.  Breakthrough  curves  for  interrupted  water  flow  onto  unfired  and  fired 
0.50-cal.  (WC  860)  propellants.  The  solid  vertical  line  at  ~13  pore  volumes  indicates 
when  the  flow  was  interrupted.  The  solid  vertical  line  at  ~30  pore  volumes  shows 
when  the  propellant  was  removed  from  the  soil  surface. 


The  M31  propellant,  which  contains  19.5%  NG  and  55%  NQ,  initially 
showed  more  NQ  in  the  outflow  than  NG  as  expected  given  its  higher  con¬ 
tent  in  the  M31  grain.  NQ,  however,  had  a  lower  steady-state  concentra¬ 
tion  than  NG  (Fig.  43),  suggesting  that  less  NQ  leaves  the  propellant  than 
NG.  The  ratio  of  NG  to  NQ  in  the  outflow  was  4.57  for  0.01  mm  tf1  flow 
rate  and  2.52  for  0.02  mm  tf1  flow  rate. 
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Figure  43.  Breakthrough  curves  for  interrupted  water  flow  onto  unfired  M31 
propellant.  The  second  plot  shows  a  detailed  view  of  the  initial  decrease  in 
concentration  for  NG  after  flow  interruption.  No  decrease  is  observed  for  the  NQ, 
indicating  it  does  not  degrade  or  transform  in  soil.  The  vertical  line  at  ~45  pore 
volumes  indicates  when  the  flow  was  interrupted.  The  solid  vertical  line  at  ~70  pore 
volumes  shows  when  the  propellant  was  removed  from  the  soil  surface. 


By  multiplying  the  degradation  rate  by  the  pore  volume  and  then  dividing 
by  the  flow  rate,  we  calculated  the  fraction  of  NG  that  degraded  as  the  NG- 
laden  water  moved  through  the  soil  column.  We  expected  and  saw  a  larger 
fraction  of  the  NG  degrade  for  the  lower,  o.oi  mL  min1,  flow  rate  tests  be¬ 
cause  it  took  longer  for  one  pore  volume  (PV)  to  get  through  the  soil.  We 
saw  this  effect  in  both  the  Sassafras  and  Plymouth  soils,  although  the  deg¬ 
radation  did  not  significantly  change  the  total  NG  in  solution  (Tables  10 
and  n). 

The  soil  column  tests  showed  that  NQ  experienced  limited  adsorption  in 
the  studied  soils,  while  NG  had  larger  adsorption  coefficients,  indicating 
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natural  attenuation  of  this  compound  in  soils  (Fig.  44),  results  similar  to 
our  batch  studies  (Fig.  37,  Table  7).  NG  experienced  retardation  due  to  ad¬ 
sorption,  as  demonstrated  by  delay  in  NG  breakthrough,  while  there  was 
little  delay  in  NQ  breakthrough  relative  to  the  conservative  tracer,  Br- 
(Fig.  44).  NG  breakthrough  also  happened  later  in  Plymouth  soil  than  in 
Sassafras  soil,  consistent  with  their  measured  adsorption  coefficients  (1.41 
cm3  g“i  for  Plymouth  and  0.26  cm3  g“i  for  Sassafras)  (Table  7). 


Sassafras,  M31  unfired,  CF,  0.02  mL  min-1 
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Figure  44.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.02  mL  min-1  rate  onto  unfired  M31  propellant  in  Sassafras  and  Plymouth  soils. 
Breakthrough  of  NG  happened  later  in  Plymouth  compared  to  Sassafras  soil 
indicating  greater  adsorption  of  NG. 
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p 


Soil 

w 

F 

CF/I 

F 

Set 

Col. 

Cum.  NG 
outflow 

m 

PV 

(mL) 

Flow  rate 

(mL  mirri) 

Fraction 

trans¬ 

formed 

Cum,  NG* 
outflow 

04) 

Gum. 
sola 
flew  PV 

WC860 

mass 

(4) 

NG 

(%) 

NG 

(Hfl 

Loss 

Less* 

<*) 

O.OlnnLmirr1- 

Sassafras 

u 

IF 

III 

cl 

606,1 

2,95 

0.0093 

0.1023 

675.2 

61,8 

0.1021 

9.7 

9904 

612 

6,82 

Sassafras 

LI 

OF 

III 

c2 

568.1 

277 

0.0087 

0.1023 

632.8 

498 

0.1009 

9.7 

0787 

5.80 

6.47 

Plymouth 

U 

IF 

1 

c3 

1360,6 

3.49 

0,0088 

00777 

1475.1 

14.7 

0.1000 

9.7 

9700 

14,03 

15.21 

Plymouth 

U 

CF 

1 

04 

854,2 

3.44 

0.0083 

0.0810 

929.4 

35.1 

0.1000 

9.7 

9700 

881 

9.58 

Average  ±  SE 

847  ±  IB  2 

92 8 ±194 

40±10 

9±2 

9±2 

Sassafras 

F 

IF 

III 

C3 

42.3 

2.87 

0.0091 

0.1012 

46.9 

59.5 

0.0165 

7.8 

1287 

3,27 

3.64 

Sassafras 

F 

CF 

III 

C4 

13  5 

288 

00090 

01026 

216 

505 

00167 

78 

1303 

150 

167 

Plymouth 

F 

IF 

1 

05 

M5 

a  55 

0.0093 

0.0744 

39.4 

335 

0.0186 

7.8 

1295 

2,82 

3.04 

Plymouth 

F 

CF 

1 

c6 

44.0 

a  57 

0,0087 

0.0797 

47,8 

345 

0.0166 

7,8 

1295 

3,40 

3.60 

Average  ±  SE 

35±6 

39±6 

45  ±6 

3±0.4 

3  ±0.5 

0.02  mL  min1 

Sassafras 

U 

IF 

IV 

Cl 

385.1 

3.03 

0.0171 

0.0571 

408.4 

465 

0.1000 

9.7 

9700 

3,97 

4.21 

Sassafras 

u 

CF 

IV 

c2 

4018 

301 

00180 

00538 

424.6 

514 

0,1003 

9,7 

9729 

4,13 

4.36 

Plymouth 

0 

IF 

II 

cl 

718.0 

296 

0.0181 

0.0319 

741.7 

41,2 

0.1004 

9.7 

9739 

737 

7.62 

P^mcuth 

u 

CF 

II 

c2 

687.2 

265 

00197 

00263 

705.7 

52,4 

0,1002 

97 

9719 

7.07 

7.26 

Average  ±  SE 

54S±39 

570±89 

48±3 

6±0.9 

6±0.9 

Sassafras 

F 

CF 

IV 

c3 

34.9 

276 

0.0177 

00502 

367 

55.0 

00166 

78 

1295 

269 

2.84 

Plymouth 

F 

CF 

II 

c3 

55.7 

2,84 

0.0183 

0.0302 

57.5 

393 

0.0157 

78 

1225 

455 

4.69 

Plymouth 

F 

IF 

II 

04 

500 

285 

0,0177 

00315 

516 

496 

00160 

78 

1248 

4.00 

4,13 

Average  ±  SE 

47  ±6 

49±6 

48±5 

4±0.5 

4±0,5 

U=Unfirea  F=  Fired:  SE  =  Standard  error:  CF-oontinuou&  flow;  I F=  interrupted  flow;  PV  =  pore  volume;  *=AcoounKtipg  for  transformation,  Tia  reformat  ion  rateffr)  used  to 
calculate  this  value  was  determined  in  batch  experiments  and  equaled  0=019 3  hfi  in  Sassafras  soil,  and  0,0117  hr1  in  Plymouth  soil  (Table  7) 


Table  10.  NG  loss  from  unfired  and  fired  WC860  propellant  tested  in  column 
experiments  at  0.01  and  0.02  mL  min-1. 


For  experiments  using  double  and  triple-based  propellants  we  observed  a 
set  of  common  features.  Both  NG  and  NQ  had  high  initial  dissolution  fol- 
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Set 
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Cum  IMG  or 
NQ  outflow 
£P0 

PV 

Flow  rate 
(mL 
min-1) 

Fraction 

trans¬ 

formed 

Cum, 

NG* 

outflow 

(Itf) 

Cum. 

aoln. 

flow 

PV 

NtVNQ 

NG*/ 

NQ 

M31  (g) 

NG 

or 

NQ 

(%) 

NS  Of 
NQ(jig) 

Loss 

(*} 

Loss* 

£%) 

0.01  ml  mirr1  NG 

Sassafras 

IF 

III 

C5 

7990 

zm 

0.0096 

0.0945 

8824 

63.3 

4.46 

4.92 

0.0918 

19 

U442 

45.81 

50.59 

Sassafras 

IF 

III 

efi 

6250 

2.85 

0.0097 

0.0949 

6905 

61,3 

423 

4.67 

0.0SS5 

19 

16815 

37,17 

41,07 

Plymouth 

IF 

III 

c7 

5169 

2,85 

0,0093 

0,0598 

5497 

52.0 

4,42 

4.70 

0,0884 

19 

16796 

30,77 

32.73 

Plymouth 

OF 

III 

C8 

7662 

2,84 

00096 

00579 

8133 

682 

877 

4.00 

0,0935 

19 

17765 

4113 

45-78 

Average  ±  SE 

6768  ±  653 

7340  ±  731 

6l±3 

4±0.2 

5±0.2 

39±3 

43±4 

0.01  mL  mirr1  NQ 

Sassafras 

IF 

III 

c5 

1793 

54.7 

50215 

3.57 

Sassafras 

IF 

III 

c6 

1477 

54.7 

48410 

3.05 

Plymouth 

IF 

III 

c7 

1169 

54  J 

48355 

242 

Plymouth 

CF 

III 

c8 

2032 

54.7 

51145 

3,97 

Average  ±  SE 

1613 ±1 S3 

3  ±0.3 

0.02  mL  mirr1  NG 

Sassafras 

IF 

IV 

cS 

2735 

2.77 

0.0197 

0.0451 

2864 

60.2 

275 

2.87 

0.0905 

19 

17195 

15.91 

16,66 

Sassafras 

CF 

IV 

c6 

3400 

2.70 

0,0183 

0,0474 

3569 

59,4 

2.39 

2,51 

0,0929 

19 

17651 

19,26 

20,22 

Plymouth 

IF 

IV 

c7 

2840 

2.43 

0.0182 

0.0261 

2916 

620 

204 

2.09 

0.0896 

19 

17024 

1668 

1713 

Plymouth 

CF 

IV 

c8 

2676 

2.59 

0.0188 

0.0269 

2750 

621 

254 

261 

0.0890 

19 

16910 

15,82 

16.26 

Average  ±  SE 

2913±  166 

3025±185 

61±0l7 

2±Q.l 

3±Q,2 

17  ±0,8 

18±Ql9 

0.02  mL  mirr1  NQ 

Sassafras 

CF 

IV 

c5 

996 

54.7 

49504 

201 

Sassafras 

CF 

(V 

c6 

1423 

54.7 

50816 

2.80 

Plymouth 

IF 

IV 

c7 

1392 

54.7 

49011 

2.84 

Plymouth 

CF 

IV 

cS 

1052 

54.7 

48683 

2,16 

Average  ±  SE 

12161111 

2,5±CL2 

SE  -  Standard  error;  CF-con  tin  nous  flow;  IF- interrupted  flew;  PV  -  pore  volume;  * -Accounting  fort  ran  sfo  rotation-  NG  transformation  rate  {k)  used  to  calculate  this  value 
was  determined  in  batch  experiments  and  equated  0,0193  h_1for  Sassafras  soil  and  0,0117  h  1  for  Pfyirrouth  soil.  Transformation  rate  for  NQ  was  not  significantly 
different  from  zero  in  both  soils  (Table  IS). 


Table  11.  NG  and  NQ  loss  from  unfired  M31  propellant  in  column  experiments 
conducted  at  0.01  and  0.02  mL  min-1  flow. 
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lowed  by  near  exponential  decrease  in  concentration  (Fig.  41,  42,  43,  and 
Appendix  C).  The  peak  concentrations  were  always  lower  than  the  solubili¬ 
ties  of  these  compounds  (Table  2).  For  experiments  where  the  flow  was 
interrupted  (a  period  of  24  h  of  no  flow),  NG  concentration  in  the  leachate 
decreased  after  flow  interruption.  We  think  that  NG  continued  to  trans¬ 
form,  and  that  the  increased  residence  time  during  no  flow  resulted  in 
lower  leachate  concentrations.  For  NQ,  there  was  no  decrease  when  the 
flow  was  stopped,  confirming  no  NQ  degradation  in  our  soils.  After  flow 
was  resumed  and  one  or  more  PV  passed  through  the  column,  however, 
both  NG  and  NQ  concentrations  increased  to  levels  similar  or  higher  than 
before  flow  interruption. 

The  energetic  dissolution  patterns  observed  from  propellants  provide  clues 
about  the  mechanisms  involved  in  dissolution.  Decreased  energetic  con¬ 
centrations  over  time  (as  well  as  increases  in  dissolution  following  flow  in¬ 
terruption)  indicated  that  dissolution  was  limited  by  the  rate  of  diffusion 
of  these  compounds  through  nitrocellulose  (Dontsova  et  al.  2009).  It  took 
longer  for  NG  or  NQ  to  reach  the  surface  as  the  distance  between  the  sur¬ 
face  of  the  particle  and  energetic  front  increased  over  time.  Interrupting 
the  flow  provided  additional  time  for  the  energetics  to  reach  the  surface 
and  be  dissolved  by  the  water.  This  idea  was  supported  by  the  time-  but 
not  flow-dependent  dissolution  rate  of  both  NG  and  NQ  in  M31  propellant 
that  was  seen  both  in  the  drip  and  the  soil  column  tests.  In  M31,  NG  con¬ 
centrations  in  the  outflow  were  higher  at  the  0.01  mL  min”1  flow  rate  than 
at  the  0.02  mL  min”1  (Fig.  45)  and  more  NG  mass  was  dissolved  (Table  11). 
As  the  same  cumulative  amount  of  water  was  supplied  in  both  flow  rate 
treatments,  about  2.4  times  more  NG  was  dissolved  at  0.01  mL  min”1  flow 
rate.  A  similar  trend  was  observed  for  NQ  but  the  differences  were  not  sta¬ 
tistically  significant  at  the  95%  confidence  level.  If  the  measured  outflow 
concentrations  (adjusted  for  transformation)  were  the  same,  independent 
of  the  flow  rate,  we  would  conclude  that  propellant  dissolution  was  limited 
by  removal  of  dissolution  products  with  the  flow,  as  observed  for  explo¬ 
sives  (Taylor  et  al.  2009b).  However,  because  concentrations  decreased 
with  increase  in  the  flow,  it  is  likely  that  dissolution  was  limited  by  the 
supply  of  NG  from  the  propellant. 

In  addition,  we  observed  that  the  percentage  of  NG  dissolved  was  lower  for 
the  WC  860  compared  to  M31.  For  similar  flow,  the  unfired  WC  860  lost 
5.9-9.596  NG,  while  the  M31  lost  17.6-42.596  (Table  11),  similar  to  what  ob¬ 
served  in  drip  dissolution  studies  (Fig.  27).  The  difference  in  mass  loss  be- 
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tween  the  two  propellants  cannot  be  explained  by  the  total  surface  area 
available  for  dissolution,  as  WC  860  are  small  particles  (Fig.  40)  and  have 
greater  surface  area  per  unit  mass.  However,  NG/NC  ratio,  the  parameter 
used  to  normalize  dissolution  in  drip  studies,  can  explain  the  differences. 
M31  has  equal  amounts  of  NG  and  NC,  each  at  about  20%  (Table  9),  while 
in  WC  860  has  about  10  NG%  and  90%  NC.  If  NG  is  being  adsorbed  and 
held  by  NC,  as  we  concluded  from  the  drip  studies,  then  an  increased  NG 
to  NC  ratio  would  result  in  greater  NG  dissolution,  as  observed  here. 
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Figure  45.  Breakthrough  curves  for  interrupted  (IF)  water  flow  onto  unfired  M31 
propellant  in  Sassafras  soil  at  two  different  flow  rates,  0.01  and  0.02  mL  min-1. 
Greater  outflow  concentrations  of  both  NG  and  NQ  were  observed  at  the  slower  flow 
rate. 
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3. 5.2.2  Mass  balance 

At  the  end  of  the  column  tests,  we  extracted  the  propellants  and  the  soils 
(Table  12).  We  recovered  about  55%  of  2,4-DNT  in  the  Ml  residues  and 
93%  from  the  Ml  unfired  propellants  (Table  12).  Differences  between  fired 
and  unfired  particles  were  similar  to  ones  observed  in  the  drip  tests  (Fig. 
26).  We  were  able  to  recover  only  small  fraction  of  the  original  2,4-DNT 
(up  to  2%  for  some  columns)  when  we  extracted  the  soil  (Table  12  and  Ta¬ 
ble  B2).  The  2,4-DNT  concentrations  were  highest  in  the  top  layers  of  the 
soil  and  decreased  with  depth  (Table  B2).  The  unfired  WC860  propellant 
lost  21%  of  its  original  NG  and  the  fired  residues  18%  (0.02  mL  min”1  flow 
rate).  Of  the  NG  originally  in  the  M31,  between  49  and  61%  was  not  in  the 
extracted  propellants  and  larger  losses  occurred  at  lower  flow  rates. 

When  we  added  the  mass  of  energetic  extracted  from  the  propellants  to 
the  mass  analyzed  in  solution,  we  did  not  recover  the  original  mass  of  en¬ 
ergetics.  This  was  true  even  when  transformation  in  the  soil  was  taken  into 
account  and  when  the  energetics  remaining  in  the  soil  were  included.  Alt¬ 
hough  soil  extractions  indicated  that  there  was  some  2,4  DNT  and  NG  in 
the  soil  even  if  these  were  absent  in  the  outflow  (Table  12,  Appendix  B), 
the  soil  concentrations  were  small  and  did  not  close  the  mass  balance.  No 
NQ  was  detected  in  the  extracted  soils. 

Less  2,4-DNT  was  recovered  in  the  fired  (54.2%)  than  in  the  unfired 
(93.3%)  Ml  propellants.  This  can  be  explained  by  non-linear  irreversible 
adsorption  of  2,4-DNT.  If  adsorption  decreases  with  increase  in  concen¬ 
tration,  fired  propellant  with  lower  dissolution  rates  would  result  in  more 
adsorption  and  smaller  2.4-DNT  concentrations.  Non-linear  adsorption 
can  happen  if  there  are  a  limited  number  of  sorption  sites  in  the  soil.  En¬ 
ergetic  compounds  can  have  liner  sorption  over  most  concentrations,  but 
exhibit  decreased  sorption  during  initial  peak  concentrations  in  propel¬ 
lants.  For  the  WC860  propellant,  the  total  NG  recovery  was  83-5±3.o% 
and  was  similar  for  unfired  and  fired  particles  and  for  both  flow  rates.  This 
was  within  the  range  of  what  was  observed  in  the  drip  studies  (86-112%) 
(Table  4). 

In  drip  studies,  NG  recoveries  for  M31  were  83±5%  (Table  5),  while  in  the 
soil  column  studies  they  were  81.6%  (0.01  mL  min”1  flow)  and  68.7%  (0.02 
mL  min”1  flow). 
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Table  12.  Enrgetic  compounds  recovered  in  leachate,  extracted  from  the  soil  and  extracted  from  propellants  at  end  of 
column  experiments:  data  for  2,4-DNT  in  the  Ml,  NG  in  the  WC860,  and  NG  the  unfired  M31  propellant. 


Soil 

U/F 

CF/ 

IF 

Set 

Col. 

%Orig. 

cone. 

%  in 

leach.* 

%in 

soil 

%Extracted 

propellant 

%Energetics 

recovered 

MlO.OlmLmirr1 

Sassafras 

u 

IF 

V 

cl 

9.7 

ND 

0.01 

87.63 

87.64 

Sassafras 

u 

CF 

V 

c2 

9.7 

ND 

ND 

92.78 

92.78 

Plymouth 

u 

IF 

V 

C5 

9.7 

ND 

0.07 

Plymouth 

u 

CF 

V 

C6 

9.7 

ND 

0.08 

99.38 

99.46 

AverageiSE 

93.3i3.4 

93.3i3.4 

Sassafras 

F 

IF 

V 

c3 

9 

ND 

ND 

54.44 

54.44 

Sassafras 

F 

CF 

V 

c4 

9 

ND 

ND 

52.00 

52.00 

Plymouth 

F 

IF 

V 

C7 

9 

ND 

ND 

57.78 

57.78 

Plymouth 

F 

CF 

V 

C8 

9 

ND 

2.18 

50.33 

52.51 

AverageiSE 

53.6il.6 

54.2il.3 

WC860  0.01  mL  mirr* 

Sassafras 

U 

IF 

III 

cl 

9.7 

6.82 

0.44 

70.21 

77.47 

Sassafras 

u 

CF 

III 

c2 

9.7 

6.47 

0.44 

77.01 

83.92 

Plymouth 

u 

IF 

1 

c3 

9.7 

15.21 

NA 

Plymouth 

u 

CF 

1 

c4 

9.7 

9.58 

NA 

AverageiSE 

9.5i2.0 

73.6i3.4 

80.7i3.2 

Sassafras 

F 

IF 

III 

c3 

7.8 

3.64 

5.89 

69.49 

79.02 

Sassafras 

F 

CF 

III 

c4 

7.8 

1.67 

2.65 

74.87 

79.19 

Plymouth 

F 

IF 

1 

c5 

7.8 

3.04 

NA 

Plymouth 

F 

CF 

1 

c6 

7.8 

3.69 

NA 

AverageiSE 

3.0i0.5 

72.2i2.7 

79.1i0.1 

WC860  0.02  mL  mirr* 

Sassafras 

U 

IF 

IV 

cl 

9.7 

4.21 

ND 

81.96 

86.17 

Sassafras 

U 

CF 

IV 

c2 

9.7 

4.36 

ND 

82.89 

87.25 

Plymouth 

u 

IF 

II 

cl 

9.7 

7.62 

ND 

67.53 

75.15 

Plymouth 

u 

CF 

II 

c2 

9.7 

7.26 

ND 

81.86 

89.12 

AverageiSE 

5.9i0.9 

78.6i3.7 

84.4i3.2 

Sassafras 

F 

CF 

IV 

c3 

7.8 

2.84 

ND 

78.97 

81.81 

Plymouth 

F 

CF 

II 

c3 

7.8 

4.69 

ND 

103.85 

108.54 

Plymouth 

F 

IF 

II 

c4 

7.8 

4.13 

2.69 

63.59 

70.41 

AverageiSE 

3.8i0.5 

82.1ill.7 

86.9ill.3 

M31 0.OlmLmirr1 

Sassafras 

U 

IF 

III 

c5 

19 

50.59 

ND 

31.89 

82.48 

Sassafras 

U 

IF 

III 

c6 

19 

41.07 

ND 

40.89 

81.96 

Plymouth 

u 

IF 

III 

c7 

19 

32.73 

0.16 

47.95 

80.84 

Plymouth 

u 

CF 

III 

c8 

19 

45.78 

0.17 

35.11 

81.06 

AverageiSE 

42.5i3.8 

39.0i3.5 

81.6i0.4 
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Soil 

U/F 

CF/ 

IF 

Set 

Col. 

%0rig. 

cone. 

%  in 

leach.* 

%in 

soil 

%Extracted 

propellant 

%Energetics 

recovered 

M31 0.02  mL  mirr1 

Sassafras 

u 

CF 

IV 

c5 

19 

16.66 

ND 

49.79 

66.45 

Sassafras 

u 

CF 

IV 

c6 

19 

20.22 

ND 

46.11 

66.33 

Plymouth 

u 

IF 

IV 

c7 

19 

17.13 

0.09 

58.95 

76.17 

Plymouth 

u 

CF 

IV 

c8 

19 

16.26 

0.10 

49.63 

65.99 

Average±SE 

17.6±0.9 

51.1±2.7 

68.7±2.5 

ACN=acetonitrile;  SE  =  Standard  error;  CF=continuous  flow;  IF=interrupted  flow;  PV  =  pore  volume;  *=Accounting  for  transformation, 
NA=not  analyzed,  ND=not  detected.  The  NG  transformation  rate  (k)  was  calculated  from  batch  experiments  and  equaled  0.0193  hr1 
for  Sassafras  soil  and  0.0117  hr1  for  Plymouth  soil  (Table  7). 


One  reason  for  incomplete  mass  balance  may  be  irreversible  adsorption  or 
degradation  in  the  soil  not  accounted  for  by  batch-determined  degradation 
and  adsorption  parameters.  However,  should  adsorption  or  degradation 
contribute,  they  would  influence  concentrations  at  slower  flow  rates  and 
we  see  the  opposite. 

3. 5.2. 3  Risk  to  groundwater 

Because  2,4-DNT  was  not  detected  in  the  column  outflow,  we  think  it  un¬ 
likely  that  2,4  DNT  will  reach  groundwater  unless  the  source  concentra¬ 
tions  are  high  (1-10  mg/L)  and  the  depth  to  groundwater  small,  less  than 
3  m.  This  agrees  with  studies  of  training  ranges  where  no  2,4-DNT  was 
found  in  groundwater  even  if  it  was  detected  in  firing  point  soils  (Clausen 
et  al.  2004).  NG  experienced  both  adsorption  and  transformation  in  the 
soils,  resulting  in  retardation  of  the  breakthrough  curve  and  decrease  in 
eluted  concentrations.  Given  the  rate  of  transformation  over  short  flow 
paths,  we  think  that  the  majority  of  NG  will  be  transformed  before  reach¬ 
ing  groundwater.  As  NQ  does  not  significantly  degrade  in  the  soil,  it  is  like¬ 
ly  to  persist  and  eventually  reach  groundwater.  Fortunately,  it  appears  to 
be  less  toxic  than  2,4-DNT  or  NG  (Nipper  et  al.  2009).  Triple-base  propel¬ 
lants  are  also  not  as  widely  used  as  double-based  propellants  and  the  field 
test  we  witnessed  deposited  no  measurable  residues.  Because  NQ  is  not 
one  of  the  analytes  quantified  by  Method  8330B,  it  has  not  been  routinely 
sought  in  the  groundwater.  We  know  of  no  instances  of  NQ  in  groundwa¬ 
ter. 
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3.6  HYDRUS  modeling 

Computational  models  such  as  HYDRUS-iD,  a  widely  used,  one¬ 
dimensional  water  flow  and  solute  transport  code,  use  transport  parame¬ 
ters,  such  as  the  adsorption  coefficients,  Kd,  and  the  transformation  rates, 
k,  to  predict  transport  of  solutes  through  soils  and  into  groundwater.  In 
inverse  mode,  HYDRUS-iD  fits  the  transport  equations  to  the  experi¬ 
mental  breakthrough  curves  like  ones  obtained  in  our  column  experiments 
to  determine  values  for  the  transport  parameters  (Simunek  et  al.  2005). 

HYDRUS-iD  has  been  modified  by  Dr.  Simunek,  its  developer,  to  also  es¬ 
timate  dissolution  rate  of  energetics  in  the  solid  form,  such  as  explosive 
particles  (Dontsova  et  al.  2006)  or  propellant  grains  (Dontsova  et  al. 

2008,  2009).  To  predict  the  dissolution  of  2,4-DNT  from  Ml  propellant, 
HYDRUS-iD  used  an  empirical  model  that  assumed  exponential  decrease 
in  energetic  dissolution  with  time.  This  resulted  in  very  good  agreement 
with  experimental  dissolution  curves  for  Ml;  however,  it  did  not  work  as 
well  for  NG  and  NQ  dissolution  from  M9  and  M30  propellants  (Dontsova 
et  al.  2008).  Observed  dissolution  rates  for  NG  and  NQ  did  not  decrease  to 
zero  but  approached  a  steady  rate  that  continued  with  little  change.  There¬ 
fore,  in  the  current  study,  the  HYDRUS  was  expanded  to  include  a  steady- 
state  dissolution  rate.  This  version  was  used  to  analyze  all  the  break¬ 
through  curves  for  the  soil  column  studies.  However,  even  with  the  chang¬ 
es  implemented,  the  model  did  not  explicitly  address  the  processes  that 
cause  the  exponential  decrease  in  dissolved  energetics. 

To  describe  dissolution  more  mechanistically,  and  to  link  dissolution  ex¬ 
periments  conducted  with  and  without  soil,  we  also  performed  2-D  simu¬ 
lations  of  NG  and  NQ  dissolution  from  M31  propellant  in  the  drip  studies 
and  created  a  framework  for  modeling  both  dissolution  and  transport  of 
energetics  in  soil. 

3.6.1  1-D  HYDRUS  studies  of  propellant  dissolution  and  transport  in  the 
soil  columns 

3.6. 1.1  Dissolution 

For  propellants,  the  NG  and  NQ  dissolution  rate,  y  [pg  g_1  h_1]  was  defined 
as  follows: 


Yw  =  max  ,  ai;  ) 


(12) 
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Where 

4  =  initial  dissolution  rate  [jag  L_1  h_1] 

X  =  decay  constant  [h1] 

a  =  dimensionless  constant  that  defines  the  minimum  or  steady- 
state  dissolution  rate  (a£)  as  a  fraction  of  the  initial 
(maximum)  dissolution  rate  (£) 
t  =  time  [h]. 

If  a  is  larger  than  1,  a  constant  dissolution  rate  is  used.  When  presenting 
the  results,  dissolution  rate  was  recalculated  relative  to  mass  of  propellant 
added  to  the  column  [pg  g_1  h-1]. 

3. 6.1.2  Transport 

The  transport  of  the  Br  1  tracer  can  be  described  using  the  convection- 
dispersion  equation  for  constant  water  content,  flux  density,  and  disper¬ 
sion  coefficient: 


0 


dc 

~dt 


(13) 


Where 

c 

q 

e 

z 

t 

D 


solute  concentration  (pg  cm“3) 
convective  flux  (cm  h  x) 

water  content  (cm3  cm“3)  equal  to  porosity  for  saturated 

experiments 

spatial  coordinate  (cm) 

time  (h) 

dispersion  coefficient  (cm2  h_1)  equal  to  the  product  of  the 
longitudinal  dispersivity,  A  (cm),  and  q  divided  by  6. 


3.6.1.3  Equilibrium  sorption  model  (with  decay) 

The  full  model  for  soil  transport  implemented  in  HYDRUS  includes  two- 
site  sorption  (Selim  et  al.  1976;  van  Genuchten  and  Wagenet  1989)  and 
non-equilibrium  adsorption-desorption  reactions.  In  the  simulations  pre¬ 
sented  in  this  report,  we  used  only  equilibrium  sorption.  Here,  sorption,  s 
(pg  g_1)  is  assumed  to  be  instantaneous 


s  =  Kdc 


(14) 
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where  Kd  is  the  linear  distribution  coefficient  (cm3  g1),  and  c  is  the  solute 
concentration  (jag  cm~3).  The  governing  transport  equation  is  then  as  fol¬ 
lows  for  constant  water  content,  flux  density,  dispersion  coefficient,  bulk 
density,  and  distribution  coefficient: 

?+?  =  iK'f )-f-n»f>c-ws  +  Ywe  (15) 


where 

p  =  bulk  density  (g  cm  3) 

pw  =  first-order  rate  coefficient  in  the  liquid  phase  (IT1) 

y  =  dissolution  rate  as  described  in  Eq.  12  (pg  L_1  h_1). 

First-order  rate  coefficients  in  liquid  and  solid  phase  were  set  to  be  equal. 
3.6.1.4  Parameter  estimation 

The  numerical  analysis  of  experimental  data  was  carried  out  as  follows. 
First,  the  conservative  tracer  breakthrough  curves  were  used  to  estimate 
the  longitudinal  dispersivity,  A  (cm).  Then  the  longitudinal  dispersivity 
was  fixed  at  a  value  determined  for  the  tracer  and  the  chemical  equilibri¬ 
um  model  was  then  used  to  analyze  energetic  breakthrough  curves.  While 
it  has  been  shown  before  that  chemical  non-equilibrium  may  play  a  role  in 
adsorption  and  transport  of  explosives  in  soils  at  similar  fluxes  and  for  the 
same  soils  as  used  in  this  study  (Dontsova  et  al.  2006),  kinetic  sorption 
has  been  not  as  important  for  NG  and  NQ  (Dontsova  et  al.  2008).  The  fol¬ 
lowing  parameters  were  estimated  for  energetics:  adsorption  coefficient, 
Kd  (g  1  cm3),  and  three  parameters  that  describe  dissolution:  initial  disso¬ 
lution  rate  £  [pg  L1  h-1],  decay  constant  %  [h1],  and  constant,  a,  that  mul¬ 
tiplied  by  Ogives  the  steady-state  dissolution  rate  [pg  L1  h1].  Transfor¬ 
mation  rate  coefficients  (pw,  same  as  k)  determined  in  batch  studies  were 
used  to  account  for  transformation  of  solute  in  the  soil. 

The  R 2  values  and  confidence  intervals  for  fitted  parameters  were  obtained 
by  analyzing  the  correspondence  between  measured  and  fitted  break¬ 
through  concentrations  and  the  behavior  of  the  objective  function  around 
its  minimum,  respectively  (Simunek  and  Hopmans  2002).  Parameter  es¬ 
timates  were  considered  significant  if  they  were  different  from  zero  (confi¬ 
dence  intervals  did  not  intersect  with  zero).  For  comparison  between 
treatments,  differences  were  considered  significant  when  greater  than  the 
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sum  of  standard  errors  of  the  means  multiplied  by  1.96  (for  95%  probabil¬ 
ity). 


3.6.2  2-D  HYDRUS  simulations  of  energetic  dissolution  from  M31 

propellants 

To  model  dissolution  of  energetics  from  propellant  particles  (a  necessary 
step  needed  to  compare  experiments  that  were  conducted  with  and  with¬ 
out  soil  and  to  be  able  to  model  dissolution  in  soil  columns  mechanistical¬ 
ly),  we  performed  HYDRUS-2D  simulations  for  dissolution  curves  from 
the  drip  studies  of  M31  propellants  (Fig.  30).  The  two  kinetic  sorption  sites 
model  was  applied  to  the  data.  This  model  was  similar  to  the  one  used  for 
soil  transport,  but  looked  at  the  condition  where  there  is  no  flow  within 
the  particle  and  the  only  way  for  energetics  to  reach  the  surface  is  through 
diffusion.  The  model  was  run  in  axi-symmetrical  2-D;  equations  here  are 
shown  in  l-D  to  simplify  them: 


which  leads  to 


Si  =  Kdc 


P^  =  dka2c-kd2ps2 


P^§  =  dka3c-kd3ps3 


ddRc 

St 


R  =  1  + 


pKD 


(16) 

(17) 

(18) 

(19) 
(16) 

(20) 
(21) 


There  was  an  equilibrium  phase  (characterized  by  the  retardation  factor  R) 
and  two  kinetic  sorption  sites,  one  with  faster  release  and  one  with  slower 
release  (basically  irreversible).  For  the  kinetic  sites  with  rates  of  attach¬ 
ment  and  detachment,  the  equilibrium  state  is  characterized  as: 
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Tf  _  @kg 

ad*  — 

pkd 


(22) 


This  model  assumes  that  dissolution  of  energetics  from  propellant  parti¬ 
cles  is  influenced  by  diffusion  and  equilibrium,  kinetic  and  irreversible  ad¬ 
sorption  of  NG  to  the  NC  matrix.  The  model  assumes  presence  of  the  water 
phase  inside  the  propellant  particle. 

3.6.3  Axi-symmetrical  2-D  HYDRUS  simulations  of  propellant  in  soil 
column 


Inflow 

A 


Soil  column,  L 


▼ 

BTC 


Figure  46.  HYDRUS-2D  axi-symmetrical  domain  showing 
modeled  processes  for  dissolution  of  energetic  from 
propellant  particles  and  their  transport  in  soil.  The  image 
is  not  done  to  scale  to  show  detail. 


The  two  models  discussed  above,  one  for  energetic  movement  through 
propellant  particle  interior  to  the  surface  and  another  for  its  reactive 
transport  in  the  soil,  were  then  combined  to  characterize  energetic  dissolu¬ 
tion  and  transport  in  the  column.  We  set  up  the  experimental  domain 
(propellant  particle  on  the  surface  of  soil  column)  in  2-D  to  simulate  axi- 
symmetrical  (radial)  3-D.  We  multiplied  each  element  by  2ttxR  (circum¬ 
ference)  (i.e.,  distance  from  the  center  of  symmetry)  to  achieve  axi- 
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symmetrical  3-D  and  relate  this  to  propellant  flux  in  the  column  outlet. 
HYDRUS-2D  simultaneously  solved  the  transport  within  the  grain  (on  the 
scale  of  1  mm  or  less,  R )  and  then  within  the  column  (on  the  scale  of  about 
10  cm,  L)  (Fig.  46).  We  tested  this  method  using  one  sample  dataset. 

3.7  HYDRUS  model  and  experiment  comparisons 

3.7.1  Propellant  dissolution  and  transport  in  the  soil  columns 

HYDRUS-iD  simulations  using  a  steady-state  dissolution  rate  (eq  12)  re¬ 
sulted  in  good  agreement  between  observed  and  fitted  breakthrough 
curves  (Fig.  47,  Table  13).  In  the  majority  of  the  cases,  parameter  a,  which 
defines  steady-state  dissolution  rate  as  a  fraction  of  original  dissolution 
rate,  was  significantly  contributing  to  the  objective  function.  Generally, 
better  fit  was  produced  for  NG  than  for  NQ  (Fig.  47)  and  determined  dis¬ 
solution  and  transport  parameters  for  NQ  were  largely  non-significant. 
The  reason  the  new  model  did  not  accurately  describe  NQ  breakthrough 
curves  is  that  NQ  generally  had  a  high  initial  peak  followed  by  a  fast  de¬ 
crease  and  then  a  slower  decrease  in  dissolution  rate  over  a  long  time  peri¬ 
od.  An  exponential  decrease  in  dissolution  rate  that  we  modeled  could  de¬ 
scribe  the  sharp  peak,  but  not  gradual  decrease  at  the  same  time,  so 
HYDRUS  fitted  a  steady-state  dissolution  rate  to  the  second  half  of  the 
breakthrough,  when  concentration  in  reality  changed  over  time.  There 
were  some  potential  issues  with  modeling  NG  dissolution,  as  well.  In  par¬ 
ticular,  the  behavior  of  the  actual  and  modeled  breakthrough  cures  during 
flow  interruption  indicated  that  transformation  rates  determined  in  the 
batch  experiments  might  have  underestimated  transformation  rates  of  NG 
in  soils.  If  this  is  true,  it  would  be  consistent  with  irreversible  adsorption 
disproportionally  contributing  to  transformation  in  columns  studies  and 
would  also  help  close  the  mass  balance.  However,  there  was  not  enough 
evidence  to  support  deriving  transformation  rates  from  interrupted  flow 
experiments.  While  increasing  the  transformation  rate  improved  fit  for 
several  simulations,  this  was  not  consistent  across  experiments,  and  for 
some  it  decreased  the  R2. 
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Figure  47.  Observed  and  HYDRUS-1D  generated  breakthrough  curves  for  NG  (a)  and  NQ  (b 
and  c).  (a)  Unfired  WC  860  propellant  in  Plymouth  soil,  interrupted  flow  experiment,  0.02  mL 
min-1  flow  rate,  set  II,  ci,  obtained  parameters:  x  =0.026±0.001;  £=189±  6  pg  g_1  h_1;  a  = 
0.136±0.017;  Kd  =  0.79±0.02  cm3  g_1;  R2= 0.99.  (b)  and  (c)  (detail)  unfired  M31  propellant  in 
Sassafras  soil  at  continuous  0.02  mL  min-1  flow  rate,  set  IV,  C6,  x  =0.56±0.06;  £, 
=3376±233  pg  g"1  h"1;  a  =0.018±0.004;  Kd  =  0.16±0.01  cm3  g"1;  R2  =0.98. 
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dissolution  rate,  ngg1  h1;  a  -  constant,  dimensionless;  a£-  minimum  (steady- state)  dissolution  rate,  pgg  1  Ir1.  NG  transformation  rate  (fr)  was  determined  in  batch 
experiments  and  equaled  0,0193  h  1  in  Sassafras  soil,  and  0.0117  h  1  in  Plymouth  soil  (Table  7),  Shaded  parameter  estimates  were  not  statistically  significant. 


io 


Table  13.  Solute  transport  parameters  obtained  by  HYDRUS-1D  for  column  saturated  flow 
experiments  involving  Br  1  and  WC  860  and  M31  propellants  in  Plymouth  and  Sassafras 
soils.  X  was  estimated  from  Br1  breakthrough  and  characterizes  water  flow  in  the  column, 
while  £ ,  a,  and  Kd  were  determined  from  NG  breakthrough.  NQ  transport  parameters  were 
not  significantly  different  from  zero  and  are  not  presented. 
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Much  higher  NG  dissolution  rates,  both  maximum  (1216  ±  192  jag  g-1  h_1 
vs.  226  ±  35  pg  g_1  h_1)  and  steady-state  (83  ±  15  jag  g_1  h_1  vs.  17  ±  6  pg  g_1 
h_1)  were  calculated  for  M31  than  for  WC860  propellant,  in  agreement 
with  drip  studies.  In  M31  propellant,  steady  state  rates  were  higher  for 
higher  flow  velocity  (99  vs.  67  pg  g_1  h_1)  but  the  differences  were  not  sig¬ 
nificant  at  95%.  Flow  rate  did  not  significantly  influence  WC860  dissolu¬ 
tion.  Soil  did  not  affect  dissolution  rates  for  either  propellant,  as  expected. 
Determined  adsorption  coefficients  were  0.92  ±  0.02  cm3  g-i  for  Plymouth 
and  0.35  ±  0.05  cm3  g-i  for  Sassafras,  in  general  agreement  with  batch 
numbers  (1.41  cm3  g-i  and  0.26  cm3  g-i  for  Plymouth  and  Sassafras,  re¬ 
spectively).  High  initial  concentrations  of  energetics  dissolving  from  pro¬ 
pellant  particles  can  result  in  non-linear  adsorption  and  lower  measured 
Kds. 

Initial  NG  dissolution  rate,  calculated  for  the  unfired  WC860  propellant 
(226  pg  g_1  h_1),  was  lower  than  one  previously  calculated  for  the  similarly 
sized  M9  (1583  pg  g_1  h_1)  (Dontsova  et  al.  2008).  The  difference  can  be 
explained  by  the  lack  of  graphite  on  the  surface  of  M9  but  more  likely  by 
the  much  higher  NG  to  NC  ratio  in  this  propellant  (0.667  for  M9  and  0.125 
for  WC860)  further  supporting  conclusions  made  in  drip  studies  about  the 
importance  of  how  NG  adsorbs  to  NC  and  the  NG/NC  ratio. 

At  the  same  time,  M31  had  similar  NG  to  NC  ratios  to  previously  studied 
M30  (1  and  0.8,  respectively),  but  much  higher  initial  dissolution  (1216  pg 
g-i  foi  jn  M31  and  60  in  pg  g_1  h1  in  M30).  A  lower  dissolution  rate  of  NG 
from  M30  propellant  may  be  explained  by  the  very  large  size  of  M30  parti¬ 
cles  (2.2  cm  in  length  and  1.05  cm  in  diameter)  and  their  smaller  surface 
area  to  mass  ratio. 

3.7.2  2-D  HYDRUS  simulations  of  energetic  dissolution  from  M31 

propellants 

HYDRUS-2D  simulations  of  NG  and  NQ  dissolution  from  M31  propellant 
in  drip  studies  indicated  that  diffusion  only  did  not  adequately  describe 
experimental  results  (Fig.  48a).  This  agrees  with  diffusion  modeling  of  NG 
dissolution  from  small  arms  propellants  done  in  this  project.  To  improve 
the  fit,  we  included  an  interaction  term  between  NG,  NQ,  and  NC  matrix 
in  our  simulations.  Adding  equilibrium  absorption  (eq  17)  did  not  improve 
the  agreement  between  experiment  and  simulation  (not  shown).  However, 
adding  kinetic  adsorption  (eq  18  and  19)  significantly  improved  the  fit 
(Fig.  48b). 
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Figure  48.  Observed  and  HYDRUS-1D  generated  cumulative  dissolution  curves  for  NG  from 
M31  propellant  particle  in  drip  studies  without  soil,  (a)  Diffusion  only:  D  =  1.81x10  10  ± 
4.82xl0"5,  cm2  s1;  R2  =  0.91;  (b)  Diffusion  and  adsorption  on  two  kinetic  sites  in  NO:  D  = 
6.63xl0"9  ±  1.06  xlO"9  cm2  s_1;  adsorption  rate  for  fast  sites,  ka i  =  0.1775  ±  0.0172  day-1; 
desorption  rate  for  fast  sites,  kdi  =  0.1163  ±  0.0110  day -1;  adsorption  rate  for  slow  sites  ka 2 
=  0.0050  ±  0.0008  day-1;  desorption  rate  for  slow  sites  kd 2  =  0.0024  ±  0.0003  day1;  R2  = 
1.00. 

For  NG,  adsorption  on  two  types  of  kinetic  sites  (one  slow  and  one  fast) 
better  described  the  data  than  one  adsorption  site  only,  but  for  NQ  one  ki¬ 
netic  site  was  sufficient.  In  fact,  adsorption  and  desorption  rates  for  a  se¬ 
cond  site  were  not  significant  (Table  14),  and  removing  this  second  kinetic 
site  did  not  change  other  parameters.  For  NG  one  kinetic  site  resulted  in  a 
good  fit,  though  not  as  good  as  two  sites.  The  diffusion  coefficient,  howev¬ 
er,  was  an  order  of  magnitude  lower  if  only  one  sorption  site  was  modeled. 
This  indicates  that  both  sites  were  contributing  to  NG  adsorption.  The  fact 
that  kinetic  sorption  is  needed  to  describe  dissolution  of  energetics  from 
propellant  particles  indicates  that  slow  release  of  solute  from  the  sorption 
sites  slows  NG  and  NQ  dissolution. 
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Table  14.  HYDRUS  2D-determined  parameters  for  dissolution  of  NG  and  NQ  from  M31  particles  in 
drip  studies  in  the  absence  of  soil:  diffusion  coefficient,  adsorption  and  desorption  rates  for  two 
different  kinetic  adsorption  sites,  adsorption  coefficient  calculated  for  these  sites. 


Parameter 

NG 

NQ 

Average 

SE  (n=8) 

Average 

SE  (n=8) 

Two  kinetic  adsorption  sites 

D,  cm2  s-1 

2.09x10-8 

4.39x10-9 

1.78x10-9 

3.74x10-10 

ka  1,  day-1 

0.7519 

0.1261 

0.0203 

0.0200 

kdi,  day-1 

0.2953 

0.0469 

0.2511 

0.0193 

ka2,  day-1 

0.0051 

0.0007 

0.0066 

0.0005 

kd2,  day-1 

0.0013 

0.0003 

0.0005 

0.0001 

Kdi*,  cm3  g_1 

0.88 

0.16 

0.04 

0.04 

Kd2*,  cm3  g-1 

1.42 

0.17 

4.33 

0.62 

One  kinetic  adsorption  site 

D,  cm2  s-1 

2.55x10-9 

3.24xl0-10 

1.89x10-9 

3.31x10-!° 

ka2,  day-1 

0.0014 

0.0006 

0.0072 

0.0014 

kd2,  cm3  g-1 

0.0011 

0.0003 

0.0005 

0.0001 

Kd2*,  cm3  g-1 

0.37 

0.04 

4.07 

0.37 

SE  =  standard  error  of  the  mean;  Shaded  numbers  were  not  statistically  significant.  D  =  diffusion  coefficient, 
cm2  s-1;  kai=  adsorption  rate  for  fast  sites,  day-1;  kdi=  desorption  rate  for  fast  sites,  day-1;  ka 2  =  adsorption  rate 
for  slow  sites,  day-1;  kdi=  desorption  rate  for  slow  sites,  day-1;  Koi*=adsorption  coefficient  calculated  for  fast 
sites,  cm3  g-1;  Ko2*=adsorption  coefficient  calculated  for  slow  sites,  cm3  g-1. 


If  no  adsorption  was  modeled,  fitted  diffusion  coefficients  were  lower  than 
determined  in  simulations  that  included  kinetic  adsorption,  as  slow  re¬ 
lease  was  attributed  to  diffusion.  If  we  included  sorption,  slow  release  was 
explained  partially  by  slow  desorption.  The  diffusion  coefficients 
(2.09x10  8  ±  4.39x10  9  cm2  s  1  for  NG  and  1.78x10  9  ±  3.74x10  10  cm2  s  1 
for  NQ)  were  higher  than  found  for  small  arms  propellants  in  this  study 
but  in  general  agreement  with  numbers  reported  by  Levy  (1955)  for  nitro¬ 
glycerin  in  cellulose  acetate,  5.2x10  9  cm2  s”1.  Lack  of  contribution  from 
fast  sites  to  NQ  adsorption  was  also  confirmed  by  the  low  adsorption  coef¬ 
ficient,  Kdi*,  for  these  sites,  0.04  cm3  g"i  (Table  14).  The  fact  that 
HYDRUS-2D  simulations  require  adsorption  term  to  describe  energetic 
dissolution  from  propellants,  supports  observations  by  Yazici  and  Kalyon 
(1998)  that  NG  in  propellant  particles  is  adsorbed  and  held  by  NC  prevent¬ 
ing  its  movement  within  the  particle. 
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3.7.3  Axi-symmetrical  2-D  HYDRUS  simulations  of  propellant  in  soil 
column 

A  domain  was  set  up  to  characterize  dissolution  of  NQ  from  M31  propel¬ 
lant  and  transport  in  Plymouth  soil.  The  model  was  run  in  direct  mode  us¬ 
ing  parameters  that  were  determined  for  M31  particles.  The  simulated 
breakthrough  curve  was  similar  to  the  experimental  curve.  Further  work 
will  be  needed  to  develop  this  capability,  but  the  initial  results  indicate  the 
potential  for  using  this  method  to  mechanistically  describe  energetics  dis¬ 
solution  and  transport  in  soil.  The  main  challenge  now  is  to  describe  mul¬ 
tiple  particles  using  axi-symmetrical  3-D. 
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4  Conclusions  and  Implication  for  Future 
Research 

We  collected  single-  and  double-base  propellant  residues  from  live  fire 
training  and  found  that  they  had  energetic  contents  similar  to,  but  ~20% 
lower  than,  their  unfired  parent.  The  appearance  of  the  residue  was  dictat¬ 
ed  by  the  shape  of  the  original  propellant  grain;  un-perforated  grains 
yielded  smaller  versions  of  themselves,  single  perforated  grains  yielded 
rings,  and  multiple  perforated  grains  deposited  slivers.  No  residues  were 
found  when  triple-base  propellants  were  fired.  This  was  also  true  for  col¬ 
lections  made  when  155-mm  mortars  were  fired  using  single-perforation 
Mi  propellant. 

Microscopically,  we  were  able  to  distinguish  between  NC  and  crystalline 
NQ  but  were  not  able  to  map  the  distribution  of  2,  4-DNT  or  NG  in  the  ni¬ 
trocellulose  matrix,  information  important  for  understanding  how  these 
compounds  dissolve  out  from  the  nitrocellulose.  The  different  compounds 
within  propellants  are  difficult  to  image  because  they  are  composed  of 
amorphous  organic  compounds  with  very  similar  compositions.  Also,  be¬ 
cause  they  are  organic,  the  carbon,  nitrogen,  and  oxygen  energy-loss  peaks 
are  weak  and  diffuse  instead  of  sharp,  intense  peaks  seen  in  inorganic 
compounds.  Nevertheless,  our  images  and  data  suggest  that  NG  remains  a 
liquid  within  the  NC  matrix.  We  know  NG  is  added  as  a  liquid  during 
manufacture,  and  others  have  noted  that  NG  migrates  toward  the  surface 
when  present  in  concentrations  >20%.  For  propellants  with  high  fractions 
of  NC  (80-90%),  we  saw  no  change  in  the  propellant  grain  size  or  the  size 
of  the  residues  when  these  were  dripped  on  or  immersed  in  water.  As  the 
percentage  of  NG  increased,  however,  we  noted  that  the  structure  of  the 
propellant  did  change.  This  was  especially  noticeable  for  propellants  with 
<6o%NC  where  the  grain  collapsed  leaving  only  NC. 

The  dissolution  versus  time  curves  for  2, 4-DNT,  NG,  and  NQ  from  unfired 
propellants  show  different  patterns.  Loss  of  2, 4-DNT  is  slower  than  that 
for  NG  or  NQ.  Furthermore,  the  shapes  of  the  mass  loss  versus  time  curves 
are  linear  for  unfired  single-base  propellants  containing  2, 4-DNT,  whereas 
both  the  NG  and  NQ  in  double-  and  triple-base  propellants  show  a  fast  ini¬ 
tial  mass  loss  rate  followed  by  a  slower  dissolution  rate.  We  found  that  the 
amount  of  NG  dissolved  is  related  to  the  NG/NC  ratio.  We  think  droplets 
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of  NG  near  the  surface  are  readily  dissolved,  followed  by  internal  diffusion 
of  NG  through  the  NC  matrix.  In  some  cases,  such  as  in  triple-base  propel¬ 
lants  where  crystals  of  NQ  are  present,  dissolution  of  the  crystals  may 
open  channels  within  the  propellant. 

Acetonitrile  is  an  excellent  solvent  for  extracting  2,4-DNT  and  NG  from 
nitroglycerin  based  propellants.  NQ,  however,  is  less  soluble  in  ACN  and 
our  low  mass  balance  for  NQ  in  our  unfired  triple-based  propellants  high¬ 
lighted  this  problem.  We  suggest  that  a  two-step  extraction  method  might 
work  better.  If  the  propellants  were  crushed  and  placed  in  water,  the  NQ 
would  quickly  dissolve  in  the  water  and  an  aliquot  could  be  used  to  quanti¬ 
fy  the  NQ.  The  crushed  propellant  and  an  aliquot  of  the  water  could  then 
be  extracted  in  ACN  to  quantify  the  NG. 

We  implemented  a  series  of  simple  diffusion  models  to  understand  the  ef¬ 
fects  of  geometry,  concentration  dependence,  and  wet/dry  cycling  on  dis¬ 
solution  rate  from  propellants.  Fitting  the  initial  slopes  of  the  dissolution 
curves  did  not  yield  similar  diffusivities  for  the  propellants,  even  if  concen¬ 
tration  dependent  diffusion,  different  grain  shapes,  and  wet/dry  cycling 
were  taken  into  account.  Fitting  the  tail  slopes  of  the  dissolution  curves, 
however,  gave  the  same  diffusivity  rates  for  different  data  sets  of  the  same 
propellant  and  varied  by  less  than  a  factor  of  four  among  the  propellant 
types  (3.2x10-14  to  1.2x10-13  cm2  s-1)-  We  think  these  values  measure 
the  long-term  diffusion  rates  of  these  materials  and  that  a  different  pro¬ 
cess,  dissolution  of  near  surface  NG  droplets,  is  responsible  for  the  initial 
high  NG  loss.  Once  this  outer  layer  of  droplets  is  depleted,  however,  NG 
would  need  to  diffuse  through  the  solid  NC  matrix  to  reach  the  water,  with 
considerably  lower  diffusivity  (~io-14  cm2  s-1)-  Late  time-dissolution 
would  thus  be  limited  by  molecular  diffusion. 

This  concept  also  qualitatively  accounts  for  the  much  lower  dissolution 
rates  observed  for  fired  grains.  All  dissolution  of  NG  from  fired  grains 
would  thus  be  limited  by  molecular  diffusion.  Indeed,  the  test  data  for 
fired  grains  do  not  show  an  abrupt  roll-off  characteristic  of  unfired  grains 
and  late-time  dissolution  rates  as  similar  to  those  of  unfired  grains.  Our 
microscopy  results  suggest  that  NG  droplets  are  being  directly  dissolved 
from  the  surface,  a  process  that  would  have  to  be  included  in  future  mod¬ 
eling. 
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Soil  column  studies  using  Ml,  WC  860,  and  M31  unfired  propellants  and 
the  fired  residues  of  Ml  and  WC  860  showed  some  of  the  same  dissolution 
trends  as  the  drip  and  batch  tests  where  no  soil  was  present.  For  WC  860 
and  M31,  both  NG  and  NQ  concentrations  in  the  outflow  were  initially 
high  and  decreased  with  time.  The  fired  WC  860  residues  released  a  lower 
percentage  of  their  NG,  supporting  microscopic  observations  that  the  rim 
of  the  residues  had  burned  and  were  depleted  in  NG. 

Different  interactions  between  2,4-DNT,  NG,  NQ  and  the  soils  were  seen 
in  both  the  soil  batch  and  column  studies.  The  2,4-DNT  interacted  strong¬ 
ly  with  soils  and  had  the  highest  adsorption  and  transformation  rates 
measured.  There  was  also  a  strong  correlation  between  soil  adsorption  and 
organic  matter  content  in  the  soil.  As  a  result,  there  was  no  2,4-DNT  in  the 
column  outflow.  This  result  is  consistent  with  lack  of  2,4-DNT  in  ground- 
water  under  firing  points  even  when  it  is  deposited  and  detected  in  soils. 

As  was  found  in  the  drip  tests,  less  NQ  than  NG  was  dissolved  from  triple¬ 
base  propellants,  despite  a  higher  NQ  content  in  the  particles.  We  found 
that  NQ  has  a  low  potential  for  soil  adsorption,  with  adsorption  coeffi¬ 
cients  for  a  variety  of  soils  all  being  less  than  one;  and  no  degradation  or 
transformation.  We  would  expect  that  NQ  dissolved  from  propellant  parti¬ 
cles  would  reach  groundwater.  Fortunately,  NQ  appears  to  be  less  toxic 
than  2,4-DNT  and  NG  (Nipper  et  al.  2009;  Table  15). 


Table  15.  Summary  of  important  characteristics  determined  for  single-,  double-  and  triple-based  propellants. 


2,4-DNT 

NG 

NQ 

Propellant 

Single-base 

Double-  and  triple-base 

Triple-base 

Solubility 

at  25  °C 

270  mg/L 

1500  mg/L 

2600  mg/L 

Phelan  and  Barnett 
(2001) 

Yinon  (1999) 

Haaget  al.  (1990) 

Fate/Behavior 

Strongly  sorbed  to 
organic  matter 

Moves  through  soil  but  subject  to 
rapid  aerobic  biodegradation 

Some  photo  degradation.  Little 
attenuation,  no  biodegradation 

Drinking  water 

0.2  pg/L 

1.5  pg/L 

1.6x103  pg/L 

screening  levels 

EPA-  Region  9 

EPA-  Region  9 

EPA-  Region  9 

Log  KOC 

0.31  (-0.61  to  0.79) 

-0.16  (-.062  to  0.64) 

-0.73  (-0.99  to  -0.30) 

Half  life,  h* 

85  (29  to  408) 

8  (2  to  347) 

79  (15  to  56) 

Toxicity 

Carcinogenic 

EPA-  Region  9, 

Quigley  (1994) 

Non  carcinogen,  EPA-  Region  9 

Affects  the  cardiovascular  system, 
blood,  and  nervous  system, 

Hathaway  et  al.  (1991) 

Non  carcinogen;  non-  endocrine 
disrupter 

EPA-  Region  9,  Hiatt  et  al.  (1988) 

*  Data  summarized  from  Tables  A2  and  A3;  average  and  range  listed. 
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5  Case  Study:  Sampling  for  NG  at  the 
Automatic  Record  Fire  (ARF)  Range, 
Florence,  AZ 

5.1  Introduction 

The  Automatic  Record  Fire  (ARF)  range  (Fig.  49)  at  the  Florence  Military 
Reservation  was  opened  for  training  in  October  of  2010.  The  area  had  not 
been  previously  used  for  military  training  and  hence  the  soils  were  devoid 
of  energetic  compounds  deposited  by  training.  The  soils  at  the  AFR  are  the 
Gunsight-Pinamt  complex  (marked  on  Fig.  50  as  26)  and  the  Laveen  loam 
(marked  on  Fig.  50  as  28  and  60)  both  grouped  as  Aridisols  formed  from 
alluvial  material  (Soil  Survey  Staff  2012).  They  have  gravely  loam  to  sandy 
loam  textures  and  are  well  drained,  with  hydraulic  conductivities  between 
14.4  and  151.2  mm  h_1.  These  soils  have  a  calcic  horizon  (commonly  called 
caliche),  are  alkaline  (pH  of  7.9-84),  low  in  organic  matter  (less  than 
0.5%  throughout),  and  contain  about  14-28%  clay.  Total  rainfall  between 
the  time  the  range  was  opened  and  when  we  sampled  was  116.33  mm  or 
4.58  in.  (Fig.  51). 


Figure  49.  Automatic  Record  Fire  range  at  the  AZ  National  Guard 
Florence  Military  Reservation. 
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Figure  50.  National  Cooperative  Soil  Survey,  Web  Soil  Survey  map 
(Soil  Survey  Staff  2012)  of  the  Florence  Military  Reservation  area 
that  includes  the  Automatic  Record  Fire  range.  Black  lines  show 
outline  of  the  range  itself.  Numbers  indicate  soil  groups:  26— 
Gunsight-Pinamt  complex,  28— Laveen  loam. 


Figure  51.  Rainfall  record  (October  2010  and  September  2011)  for 
Florence  Military  reservation. 


Since  it  opened,  the  ARF  range  has  been  used  to  fire  small  arms,  including 
rifles,  machine  guns,  and  pistols.  The  majority  of  rounds  fired  were  5.  56- 
mm,  but  0.45-  and  0.50-cal.,  and  9-mm  cartridges  were  also  used  (Table 
16).  Because  the  exact  number  of  rounds  fired  at  the  site  is  known,  sam¬ 
pling  the  range  soil  provides  the  first  opportunity  to  relate  range  usage  to 
energetic  compound  concentrations  in  the  soil.  The  propellants  WC844 
(used  in  5.56-mm  rounds)  and  WC  860  (used  in  50-cal.  rounds)  are  dou¬ 
ble-base  and  contain  8  to  11%  nitroglycerin  (NG)  in  an  85%  nitrocellulose 
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matrix  (Defense  Ammunition  Center  2006).  All  of  the  propellants  fired  at 
this  range  contained  NG. 


Table  16.  Record  of  ammunition  fired  on  the  Automatic  Record  Fire  (ARF)  range  between  1  October  2010  and  9 
August  2011. 


Munition 

(Mil/DODIC) 

Weapon 

Propellant 

Constituent 

No.  rounds 
fired 

M855/AA33 

Ball 

M16A2  Automatic  Rifle  (5.56-mm) 

WC844 

NG 

34,198 

M2  or  M33/ A555 

Ball  Linked 

M2  Machine  Gun  (.50-cal.)* 

WC860 

NG 

250 

M1911/A475 

Ball 

Pistol  (.45-cal)** 

SR7970 

NG 

51 

M882/A363 

Ball 

M9  Pistol  (9-mm) 

WPR289 

NG 

23 

M856/A063 

Tracer 

M16A2  Automatic  Rifle  (5.56-mm) 

WC844T 

NG 

1,200 

M855/A062 

Ball  Linked 

M249  Machine  Gun  (5.56-mm) 

WC844 

NG 

2,400 

M855/A059 

Ball 

M16A2  Automatic  Rifle  (5.56-mm) 

WC844 

NG 

108,276 

Total 

146,398 

*  12.7-mm,  **  11.4-mm,  (Defense  Ammunition  Center  2006) 


5.2  Materials  and  methods 

The  ARP  range  is  240  by  300  m  and  has  12  firing  lanes  along  its  short  side 
(Fig.  52).  We  sampled  soils  from  decision  units  that  were  8  m  wide  and 
approximately  10  m  long,  about  2  m  in  front  of  firing  positions  3  through 
12.  The  first  two  firing  lanes  were  not  sampled  as  they  were  difficult  to  ac¬ 
cess,  were  overgrown  with  vegetation,  and  were  not  in  use.  We  collected 
approximately  96  increments,  using  a  2.5-cm-diameter  core  pushed  2  cm 
into  the  soil,  and  combined  these  into  one  multi-increment  sample  of 
about  l-kg  size  for  each  decision  unit.  For  three  of  the  decision  units  (lanes 
3,  5,  and  9),  we  collected  three  multi-increment  samples.  In  addition,  we 
excavated  six  soil  pits  in  two  of  the  decision  units  and  collected  single  pro¬ 
file  samples  to  determine  if  NG  was  present  in  the  subsoil. 

The  sizes  and  locations  of  the  decision  units  were  determined  based  on  the 
residue  distribution  pattern  found  by  firing  tests  of  the  same  munitions 
(Walsh,  M.E.  et  al.  2007).  Walsh  et  al.  found  that  the  majority  of  the  resi¬ 
dues  were  located  in  a  42-m2  area  immediately  in  front  of  the  firing  point. 
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We  placed  our  decision  units  2  m  in  front  of  the  firing  points,  as  these 
were  2  m  above  ground  (Fig.  53),  causing  the  residues  to  travel  a  longer 
horizontal  distance  before  being  deposited.  Background  samples  in  areas 
away  from  the  firing  range  were  also  collected.  The  average  area  of  the  de¬ 
cision  units  was  74  m2. 


Figure  52.  A  Google  Earth  image  showing  the  Automatic  Record  Fire 
(ARF)  range  at  National  Guard  Florence  Military  Reservation,  AZ. 
Firing  positions  are  numbered,  1  through  12. 


Figure  53.  Firing  position  7. 
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To  measure  any  NG  moving  into  the  subsoil  we  excavated  soil  pits  in  two 
of  the  decision  units  (Fig.  54),  three  each  from  firing  lanes  5  and  6.  We  col¬ 
lected  discrete  soil  samples  from  the  surface  0-1  cm  and  at  1-3,  3-5,  5-10, 
10-15, 15-20,  and  20-25  cm  depth  intervals.  An  impermeable  caliche  lay¬ 
er  limited  sampling  depth  (Fig.  54b).  Pits  were  located  about  5  to  6  m  from 
the  firing  point  in  the  middle  of  the  lanes.  These  lanes  were  selected  be¬ 
cause  they  were  close  to  the  observation  tower  and  are  used  preferentially 
(Fielding  2012),  so  would  contain  the  most  propellant  residues  and  a 
greater  potential  for  subsurface  transport. 


Figure  54.  Sampling  a  soil  pit  (a)  and  view  of  pit  showing  a  low  permeability  white  caliche 
layer  (b). 


The  multi-increment  samples  were  processed  and  analyzed  according  to 
standard  procedures  (U.  S.  Environmental  Protection  Agency  2006). 
Briefly,  samples  were  air-dried,  sieved,  and  the  <2-mm  size  fraction 
ground  for  90  s,  homogenized,  and  10-g  subsamples  were  extracted  over¬ 
night  (18  h)  using  acetonitrile  at  1:2,  soihacetonitrile,  ratio.  The  pit  sam¬ 
ples  were  dried  and  sieved  and  the  entire  sample  extracted.  To  make  the 
sample  mass  manageable,  we  subdivided  them  into  three  subsamples,  and 
extracted  each  subsample  using  acetonitrile  in  1:2  soil  to  acetonitrile  ratio, 
similar  to  the  multi-increment  samples.  Acetonitrile  extracts  were  diluted 
1:3  with  water,  filtered  through  a  PTFE  filter,  and  analyzed  for  NG  as  de¬ 
scribed  in  the  Methods  and  Materials  section  (5.2). 
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5.3  Results 

Chemical  analyses  of  the  multi-increment  samples  show  no  NG  in  the 
background  samples  and  surface  soil  NG  concentrations  that  varied  be¬ 
tween  0.31  and  2.92  mg  kg”1.  The  central  lanes,  6  and  7,  had  the  highest 
NG  concentrations  (2.92  and  2.80  mg  kg”1,  respectively,  Table  17.)  con¬ 
sistent  with  the  observation  that  these  lanes  are  used  most  frequently.  Var¬ 
iability  between  three  multi-increment  samples  taken  of  the  same  decision 
unit  was  low,  with  relative  standard  deviation  (RSD)  between  triplicates  of 
2, 12,  and  25%.  The  highest  %RSD  was  for  the  sample  with  the  lowest  con¬ 
centrations.  These  low  RSD  indicate  that  the  multi-increment  samples  ac¬ 
curately  measured  the  average  NG  concentration  present  in  the  decision 
units  sampled. 


Table  17.  Concentration  of  NG  in  multi-increment  surface 
samples  collected  at  the  Automatic  Record  Fire  range. 


Firing  point 

Measured 

mg  kg-1  soil 

Ave  ±  SE 

(n=3) 

RSD 

% 

3 

0.40 

0.31  ±  0.04 

25.3 

0.25 

0.28 

4 

0.85 

5 

1.36 

1.36  ±  0.02 

2.1 

1.33 

1.38 

6 

2.92 

7 

2.80 

8 

0.83 

9 

1.41 

1.43  ±  0.1 

12.3 

1.62 

1.27 

10 

0.52 

11 

1.21 

12 

1.16 

Average 

1.34  ±  0.28 

Background 

B.D. 

B.D. 

SE  =  SD/n2,  standard  error,  RSD  =  relative  standard  deviation,  B.D.  = 
below  detection,  Control  recovery  100%,  blank  =  0  mg  L-1 


The  top  surface  samples  from  the  soil  pits  had  NG  concentrations  that 
generally  agreed  with  the  results  for  the  multi-increment  samples  (Table 
18,  Fig.  55).  These  discrete  samples  had  higher  variability  (29-60%  RSD) 
than  the  multi-increment  samples.  No  NG  was  detected  in  the  subsurface 
samples,  indicating  no  measurable  transport.  We  think  the  lack  of  rain  is 
primarily  responsible.  The  total  rainfall  for  the  month  before  we  sampled 
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was  about  10  mm,  a  quantity  unlikely  to  cause  downward  transport  of  any 
dissolved  NG.  The  maximum  monthly  rainfall  measured,  30  mm,  would 
also  keep  any  dissolved  NG  in  the  surface  layers  of  the  soil.  The  soil  char¬ 
acteristics  may  also  play  a  role  but  the  role  is  unclear.  The  low  organic 
matter  content  suggests  little  attenuation  of  NG  while  the  clay  should  ad¬ 
sorb  and  transform  the  NG.  Both  absorption  and  transformation  of  NG  are 
positively  correlated  with  the  clay  content  (Brannnon  and  Pennington 
2002),  Subsurface  samples  collected  at  other  ranges,  e.g.,  Fort  Richardson, 
Alaska,  contain  high  NG  concentrations  at  the  surface  that  diminished 
with  depth  (Jenkins  et  al.  2008,  chapter  8). 


Table  18.  Concentration  of  NG,  mg  kg-1  soil,  from  the  0-1  cm  layer  in  the 
soil  pits  at  the  Automatic  Record  Fire  range. 


Firing 

point 

Pit 

0-1  cm  soil 

cone. 

Ave.  pit  ±  SE 
(n=3) 

%RSD 

Ave.  firing  point  ±  SE 
(n=3) 

%RSD 

5 

1 

1.20 

1.66±0.32 

32.8 

2.68±0.94 

60 

1.53 

2.27 

2 

1.71 

1.83±0.19 

18.3 

1.57 

2.21 

3 

4.82 

4.55±0.16 

6.0 

4.28 

4.55 

6 

4 

1.92 

2.11±0.32 

26.4 

2.38±0.39 

29 

1.68 

2.74 

5 

2.84 

3.15±0.18 

9.9 

3.14 

3.46 

6 

1.72 

1.87±0.28 

25.8 

1.48 

2.40 
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2.11  ±0.32  3.15  ±0.18  1.87  ±0.28  1.66  ±0.32  1.83  ±0.19  4.55  ±0.16 

0  -  1  cm 
1-  3  cm 

3  -  5  cm 

5  - 10  cm 


10  - 15  cm 


15  -  20  cm 


20  -  25  cm 


L 


mg  kg'1 


ND 


ND 


ND 


ND 


ND 


ND 


mg  kg'1 

-+ 


ND 


ND 


ND 


ND 


ND 


ND 


Pit  4  Pit  5 


mg  kg4 

-+ 


ND 


ND 


ND 


ND 


ND 


ND 


mg  kg"1  mg  kg"1 


— Np — 

- W - 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

Pit  6  Pitl  Pit  2 

_ I  I _ 


mg  kg'1 


ND 


ND 


ND 


ND 


ND 


ND 


Pit  3 


J 


Firing  point  6  Firing  point  5 

Figure  55.  Depth  and  NG  concentration  of  pit  samples  collected  at 
the  Automatic  Record  Fire  range. 
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Appendix  A:  Propellant  and  Soil  Property 
Data 


Table  Al.  Nitrocellulose-based  propellants  used  by  the  U.S.  Army.* 


Used  in 

Propellant 

Constituent 

Percent 

shape 

Single-base 

mortars 

M10 

NC 

98 

flake 

mortars 

M10 

NC 

98 

single  perforated  grain 

mortars 

M10 

NC 

98 

multi-perforated  grain 

tank  training 

M14 

DNT 

8 

multi-perforated  cylinder 

artillery 

Ml 

DNT 

10 

single  or  multi-perforated  cylinder 

artillery 

Ml 

DNT 

10 

single  or  multi-perforated  cylinde 

artillery 

Ml 

DNT 

10 

single  or  multi-perforated  cylinde 

40- to  155-mm  rounds 

M6 

DNT 

10 

multi-perforated  cylinder 

75-mm 

M6+2 

DNT 

10 

multi-perforated  cylinder 

Double-base 

? 

M38 

NG 

2 

grains 

120-mm  mortar 

M45 

NG 

10 

grain  with  single  perforation 

120-mm  mortar 

M47 

NG 

10 

grains 

90-mm  artillery  round 

M5 

NG 

15 

single-perforated  cylinder 

40  mm  grenades 

M2 

NG 

19.5 

single  perforated  grain 

76-  and  90-mm  artillery 

M2 

NG 

19.5 

multi-perforated  cylinder 

? 

M28 

NG 

23 

Cast  grain 

106-mm  round 

M26 

NG 

25 

multi-perforated  cylinder 

152-mm  round 

M26A1 

NG 

25 

multi-perforated  cylinder 

76  mm 

M7 

NG 

32 

single  perforated  stick 

5  in  round 

M7 

NG 

35.5 

single  perforated  stick 

40-,  60-  and  81-mm  rounds 

M9 

NG 

40 

flake 

60-,  80-mmand  4.2"  mortars 

M8 

NG 

43 

sheet 

120-mm  mortar 

M44 

NG 

44 

Single  perforated  stick 

small  arms 

WC844 

NG 

10 

ball  propellant 

small  arms 

WC846 

NG 

10 

ball  propellant 

small  arms 

WC860 

NG 

10 

ball  propellant 

small  arms 

WPR289 

NG 

10 

ball  propellant 
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Used  in 

Propellant 

Constituent 

Percent 

shape 

artillery 

Triple-base 

155  mm 

M30A2 

NQ&NG 

46 

7-perforated  cylinder 

155  mm 

M30A1 

NQ&NG 

47 

7-perforated  cylinder 

90-  155mm 

M30 

NQ&NG 

48 

7-perforated  cylinder 

8” 

M31A1 

NQ&NG 

54 

7-perforated  cylinder 

76-,90-,  120-mm 

M15 

NQ&NG 

55 

single  or  multi-perforated  cylinder 

120, 155-mm 

M31 

NQ&NG 

55 

single  or  multi-perforated  cylinder 

155  mm 

M31A1E1 

NQ&NG 

55 

multi-perforated  cylinder 

Rocket  motor 

M7 

NQ&NG 

55 

Single  perforated  stick 

*Data  taken  from  Propellants:  A  guide  to  recognizing  and  identifying  specific  types  of  nitrocellulose-based 
propellant-online  resource.  Contact  information  for  manual: 

mcal.dac.simac-avpersonnel@conus.arnnv.nnil 

U.S.  Army  Defense  Ammunition  Center  (DAC),  McAlester,  OK 


Table  A2.  Propellant  constituent  adsorption  coefficients  (Kd)  and  associated  soil  property  data. 


Soil  Type 

Kd 

(cm3  g-1) 

CEO 

cmol+ 

(kg-1) 

TOC 

(%) 

Clay 

(%) 

Reference 

2,4-DNT 

Sharkey/Yokena  clay 

Surface 

12.5 

38.9 

2.4 

48.7 

Pennington  et  al.,  2001 

Plymouth 

Surface 

0.28-2.01 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2009 

Plymouth 

Surface 

3.3 

9.2 

1.84 

NM 

Yamamoto  et  al.,  2004 

Yokena/Sharkey  clay 

Surface 

9.43 

38.9 

2.4 

49 

Pennington  etal.,  2003 

Picatinny 

Surface 

2.06 

9.8 

0.63 

5 

Pennington  etal.,  2003 

Grange  Hall 

Surface 

0.43 

16.7 

0.29 

10 

Pennington  etal.,  2003 

Cotto  clay 

Surface 

1.46 

25.7 

0.67 

NM 

Hernandez  et  al.,  2006 

Isabella  sand 

Surface 

0.0172 

2.21 

0.07 

NM 

Hernandez  et  al.,  2006 

Halifax  marine  sed. 

Sediment 

11.7 

14.5 

1.9 

8.25 

Yang  etal.,  2008 

K+-kaolinite 

Pure  clay 

690 

0.3 

NM 

NM 

Haderlein  et  al.,  1996 

K+-illite 

Pure  clay 

3650 

1.6 

NM 

NM 

Haderlein  et  al.,  1996 

K+-montmorillonite 

Pure  clay 

7400 

12 

NM 

NM 

Haderlein  et  al.,  1996 

2,6-DNT 

Sharkey  clay 

Surface 

5.96 

38.9 

2.4 

48.7 

Pennington  et  al.,  2001 

Plymouth 

Surface 

0.84-0.97 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2009 

Yokena/Sharkey  clay 

Surface 

3.56 

38.9 

2.4 

49 

Pennington  etal.,  2003 

Picatinny 

Surface 

1.61 

9.8 

0.63 

5 

Pennington  etal.,  2003 

Grange  Hall 

Surface 

0.07 

16.7 

0.29 

10 

Pennington  etal.,  2003 

Halifax  marine 

Sediment 

8.3 

14.5 

1.9 

8.25 

Yang  etal.,  2008 

K+-kaolinite 

Pure  clay 

10 

0.3 

NM 

NM 

Haderlein  et  al.,  1996 
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Soil  Type 

Kd 

(cm3  g-i) 

CEC 

cmol+ 

(kg-1) 

TOC 

(%) 

Clay 

(%) 

Reference 

K+-illite 

Pure  clay 

52 

1.6 

NM 

NM 

Haderlein  et  al.,  1996 

K+-montmorillonite 

Pure  clay 

125 

12 

NM 

NM 

Haderlein  et  al.,  1996 

NG 

Plymouth 

Surface 

0.17 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2007 

Adler 

Surface 

0.08 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2007 

Plymouth 

Surface 

0.165 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2008 

Adler 

Surface 

0.113 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2008 

Plymouth 

Surface 

2.8 

9.2 

1.84 

NM 

Speitel  et  al.,  2002 

Picatinny 

Surface 

3.8 

9.8 

0.634 

5 

Pennington  etal.,  2002 

NQ 

Plymouth 

Surface 

0.14 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2007 

Adler 

Surface 

0.03 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2007 

Plymouth 

Surface 

0.069 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2008 

Adler 

Surface 

0.034 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2008 

soil 

<0.1 

Haag  etal.,  1990 

Grange  Hall 

Surface 

0.15 

16.7 

0.3 

10 

Pennington  etal.,  2004 

Yokena  clay 

Surface 

0.43 

38.9 

2.4 

48.75 

Pennington  etal.,  2004 

Table  A3.  Soil  transformation  rate  coefficients  (k)  for  energetic  compounds  under  different  soil  conditions. 


Soil 

Oxid 

k  (hr*) 

Half  life  (h) 

CEC  (umol/g) 

%TOC 

%Clav 

Reference 

2,4-DNT 

LAAPML 

An 

0.017 

40 

3.5 

0.015 

5 

Pennington  etal.,  1999 

LAAP  SP-SM 

An 

0.0017 

410 

3.6 

0.015 

5 

Pennington  et  al.,  1999 

LAAP  CL 

An 

0.0021 

330 

8.1 

0.162 

15 

Pennington  etal.,  1999 

LAAP  SM 

An 

0.0021 

320 

5.5 

0.02 

7.5 

Pennington  et  al.,  1999 

Yokena  clay 

An 

Rapid 

38.9 

2.4 

48.7 

Pennington  et  al.,  2001 

Plymouth 

Aero 

0.002-0.023 

8 

347-29 

4.4 

0.49 

0.8 

Dontsova  etal.,  2009 

2,6-DNT 

LAAP  C 

An 

0.0023 

301 

6.6 

0.08 

12 

Pennington  et  al.,  2001 

LAAP  D 

An 

0.0035 

198 

15.5 

0.2 

32 

Pennington  et  al.,  2001 

Yokena  clay 

An 

0.0235 

29 

38.9 

2.4 

48.7 

Pennington  et  al.,  2001 

Plymouth 

0.0085 

82 

4.4 

0.49* 

0.8 

Dontsova  etal.,  2009 

2,4-and  2,6-DNT 

Groundwater 

Aero 

0.00011 

4320-672 

NM 

NM 

NM 

Howard  et  al.,  1991 

Groundwater 

An 

0.01441 

300-48 

NM 

NM 

NM 

Howard  et  at,  1991 

Soil 

Aero 

0.2-1.51E-3 

4320-672 

NM 

NM 

NM 

Howard  et  al.,  1991 

NG 
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Soil 

Oxid 

k  (hr1) 

Half  life  (h) 

CEC  (umol/g) 

%TOC 

%Clav 

Reference 

LAAP-D/SM 

0.002 

335 

5.5 

0.02 

7.5 

Pennington  et  al.f  2002 

Picatinny 

0.008 

84 

9.8 

0.634 

5 

Pennington  et  al.f  2002 

Yokena  clay 

0.095 

7.3 

38.9 

2.4 

48.75 

Pennington  et  al.,  2002 

Plymouth2 

0.008 

88.1 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2007 

Adler2 

NS 

NS 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2007 

Yakima  TC 

>0.029 

<24 

19.6 

2 

20.1 

Jenkins  et  al.,  2003 

Camp  Guernsey 

>0.029 

<24 

30.3 

1.6 

25.0 

Jenkins  et  al.,  2003 

Fort  Greely 

0.059 

11.8 

9.0 

1.1 

3 

Jenkins  et  al.,  2003 

LAAP-D 

Aero 

0.020 

38.3 

15.5 

0.2 

32 

Yost,  2004 

Yokena  clay 

Aero 

0.018 

38.1 

38.9 

2.4 

34 

Yost,  2004 

LAAP-D 

An 

0.310 

3.2 

15.5 

0.2 

32 

Yost,  2004 

Yokena  clay 

An 

0.263 

2.7 

38.9 

2.4 

34 

Yost,  2004 

NQ 

Yakima 

0.0449 

15.4 

23.6 

1.06 

15 

Mulherin  et  al.,  2005 

Lebanon  Landfill 

0.0123 

56.4 

NM 

0.3 

50 

Mulherin  et  al.,  2005 

Fort  Edwards  Clay 

0.0128 

54.2 

NM 

0.6 

14 

Mulherin  et  al.,  2005 

Plymouth 

NS 

NS 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2007 

Adler 

NS 

NS 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2007 

Grange  Hall 

NS 

NS 

16.7 

0.3 

10 

Pennington  et  al.,  2004 

LAAP  D  (SM?) 

NS 

NS 

5.5 

0.02 

7.5 

Pennington  et  al.,  2004 

Yokena  clay 

NS 

NS 

38.9 

2.4 

48.75 

Pennington  et  al.,  2004 

DPA 

Adler 

0.011 

62.3 

4.4 

0.49 

0.8 

Dontsova  et  al.,  2007 

Plymouth 

0.014 

48.1 

16.6 

0.29 

4.2 

Dontsova  et  al.,  2007 

Grange  Hall 

0.0027 

256 

16.7 

0.3 

10 

Pennington  etal.,  2004 

LAAP  D/SM 

NS 

NS 

5.5 

0.02 

7.5 

Pennington  et  al.,  2004 

Yokena  clay 

1.098 

0.63 

38.9 

2.4 

48.75 

Pennington  et  al.,  2004 

EC 

Grange  Hall 

NS 

NS 

16.7 

0.3 

10 

Pennington  etal.,  2004 

LAAP  D/SM 

NS 

NS 

5.5 

0.02 

7.5 

Pennington  et  al.,  2004 

Yokena  clay 

NS 

NS 

38.9 

2.4 

48.75 

Pennington  et  al.,  2004 

1.  Estimated;  2.  C-14  data;  Oxid=Oxidation  Status;  An=anaerobic;  Areo=aerobic;  NS  =  not  significantly  different  from 
zero;  NM  =  not  measured. 
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Appendix  B:  Soil  Properties  and  Soil  Column 
Test  Data 


Table  Bl.  Properties  of  soils  used  in  column  studies. 


Soil 

Sassafras 

Plymouth 

Catlin 

Property 

aver. 

stdev 

aver. 

stdev 

aver. 

stdev 

pH 

4.40 

± 

0.01 

4.23 

± 

0.01 

6.87 

± 

0.01 

EC,  pS  cm-1 

212 

± 

27 

130 

± 

6 

538 

± 

20 

CEC,  cmolc  kg1 

16.58 

± 

1.32 

18.78 

± 

2.87 

23.59 

± 

0.98 

OC,  % 

1.30 

0.02 

1.72 

0.05 

3.75 

0.15 

Clay,  % 

16.37 

± 

0.07 

14.40 

± 

0.74 

15.73 

± 

0.20 

Silt,  % 

42.27 

± 

0.15 

22.15 

± 

0.31 

71.53 

± 

0.06 

Sand,  % 

41.37 

± 

0.09 

63.46 

± 

1.05 

12.75 

± 

0.26 

Texture 

loam 

sandy  loam 

silt  loam 

pH  (1:1  soil:water);  EC  =  electrical  conductivity  in  1:1  soil: water,  (pS  crrr1);  CEC  =  cation  exchange  capacity  by  so¬ 
dium  acetate  method,  cmolc  kg1;  OC  =  total  organic  carbon  by  combustion,  %;  %  Clay  (<2  pm),  Silt  (2  -  50  pm),  and 
Sand  (50  -  2000  pm)  by  laser  diffraction. 


Table  B2.  NG  loss  from  unfired  and  fired  WC860  propellant  and  unfired  M31  propellant  in  column  experiments  conducted  at  0.01  and 
0.02  ml_  min-1  flow  measured  by  extraction  of  propellants  and  propellant  residues  after  experiment  completion. 


Soil 

U/F 

CF/ 

IF 

Set 

Col. 

HPLC 

Measured 

cone. 

(pgL-iC 

Dilut. 

Factor 

Total 

NG  in 

extract 

(pgL1) 

ACN 

(mL) 

Total 
ppb  in 
ACN 

Propellant 

weight 

(g) 

Meas. 

NG 

cone. 

in 

prop. 

(mg 

g-1) 

Original 

NG 

cone, 
in  prop. 

(%) 

Final 

meas. 

NG  left 

(%) 

NG 

lost 

(%) 

0.01  mL  min-1  WC860 

Sassafras 

u 

IF 

III 

cl 

1465 

200 

292900 

20 

5858 

0.0860 

68.12 

9.7 

6.81 

29.78 

Sassafras 

u 

CF 

III 

c2 

1883 

200 

376633 

20 

7533 

0.1009 

74.65 

9.7 

7.47 

23.04 

Plymouth 

u 

IF 

1 

c3 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

Plymouth 

u 

CF 

1 

c4 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

Sassafras 

F 

IF 

III 

c3 

2018 

20 

40360 

20 

807 

0.0149 

54.17 

7.8 

5.42 

30.55 

Sassafras 

F 

CF 

III 

c4 

2438 

20 

48753 

20 

975 

0.0167 

58.39 

7.8 

5.84 

25.14 

Plymouth 

F 

IF 

1 

c5 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

Plymouth 

F 

CF 

1 

c6 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

NE 

0.02  mL  min-1  WC860 

Sassafras 

U 

IF 

IV 

cl 

1928 

200 

385667 

20 

7713 

0.0970 

79.52 

9.7 

7.95 

18.02 

Sassafras 

U 

CF 

IV 

c2 

1931 

200 

386100 

20 

7722 

0.0960 

80.44 

9.7 

8.04 

17.07 
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Soil 

U/F 

CF/ 

IF 

Set 

Col. 

HPLC 

Measured 

cone. 

(UgL-iC 

Dilut. 

Factor 

Total 

NG  in 

extract 

(UgL1) 

ACN 

(mL) 

Total 
ppb  in 
ACN 

Propellant 

weight 

<g) 

Meas. 

NG 

cone. 

in 

prop. 

(mg 

g-1) 

Original 

NG 

cone, 
in  prop. 

(%) 

Final 

meas. 

NG  left 

(%) 

NG 

lost 

(%) 

Plymouth 

u 

IF 

II 

cl 

1344 

200 

268700 

20 

5374 

0.0820 

65.54 

9.7 

6.55 

32.44 

Plymouth 

u 

CF 

II 

c2 

1777 

200 

355333 

20 

7107 

0.0895 

79.40 

9.7 

7.94 

18.14 

Average 

21.42 

SE 

3.68 

Sassafras 

F 

CF 

IV 

c3 

2234 

20 

44677 

20 

894 

0.0145 

61.62 

7.8 

6.16 

21.00 

Plymouth 

F 

CF 

II 

c3 

1418 

40 

56720 

20 

1134 

0.0140 

81.03 

7.8 

8.10 

-3.88 

Plymouth 

F 

IF 

II 

c4 

1984 

20 

39680 

20 

794 

0.0160 

49.60 

7.8 

4.96 

36.41 

Average 

17.84 

SE 

10.17 

0.01  mL  min-i  M31 

Sassafras 

U 

IF 

III 

c5 

1159 

200 

231800 

20 

4636 

0.0765 

60.60 

19 

6.06 

68.10 

Sassafras 

U 

IF 

III 

c6 

1453 

200 

290600 

20 

5812 

0.0748 

77.70 

19 

7.77 

59.10 

Plymouth 

u 

IF 

III 

c7 

1748 

200 

349533 

20 

6991 

0.0767 

91.14 

19 

9.11 

52.03 

Plymouth 

u 

CF 

III 

c8 

1308 

200 

261633 

20 

5233 

0.0784 

66.74 

19 

6.67 

64.87 

Average 

61.03 

SE 

3.53 

0.02  mL  min-1  M31 

Sassafras 

u 

CF 

IV 

c5 

3831 

100 

383067 

20 

7661 

0.0810 

94.58 

19 

9.46 

50.22 

Sassafras 

u 

CF 

IV 

c6 

1770 

200 

354000 

20 

7080 

0.0808 

87.62 

19 

8.76 

53.88 

Plymouth 

u 

IF 

IV 

c7 

1176 

200 

235267 

20 

4705 

0.0420 

112.03 

19 

11.20 

41.04 

Plymouth 

u 

CF 

IV 

c8 

1876 

200 

375267 

20 

7505 

0.0796 

94.29 

19 

9.43 

50.37 

Average 

48.88 

SE 

2.75 

ACN=acetonitrile;  SE  =  Standard  error;  CF=continuous  flow;  IF=interrupted  flow;  PV  =  pore  volume;  NE=  not  extracted;  *=Accountingfor 
transformation.  NG  transformation  rate  (k)  used  to  calculate  this  value  was  determined  in  batch  experiments  and  equaled  to  0.0193  h^in  Sassafras 
soil  and  0.0117  h1  in  Plymouth  soil  (Table  7). 
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Table  B3.  Soil  concentrations  of  NG  in  column  tests  using  WC860  and  M31  propellants  and 
soil  concentrations  of  2,4-DNT  in  column  tests  using  Ml  propellant.  For  sets  IV  and  V,  the 
column  was  divided  into  top,  middle  and  bottom  sections  and  each  of  these  extracted 
separately. 


Propellant 

Soil 

U/F 

CF/IF 

Set 

Col. 

Position 

Soil  cone,  mg 
kg1 

WC860 

Sassafras 

u 

IF 

III 

cl 

3.47 

WC860 

Sassafras 

u 

CF 

III 

c2 

3.53 

WC860 

Sassafras 

F 

IF 

III 

c3 

6.14 

WC860 

Sassafras 

F 

CF 

III 

c4 

2.86 

WC860 

Sassafras 

U 

IF 

IV 

cl 

T,M,B 

ND 

WC860 

Sassafras 

u 

CF 

IV 

c2 

TfM,B 

ND 

WC860 

Sassafras 

F 

CF 

IV 

c3 

T,M,B 

ND 

WC860 

Plymouth 

U 

IF 

II 

cl 

ND 

WC860 

Plymouth 

u 

CF 

II 

c2 

ND 

WC860 

Plymouth 

F 

CF 

II 

c3 

ND 

WC860 

Plymouth 

F 

IF 

II 

c4 

2.76 

M31 

Sassafras 

U 

IF 

III 

c5 

ND 

M31 

Sassafras 

U 

IF 

III 

c6 

ND 

M31 

Sassafras 

u 

CF 

IV 

c5 

TfM,B 

ND 

M31 

Sassafras 

u 

CF 

IV 

c6 

T,M,B 

ND 

M31 

Plymouth 

u 

IF 

III 

c7 

2.15 

M31 

Plymouth 

u 

CF 

III 

c8 

2.42 

M31 

Plymouth 

u 

IF 

IV 

c7 

Top 

1.11 

Middle 

0.94 

Bottom 

1.75 

M31 

Plymouth 

u 

CF 

IV 

c8 

Top 

2.29 

Middle 

1.17 

Bottom 

0.84 

Ml 

Sassafras 

u 

IF 

V 

cl 

Top 

0.34 

M,B 

ND 

Ml 

Sassafras 

u 

CF 

V 

c2 

TfM,B 

ND 

Ml 

Sassafras 

F 

CF 

V 

c3 

T,M,B 

ND 

Ml 

Sassafras 

F 

IF 

V 

c4 

T,M,B 

ND 

Ml 

Plymouth 

u 

IF 

V 

c5 

Top 

1.07 

Middle 

0.61 

Bottom 

0.50 

Ml 

Plymouth 

u 

CF 

V 

c6 

Top 

1.59 

Middle 

0.58 

Bottom 

ND 

Ml 

Plymouth 

F 

IF 

V 

c7 

T,M,B 

ND 

Ml 

Plymouth 

F 

CF 

V 

c8 

Top 

1.21 

M,B 

ND 

U=Unfired,  F=  Fired;  CF=continuous  flow;  IF=interrupted  flow;  ND=not  detected 
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Appendix  C:  Breakthrough  Curves  for  Soil 
Column  Tests 


Plymouth,  WC860  unfired,  CF,  0.01  ml.  min1 
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Plymouth,  WC860  unfired,  IF,  0.01  mL  min-1 
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Figure  Cl.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.01  mL  min-1  rate  onto  unfired  0.50-cal.  propellants  in  Plymouth  soil.  The  solid 
vertical  line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Plymouth,  WC860  fired,  IF,  0.01  ml.  min  1 
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Plymouth,  WC860  fired,  IF,  0.01  mL  min  1 


Figure  C2.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.01  mL  min-1  rate  onto  fired  0.50-cal.  propellants  in  Plymouth  soil.  The  solid  vertical 
line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Plymouth,  50  cal  unfired,  IF,  0.02  mL  min 1 
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Figure  C3.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at  0.02  mL 
min-1  rate  onto  unfired  and  fired  0.50-cal  propellants  in  Plymouth  soil.  The  solid  vertical  line 
at  ~30  pore  volumes  shows  when  the  propellant  was  removed  from  the  soil  surface. 
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Plymouth,  50  cal  fired,  CF,  0.02  mL  min1 


Plymouth,  50  cal  fired,  IF,  0.02  mL  min1 


Figure  C4.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.02  mL  min-1  rate  onto  unfired  and  fired  0.50-cal.  propellants  in  Plymouth  soil.  The 
solid  vertical  line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 


ERDC/CRREL  TR-12-9 


105 


Sassafras,  50  cal  unfired,  CF,  0.01  mL  min  1 


Sassafras,  50  cal  unfired,  IF,  0.01  mL  min  1 


Figure  C5.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at  0.01  mL 
min-1  rate  onto  unfired  and  fired  0.50-cal.  propellants  in  Sassafras  soil.  The  solid  vertical  line 
indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Sassafras,  50  cal  fired,  CF,  0.01  mL  min  1 


Sassafras,  50  cal  fired,  IF,  0.01  mL  min1 


Figure  C6.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.01  mL  min-1  rate  onto  unfired  and  fired  0.50-cal.  propellants  in  Sassafras  soil.  The 
solid  vertical  line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Sassafras,  50cal  unfired  IF  0.02  mL  min1 
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Sassafras,  50cal  unfired  CF  0.02  mL  min1 


Figure  C7.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.02  mL  min-1  rate  onto  unfired  and  fired  0.50-cal.  propellants  in  Sassafras  soil.  The 
solid  vertical  line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Sassafras,  50cal  fired  CF  0.02  mL  min1 


Figure  C8.  Breakthrough  curves  for  continuous  water  flow  at  0.02  mL  min-1  rate 
onto  fired  0.50-cal.  propellants  in  Sassafras  soil.  The  solid  vertical  line  indicates 
when  the  propellant  was  removed  from  the  soil  surface. 
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Plymouth,  M31  unfired,  CF,  0.01  mL  min  1 
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Figure  C9.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow  at 
0.01  mL  min-1  rate  onto  unfired  M31  propellant  in  Plymouth  soil.  The  solid  vertical 
line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Sassafras,  M31  unfired,  IF,  0.01  mL  min1 
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Figure  CIO.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow 
at  0.01  mL  min-1  rate  onto  unfired  M31  propellant  in  Sassafras  soil.  The  solid 
vertical  line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Sassafras,  M31  unfired  CF  0.02  mL  min1 
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Figure  Cll.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow 
at  0.02  mL  min-1  rate  onto  unfired  M31  propellant  in  Sassafras  soil.  The  solid 
vertical  line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Plymouth,  M31  unfired  IF  0.02  mL  min1 


1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 


o 


U 

U* 

i- 

CQ 


Plymouth,  M31  unfired  CF  0.02  mL  min1 


o 


U 

u“ 

GO 


Figure  C12.  Breakthrough  curves  for  continuous  (CF)  and  interrupted  (IF)  water  flow 
at  0.02  mL  min-1  rate  onto  unfired  M31  propellant  in  Plymouth  soil.  The  solid  vertical 
line  indicates  when  the  propellant  was  removed  from  the  soil  surface. 
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Appendix  D:  Raman  Spectroscopy  of 
Propellants 

We  tested  the  applicability  of  Raman  spectroscopy  and  confocal  Raman 
microscopy  for  distinguishing  between  nitroglycerin  (NG),  nitroguanidine 
(NQ)  and  nitrocellulose  (NC)  phases  in  double  and  triple-based  propel¬ 
lants.  In  Raman  spectroscopy  a  laser  is  shone  onto  a  sample  surface  and  a 
detector  records  the  sample’s  emitted  radiation  at  different  wavelengths, 
generating  a  Raman  spectrum  of  the  point  with  wavelength  numbers  be¬ 
tween  o  and  4000  cm"1.  As  NG,  NC  and  NQ  all  have  distinct  spectra  (Fig. 
Di),  wavelengths  at  which  prominent  peaks  in  the  spectra  of  a  compound 
occur,  can  be  used  to  map  the  location  of  that  compound  (Confocal  Raman 
Microscopy  General  Overview  2012;  Trewartha  et  al.  2007).  Confocal  Ra¬ 
man  microscopy  merges  Raman  spectroscopy  techniques  with  a  confocal 
microscope.  The  confocal  microscope  focuses  the  laser  beam  on  a  very 
small  area  of  the  sample’s  surface,  thereby  increasing  the  spatial  resolu¬ 
tion  of  the  Raman  spectra. 


Wavelength  [cm]] 


Figure  Dl.  Raman  spectra  of  nitroguanadine,  nitroglycerine  and  nitrocellulose. 
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For  the  Raman  analyses  it  helps  to  have  flat  samples.  We  potted  a  single¬ 
base  (Mio  with  only  NC),  double-base  (M9,  0.50-cal.)  and  triple-base 
(M31)  propellant  in  both  wax  and  epoxy  and  sectioned  their  surfaces  to 
produce  a  flat  working  surface.  Both  of  these  potting  materials  are  less 
than  ideal;  the  epoxy  interacts  with  constituents  in  the  propellants  causing 
changes  to  their  compositions  (Fig.  D2);  the  wax  is  soft  and  the  propel¬ 
lants  occasionally  move  while  being  sectioned.  We  used  the  wax,  potted 
butts  to  obtain  point  and  line  spectra  from  these  samples.  For  the  confocal 
microscope  Raman  analyses,  the  areas  were  small  enough  that  we  did  not 
have  to  section  the  sample.  Samples  were  observed  using  a  20  x  objective 
lens  and  excited  by  a  laser  at  532  nm.  Raman  spectra  were  obtained  from  a 
series  of  adjacent  points  on  the  samples,  allowing  us  to  map  the  different 
chemical  phases. 


Figure  D2.  M31  propellant  in  epoxy;  note  bubbles  formed  by  the 
interaction  of  the  propellant  with  the  epoxy. 


Figure  D3  is  a  line  scan  across  a  0.50-cal.  propellant.  There  are  three 
peaks— a  doublet  centered  at  1300  cm  1  and  a  broad  peak  at  “850  cm1. 
None  of  these  peaks  are  clearly  associated  with  NG,  although  the  doublet 
might  be  a  combination  of  NG  and  NC.  Nitroguanadine,  on  the  other 
hand,  could  be  identified  in  point  spectra  (Fig.  D4)  and  it  could  be  mapped 
using  the  confocal  Raman  (Fig.  D5).  In  Figure  D5,  NQ  crystals  (bright  yel¬ 
low)  are  different  from  the  NC/NG  matrix.  These  results  show  that  5 -pm 
crystals  of  NQ  are  visible  using  Raman  but  that  the  more  intimately  mixed 
NC  and  NG  cannot  be  separated  from  one  another  using  either  regular  or 
confocal  microscope  Raman. 
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Figure  D3.  Line  scan  spectrum  across  a  section  of  a  0.50-cal.  double-base  propellant. 
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Figure  D4.  Raman  point  measurement  of  NQ  in  a  triple-base  propellant. 
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Raman  Image 
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30pm 


NC  map 


NQ  map 


Figure  D5.  M31  surface  map;  the  red  circle  on  the  optical  image  (left)  is  the  center  of  the 
NC/NG  and  NQ  maps  shown  on  the  right;  the  +  mark  the  edges  of  the  maps.  Note  that  the  NQ 
crystals  are  easily  seen  on  the  optical  image. 


Raman  and  Confocal  Raman  were  able  to  distinguish  NQ  from  the  back¬ 
ground  NC  and  NG  matrix,  but  were  not  able  to  distinguish  between  NC 
and  NG.  These  techniques  also  have  some  limitations.  First,  spectra  can 
only  be  taken  from  the  surface  of  the  samples,  requiring  sectioning  of 
samples  if  analysis  of  an  entire  particle  is  desired.  Second,  although  Con- 
focal  Raman  can  map  the  surface  at  high  resolution,  mapping  a  large  area 
of  the  surface  is  time-consuming  as  only  a  tiny  portion  of  the  surface  is  in 
focus  at  each  point.  Finally,  the  propellant  samples  fluoresce,  thereby  de¬ 
grading  the  quality  of  the  spectra  and  the  resulting  maps.  If  Raman  mi¬ 
croscopy  is  to  be  used  for  future  investigations,  the  fluorescence  should  be 
removed  using  commercially  available  background  subtraction  software 
(Cao  et  al.  2007)  and  more  work  is  needed  to  develop  better  potting  and 
sectioning  methods.  For  this  project  we  decided  to  use  more  established 
procedures  for  propellant  analysis. 
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14.  ABSTRACT 

During  firing,  propellant  residues  are  scattered  onto  the  soil  surface  where  their  energetic  compounds  can  be  dissolved  by  precipitation. 
The  residues,  like  the  unfired  propellants,  are  composed  of  nitrocellulose  imbibed  with  either  2,4-DNT  (single-base),  nitroglycerin  (NG) 
(double-base)  or  NG  and  nitroguanidine  (NQ)  (triple-base).  Although  nitrocellulose  is  insoluble,  2,4-DNT,  NG,  and  NQ  are  soluble;  and 
2,4-DNT  and  NG  are  also  toxic.  Consequently,  data  on  how  quickly  2,4-DNT,  NG,  and  NQ  are  dissolved  from  propellant  residues  are 
needed  to  determine  the  flux  of  these  compounds  to  soil.  Once  in  soil  solution,  the  partition  coefficient,  Kd,  and  degradation  rate,  k 
values  are  needed  to  predict  the  transport  of  energetics  through  the  vadose  zone  and  to  groundwater.  We  measured  the  2,4-DNT,  NG, 
and  NQ  dissolution  rates  for  different  propellants  using  laboratory  batch  and  drip  tests  where  no  soil  was  present  and  soil  column 
studies,  which  used  similar  propellant  and  residues  as  source  terms,  to  determine  partition  coefficients  and  degradation  rates.  Because 
the  surfaces  of  propellants  and  residues  may  play  an  important  role  in  dissolution  of  the  energetic  constituents,  we  studied  these  using 
both  light  and  electron  microscopy.  We  found  that  2,4-DNT  is  well  bound  to  NC  and  dissolves  out  slowly,  but  that  both  NG  and  NQ 
have  fast  initial  dissolution  followed  by  slower  mass  loss.  The  amount  of  NG  dissolved  is  a  function  of  the  NG/NC  ratio  in  the 
propellant  and  both  our  mass  loss  data  and  our  microscopy  results  suggest  that  NG  exists  as  fine  liquid  droplets  within  an  NC  matrix 
rather  than  as  dispersed  molecules. _ 
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14.  ABSTRACT  (cont’d) 


NG  droplets  near  the  grain  surface  are  quickly  dissolved  and  once  this  layer  of  liquid  NG  is  depleted,  NG 
diffuses  through  the  NC  matrix  slowly  (~10-14  cm2  s_1).  The  NQ  also  dissolves  rapidly  initially  but 
quickly  mass  loss  for  the  NQ  becomes  smaller  than  that  for  NG,  despite  higher  NQ  concentrations  in  the 
studied  triple-base  propellants.  NQ  is  added  as  a  crystal  during  manufacturing,  and  was  observed  to 
remain  solid  in  the  propellant,  so  dissolution  of  the  NQ  crystal  would  have  to  precede  its  removal  by 
water.  Both  2,4-DNT  and  NG  are  added  as  liquids  and  cannot  be  distinguished  from  the  NC  matrix. 
Therefore,  their  distribution  and  movement  within  the  nitrocellulose  matrix  is  poorly  understood, 
hampering  our  ability  to  derive  a  physically  based  dissolution  model  that  can  predict  energetic  losses 
from  a  variety  of  propellant  types.  Different  interactions  between  2,4-DNT,  NG,  NQ,  and  the  soils  were 
seen  in  both  the  soil  batch  and  column  studies.  The  2,4-DNT  interacted  strongly  with  soils  and  had  the 
highest  adsorption  and  transformation  rates  measured.  As  a  result,  no  2,-4-DNT  was  detected  in  column 
outflow.  NG  experienced  both  adsorption  and  transformation  in  the  soils,  resulting  in  retardation  of  the 
breakthrough  curve  and  decreased  concentrations  in  the  outflow.  The  short  half-life  of  NG  in  most  soils 
suggests  that  it  should  rarely  reach  groundwater.  NQ,  on  the  other  hand,  does  not  readily  adsorb  to  soil, 
and  does  not  degrade  or  transform.  We  would  expect  that  NQ  dissolved  from  propellants  would  reach 
groundwater. 


